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Foreword 


These Lecture Notes introduce the theoretical basics of solid mechanics 
to environmental engineering students. Born out of and supported by the 
European Project DEREC TEMPUS JEP Development of Environmen- 
tal and Resources Engineering Curriculum, it collects the lectures held by 
the Authors during the course of Mechanic of Solids at the University of 
Florence, Degree of Environmental Engineering and Resources. Although 
the course is extended to basic structural engineering principles, such as 
mechanics, statics, kinematics and fundamental equations of beam struc- 
tures, inertia, iso static and hyper static solution methods, these Lecture 
Notes reflect only the content of the lectures of continuum mechanics. 


Several approaches are possible to the subject depending on the con- 
cern, either mathematically or physically oriented. The volume aims to 
provide a synthesis of both approaches, presenting in an organic whole 
the classical theory of solid mechanics and a more direct engineering ap- 
proach. It is the Authors’ opinion that a top-down learning process may 
offer to the engineering students those critical and autonomy tools neces- 
sary to gain awareness of that continuous learning process that is required; 
it characterizes the cultural and technical personality of an engineer. An 
ongoing learning is all the more necessary today, where the rapid develop- 
ment of powerful computers and computer solving methods (finite element 
methods, discrete volume methods, boundary methods, etc.) have opened 
up the way to new horizons that the classical approaches were only able to 
formulate. This fast and impressive growth of computer methods seems to 
be replacing the importance of gaining a consolidated knowledge of solid 
mechanics background. On the contrary, the Authors believe that only 
a conscious knowledge of theory can be that cultural instrument through 
which an engineer can really hope to control the use of computer methods. 
With this aim, the Reader addressed by this volume is mainly the under- 
graduate student in Engineering Schools: it is organized in eight Chapters: 
Chapter 1 proposes a synthesis of the basic concepts of mathematics and 
geometry that the readers need in the following chapters. Chapter 2 and 
Chapter 3 are devoted to the elementary framework of strain and stress 
in an elastic body. The concept of finite strain and Cauchy stress state 
is introduced, together with Mohr’s representation of a general state of 
stress. Chapter 4 focuses on the classical law of linear elasticity. Chapter 
5 deals with the Principle of Virtual Works. Chapter 6 treats the energy 
principles and provides a basic introduction to the variational methods. 
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XIV FOREWORD 


Finally, Part I ends with a chapter introducing the notion of strength of 
materials. At the end of each chapter of the first part the basics of the 
tensor—based shell theory are also presented and then an application to 
some standard shell geometries is provided in appendix A. 

The second part, Chapter 8, is dedicated to De Saint-Venant’s problem 
where the classical beam theory is presented focusing on the four funda- 
mental cases: beam under axial forces, terminal couples, torsion, bending 
and shear. 

The volume, that consolidates the Lecture Notes prepared by the Au- 
thors for the second-year undergraduate students in environmental engi- 
neering, proposes a widening of the classical theories approached, giving 
a list of references used during its preparation as a possible suggestion to 


the Reader. 


The Authors wish to express their heartfelt gratitude to professor 
Marco Modugno for the inspiring discussions and stimulating suggestions. 


It is also our pleasure to thank Eng. Seymour Milton John Spence for 
kindly revising the English text. 


The publication of this book has been possible thanks to the finan- 
cial support of the European Commission (DEREC Tempus Project) and 
Ente Cassa di Risparmio di Firenze to whom the Authors are extremely 
grateful. 
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Chapter 1 
Outline of linear algebra 


This chapter briefly presents some preliminary mathematics necessary 
to understand continuum mechanics. To this end the basic concepts of 
linear algebra and tensor analysis will be introduced. At the end of the 
chapter an overview of the theory of surfaces will be exposed in order to 
make the reader familiar with some background required for the mechanics 
of shell continuums, even though the latter is not the key theme of this 
book. 

This introduction is neither exhaustive nor complete; indeed for any 
further insight the reader is warmly recommended to refer to the main 
sources from which this summary has been derived: Modugno, /4] and 


[5]; Sokolnikoff, /1]; Green-Zerna, [3]. 


1.1 Vector spaces and linear mappings 


1.1.1 Vector spaces 


We define vector space a set V equipped with the following op- 
erations 


do. (1.2) 


Elements belonging to V are named vectors and are character- 
ized by the following properties 


1. Uu+(0+0)=(0+0)+% vuvwevV 


4. vaeVi=-u€V so that ùu+(-u)=0 
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where 0 is called null vector. 
Every vector space admits the existence of a subset 


B = {b1,... bn} CV 


called the basis of V. Thus, each vector © € V can be univocally 
represented through the basis B as follows 


v= vb; t= Lg? (1.3) 


where vt € JR are the components of è related to the basis B and n 
is a number which defines the dimension of V, namely the number 
of vectors in any basis of V. 

Notice that in equation (1.3) the Einstein’s summation conven- 
tion has been used. It is a notational convenience where any term 
in which an index appears twice will stand for the sum of all such 
terms as the index assumes all of a preassigned range of values, 
hence 


v= vb; = X oh (1.4) 
i=1 


1.1.2 Linear mappings 


Functions between two vector spaces assume a crucial impor- 
tance in linear algebra. In particular, we define a linear map as a 
linear transformation between two vector spaces that preserves the 
operations of vector addition and scalar multiplication. 

Let V and V’ be two vector spaces equipped with the bases 


B = {b1,... bn}, B= {b01,..., bm} 


respectively. 
We define a linear mapping as the transformation 


f:V>V', ted (1.5) 


if the two following conditions are satisfied 


1. f+) = f(T) + f) Vu,v € V : additivity; 


2. f (Au) = Af (u) VueVerA€ R: homogeneity. 


LECTURES ON SOLID MECHANICS 5 


The set of all linear maps from V to V’, denoted by L(V,V’), 
represents a n x m-dimensional vector space, where n and m are 
the dimensions of V and V’, respectively. 


{f:V-V}=:L(V,V) (1.6) 
For linear mappings the following properties hold 
1. (f+9)(@=f(@+g(d), ViigeL(v,VW);ueV 
2. (A) = Afa), vfeL(V.,V');aeV 


Matrix representation 


Notions so far introduced allow us to assert that if f is a linear 
mapping from V to V’, then f(v) is a vector in V’. Consequently, 
by recalling the expression in components for v, (1.3), we have 


f@=f0°d; i=1,...,m (1.7) 


and accounting for the fact that v = vb; with j = 1,...,n, and 
by using the homogeneity property for linear mappings, the latter 
equation leads to 


f (vib) Vi= if 0) Vi j=1,...,n i=1,...,m. (1.8) 
In a shorter form the components of f (v) are then 
(F 0) = fie? (1.9) 
so that the m x n-dimensional matrix fi = f (b;)' is the matrix 
representation of the linear mapping f referred to the bases B e B’. 
Linear forms and the dual space 


Linear forms are a special case of linear mappings. Let V be 
a vector space and B = {b;} its basis. A linear form w is a linear 
transformation from V to a scalar field 


w:V >R (1.10) 
Hence, we define V* as the set of linear forms from V to IR 


Vi =: {w: V > R} =: L(V, R) (1.11) 
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V* and V have the same dimension. 
The dual space V* admits a basis B* = {8°} whose elements 
are linear forms operating as follows 


8 @)=8={ 5 17) (1.12) 


By the definition, we can state that the element Bi belonging to 
B*, applied to the vector &, yields a scalar that is the i-th compo- 
nent of u. In fact we write 


6 (a) = B (ubj) = ut 6 (bj) = uls; — ul (1.13) 


We highlight that, as done for vectors, each linear form, chosen 
the n-dimensional basis B*, can be written in components as follows 


ww Plot (1.14) 


Bilinear forms 
We can define a bilinear form f as a mapping 
f:VxV> R, (0,0)5A (1.15) 


where 0,0 € V and A € IR, and such that it is linear in each 
argument separately. That is 


1. f (+0,0) = f(0,0) 


Endomorphisms 


Frequently in the field of solid mechanics we will meet special 
linear mappings from a vector space into itself, i.e. f € L (V, V). 
These are defined endomorphisms 


f:V>V, ted’ oveV (1.16) 


The set of linear mappings from V into itself forms a n x n- 
dimensional vector space, where n is the dimension of V. 


{f:V >V} =: L(V, V) =: End (V) (1.17) 
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Change of basis for endomorphisms 


Let B be a fixed basis for V, we are interested in evaluating how 
the endomorphism f € End (V) changes when passing to a new 
basis B’ of V. The following transformation rules are established 


bi = alt, (1.18) 
b, = alb; (1.19) 

that, by replacing (1.19) into (1.18), yield 
bi = al at by (1.20) 


and so 
(atak = öt) by = 0 => al'af = dl (1.21) 


therefore, each change of basis for V is characterized by a square 
invertible matrix n x n. 
Likewise vectors, the following rules hold for dual elements 


6 = à p” (1.22) 
B" = aig” (1.23) 


When both bases are orthogonal, then the transformation ma- 
trices are also orthogonal, that is 


a! = ai (1.24) 

where al? = (ab), and 
al = COS (bi, b;) (1.25) 
al = cos (bn, dk.) (1.26) 


The change of basis implies a change of the vector components. 
In fact we have 
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The proof of the above equations can be easily provided. For in- 
stance, for equation (1.27) we have that a vector © can be expressed 
with respect to two basis B and B' as v= v'b; = v” bi. Hence 


vib; = v ak by => viatby — vb; = 0 > (1.29) 
viatby — vidi bp =0> (via — vist) bp = 0 (1.30) 
finally, by putting zero the coefficient in brackets, we obtain relation 


(1.27). 


Covector components change by the the following rules 


vp = av} (1.31) 


vi, = ain; (1.32) 


Furthermore, recalling equation (1.9), via some manipulations, 
we get the rule to transform the endomorphism f, that is! 


fi = aj, fia’ Ik (1.33) 
and l l 
ff = an fea; (1.34) 
Similar relationships can be found for higher order matrices, for 
instance for a mixed fourth-order tensor we have 


ij Ilm nto 

hk = alal, no ap Q% (1.35) 
and likewise 

hnj _ til lm n 

nI = atdi fianna? (1.36) 


1.2 Euclidean spaces 


A Euclidean vector space is a space which admits a Euclidean 
metric, that is a structure inducing some special relationships be- 
tween distances and angles. In particular, fixed a Cartesian coordi- 
nate system (that will be better defined later on) and its standard 
basis, in a Euclidean space the distance between two points can be 
computed by means of Pitagora’s formula. 


Often, within an engineering context, it is convenient to represent equations 
(1.33) and (1.34) in the matrix form, such as F” = RT FR and F = RF’R’, 
where R” and R are nothing but a; and al, respectively. 
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1.2.1 Euclidean metric tensor and scalar product 


Let V be a n-dimensional vector space and B = {b;} be its 
basis. We define Euclidean metric the symmetric positive definite 
bilinear mapping 

g:VxV-R (1.37) 
that, given a pair of vectors u,0 € V, gives a real number g (u, 0) 


as follows 
u: v =: g (u,v) (1.38) 


The number g (ŭu, ©) is termed scalar product. The Euclidean 
metric allows us to compute distances. Indeed, we define length (or 
modulus, or norm) of & € V the real number 


Ilol] = 4/g (0,0) > 0 (1.39) 


The angle 9 amid vectors u and ® is given by the following 
equation _ 
giju v 

VIgijuivi||gijutvi| 

To compute the components of the metric tensor, i.e. the matrix 


cos v = 


(1.40) 


representing the mapping g, given the basis B = {bi} of V, the 
following general rule is adopted 


Jij = g (bi, bj) = bi - bj (1.41) 
that in the expanded form becomes 
bib o bieba 
Jij = i (x (1.42) 
bn: by +++ da da 


In the light of the above general expression for the metric tensor, 
the scalar product between two vectors becomes 


u-v= utb- wb; = uu! bj - bj = gue (1.43) 

Expression (1.43) includes, of course, the special case when, 
fixed a Cartesian coordinate system, the metric matrix equals the 
identity matrix dj; and consequently the scalar product can be car- 
ried out multiplying component by component, i.e. è -v = ulv! + 
«+ Nu”. 
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Now we want to point out that between the n-dimensional vec- 
tor space V and its dual V* there exists an isomorphism. Note that 
we are using some special words, e.g. isomorphism, without giving 
the formal mathematical definition. This lies beyond the purpose 
of this book, so that, also in this case, we will restrict the current 
exposition to an intuitive concept. From this point of view, an iso- 
morphism is a one-to-one mapping of an algebraic structure, e.g. 
vector space, into another of the same type, preserving all algebraic 
relations. 

Thus we define the musical isomorphisms: flat and sharp, re- 
spectively, as follows 


g VoV*:tHY 1.44) 
gi: Vi Viu o 1.45) 

where 
u(v)=g(u,v), VaeV (1.46) 


The isomorphism between V and V* implies the existence of a 
metric tensor 


gi:V*xV*— IR (1.47) 


so that 
u: v= 9(U,v) := g (u,v) =u-v (1.48) 


For further details the reader is referred to [4]. 

Both g and g are particularly helpful when carrying out com- 
putations it is necessary to switch from the contravariant form to 
the covariant form (and vice versa); namely when we need to lower 
or raise the indices. 


1.2.2 Eigenvalues and eigenvectors 


Let V be a n dimensional vector space and B = {b1,...,bn} 
the vector basis. Given f € End (V), we define the eigenvector a 
nonzero vector v whose direction does not change under the effect 
of f. Formally 

f=, AER (1.49) 


When equation (1.49) holds we can also define the real number 
A as the eigenvalue for v. 
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The eigenvalues of f represent the real roots of the following 
polynomial 


Pn (A) = det (fi — Adi) (1.50) 


where pn (A) is called the characteristic polynomial of degree n and 
fi is the matrix representation of the endomorphism f. 


1.3 Tensors 


This section is devote to a short outline of tensor analysis. 

Given two vector spaces U and V it is possible to construct a 
new structure, i.e. a third vector space, called tensor product of U 
times V that is symbolically denoted by U ® V. This vector space 
is made up of elements called tensors. It is possible to demonstrate 
that if 


Bg = {ŭ1,..., Un} 
By SUIT: 
are bases for U and V, respectively, then 
Boar = {0 8 Tj}, i=1,...,n;j=1,...,m 


is a basis of the vector space U@V. Therefore, each tensor 7 € U@V 
can be univocally expressed by 


T=T" (G @ d;) (1.51) 
where again the Einstein’s summation convention has been used, in 
fact (1.51) can also be written 


n 


m 
T= ) rigo; 
i=1 j=l 


1.3.1 Tensors and linear mappings 


The definition of tensors does not alter the structure of U and 
V, and, since the dual space V* preserves the structure of a vector 
space, we can introduce tensors belonging to spaces such as U* ® V* 
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and U*®V. In other words we distinguish the following second order 
tensors: 


U*®V: mixed tensors 
U* Q V*: covarinat tensors 
U Q V: contravariant tensors 


Mixed tensors: given the n-dimensional vector spaces V and V*, 
let a € V* be a dual form and d € V a vector, then the 
tensor a®v € V* QV can be identified by the endomorphism 
a QUE End (V) =L (V, V) defined as 


a@t:VoV:tH (a@d\ti=a(uyvevVv (1.52) 
Hence, a natural isomorphism has been obtained 
V*@VEL(V,V) (1.53) 


Covariant tensors: Let a, 3 be two linear forms belonging to V*. 
We can identify the tensor a ® $ € V* @ V* by the bilinear 
form a ® 0 € L? (V, R) defined as 


a@B:VxVOR: (Gt) alaB(v)eR (1.54) 
Therefore we can realize another isomorphism, which is 
V* @V* = 1? (V,R) (1.55) 


Vectors, linear forms and tensors so far discussed can be sum- 
marized in the following scheme 


Vectors 
dev (1.56) 
Linear forms 
ae V*=L(V,R) (1.57) 
Il-order mixed tensors 
a®vEeV*@V = End (V) (1.58) 


Il-order covariant tensors 


a®bBeV*®V*=L?(V,R) (1.59) 
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1.4 Coordinate systems 


Within the three-dimensional affine Euclidean space F it is pos- 
sible to define a coordinate system through the following bijections 
X:E> R OSI Se (1.60) 
where X = (ue The injectivity of X assures the one-to- 
one correspondence between points belonging to E and their co- 
ordinates. Namely, given a point p € E there exists the triplet 
(ae) which identifies such a point. The mapping X is as- 
sumed to be differentiable as many times as required. 
The coordinate system X is made up of coordinate functions 


g:E+R 4=1;2,3 (1.61) 


Moreover, we define the coordinate curves as the following map- 
pings 


tp: ROE i=1,2,3 (1.62) 
such as 
cip (A) =X (2 (p) + A, 2° (p) , x” (p)) 
Lop (A) =XT (xt (p) , 22 (p) + à, £? (p)) 


3p (A) =XT! (x (p) ,x? (p), £? (p) +A) 
that in a shorter form become 
ci (xip(A)) = ci (p) +X peE, AER (1.63) 


Given a point p € E, there are three coordinate curves passing 
through it. 

It is possible to demonstrate that the derivatives of the coor- 
dinate curves, computed for a fixed A, are vectors forming a basis 
B= {0;} in p. 

Analogously, it can be proved that the derivatives of the coordi- 
nate functions {2°} computed in p form a covariant basis B* = {d'} 
in such a point. 

The above two bases satisfy the following relation 


di (d;) = di (1.64) 
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Figure 1.1: Contravariant and covariant bases related to a 2D curvi- 
linear coordinate system. 


See |4] for further details. 

Bases B = {1, Oo, 03} and B* = {dt,d*,d31 related to X allow 
the representation of vectors, linear forms and tensor fields. For 
example we write 


Vo:E- E (1.65) 
w=wd' Vw:E- E (1.66) 


where v and w are vector and covector fields, respectively. 


1.4.1 Linear mappings and the metric tensor 


In order to represent an endomorphism f by means of the coor- 
dinate system X we can write 


f= fid o8, Vf:E-L(EE) SERE (1.67) 
where i l 
f? =d (f (8;)) : E > R (1.68) 
and likewise, for the bilinear form we write 
f= fyt d, Vf:E-L°(ER)SE ok (1.69) 
where 
fij = f (0,,0;): E > R (1.70) 


It is straightforward now to realize that the metric tensor g is 
nothing but the following bilinear form 


g: E > L (E, R) = E* ok (1.71) 
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indeed 
g = 94 @ df (1.72) 
where 
Jij = g (9;,9;) (1.73) 
As a concluding remark of this section we point out the fact that 
once the coordinate system is fixed it is possible to find its vector 
basis, i.e. the covariant basis, and therefore the covariant expression 
of the metric tensor can be directly computed. 


1.4.2 Components of the metric tensor 


Suppose that Xe = {xt}, i = 1,2,3 is a Cartesian coordinate 
system, with the origin o € E, which describes the affine Euclidean 
space E and {é;}, i = 1,2,3 its unit normal basis. Moreover, let 
X = {x}, j = 1,2,3 be a generic curvilinear coordinate and 
{0;}, j = 1,2,3 its basis. Suppose that the functions xi and q’ 
are single-valued and continuously differentiable with respect to 
each of their variables as many times as required, we can therefore 


write 
xi = qi (a a) i=1,2,3 (1.74) 
g = a (xl, x2, x3) i= 1,2,3 (1.75) 
and the rules for changing basis (1.18) and (1.19) on page 7, become 
> ðx! . Oat h 
ði = Ori oh dre Inke (1.76) 
and 
Og 2 io Orn 
C: = —— X = d 1.77 
Ci Öri h; e Aa ( ) 


where equations (1.76) transform the covariant and contravariant 
elements of the Cartesian basis into the elements of the curvilinear 
basis while expressions (1.77) perform the vice—versa. 

Now, according to equation (1.41), it is possible to compute the 
covariant components of the metric tensor related to the curvilinear 
coordinate system 


a Oxh | Oxk _ 
sita» x? On 


~ da dai "AT Axi Oxi (9) 
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Moreover, the contravariant components are 


= pe = 1.80 
Ù ~ Oh” Bak È ey) 
Ox' Ox) pp Ox dai 
= = 1.81 
Ort dak Ort Oxh SEGA) 
and finally the mixed components of the metric tensor are 
. E dx Oak 
i — d (d,;) = — el. Ce, = 1.82 
_ as dre h_ de Ore (1.83) 
ox? dx) Ox! dx) 
Christoffel symbols 
The Christoffel?’s symbol are defined as follows 
Ty, =d"(Vid;):E- R (1.84) 


so that Vid; z T% bk. Hence, r$, is the k-th component of the 
derivative of the basis element d; along the i-th direction. Analyt- 
ically they can be computed by the following formula 


1 
ri = 59 (aighi + djIni — Ongiz) (1.85) 
where ð; denotes the partial derivatives. 
Moreover, it can be proved that 
rk, = (V,d;)" =- (Vid) (1.86) 


J 


For proofs and more details the reader is referred to |4], [5] and 
[1]. 


?Elwin Bruno Christoffel (November 10, 1829 Montjoie, now called Mon- 
schau - March 15, 1900 Strasbourg) was a German mathematician and physi- 
cist. 


Source: http://en.wikipedia.org/wiki/. 
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1.4.3 Examples of coordinate systems 
Cartesian coordinates 


The Cartesian coordinate system introduces considerable sim- 
plifications with respect to other curvilinear systems, e.g. cylindri- 
cal, spherical, hyperbolic, etc. 

Therefore, let us begin by defining a Cartesian coordinate system 
as the triplet of coordinate functions 


Xe = (x,y,z) = (xt, x°, x°) : E > R? (1.87) 


with an origin in o € E and equipped with the unit normal basis, 
also called standard basis, {€1, €2, 3}. Given p € E, the coordinate 
functions are such that 


x! (p) =: (p—0)-&; (1.88) 
The coordinate curves of a Cartesian system are 
Lip (A) = p + Ae; (1.89) 


Notice that for rectangular coordinate systems the symbols de- 
noting the bases will turn into 


DI 


Ci (1.90) 
e (1.91) 


= 
Il 


The covariant form of the metric tensor can be readily computed 

as follows 
Gig = & Ej = di; (1.92) 
Elements of the standard basis related to the Cartesian coordi- 


nate system do not vary with the point p € E. As a consequence, 
the Christoffel symbols are identically null. 


TË =0 (1.93) 


In addition to that we also highlight that here the upper or lower 
position of the indices does not influence the structure of the field 
we are dealing with. Namely, vectors and linear forms are the same 
and the unit normal basis equals its dual. 


g (Ei) =e = è (1.94) 
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For this reason whenever given two sets of numbers having the 
same dimension, they can be ordered in a row and a column, respec- 
tively, and by using the multiplication rule row—by—column a scalar 
is always yielded without taking any care whether we are dealing 


with vectors or linear forms. 


Cylindrical coordinates 
We define a cylindrical coordinate system the functions 
X=(0,9,2):E- R? (1.95) 


In this case, with the help of figure 1.2, equation (1.74) becomes 


x = psenv 

y = pcos 

Paz 

AZ 
io, 
zA 
©, 
é, > > 
a y 
tr 


Figure 1.2: Cylindrical coordinate system. 


Now, through equation (1.76), it is easy to compute the basis 
related to the cylindrical system 


5 Or. Poi gl 
gp ap? ə €3 
= ôx Oy _ Oz _ 
Og = It 397 99° 
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Hence, the covariant components of the metric tensor is 


Jepp = 1 

_ 2 
G99 = P 
gzz = 1 


Joz = Ipz = Ips = 0 


that in the matrix form can be written as follows 


1 0 0 
gj=| 0 0 0 (1.96) 
0 0 1 


The contravariant form of the metric tensor is 


E 
Pau Mr (1.97) 
ioe 


Using equation (1.85), the Christoffel symbols are 


1 
Me =p. Fi, =1 = 7 (1.98) 


Spherical coordinates 
We define a spherical coordinate system by the functions 
X =(r,0,y): E > R? (1.99) 
In this case, with the help of figure 1.3, equation (1.74) becomes 
x = rsin y cos Ŷ 
y = rsin ysin V 
Z = r COS 9 


Now, through equation (1.76), it is easy to compute the basis 
related to the spherical system 


= ôr_ ðy oz _ 
Or Ape FT ar? att dr 3 
=- OL x dy_ dz. 
IT 991 Bg? ag? 
5 = dx _ pa a 
yp = ap! ap? E) €3 
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cv 


Figure 1.3: Spherical coordinate system. 


Hence, the covariant components are 


Grr = 1 
999 = r? sin? p 
Ippo = n 


Grd = Gre = Ivy = 0 


which in the matrix form are written as follows 


1 0 0 
Gg | 0 sin yo 0 (1.100) 
0 0 r? 
The contravariant form of the metric tensor is 
1 0 0 
Psu e A (1.101) 
0 4 
ze 


Using equation (1.85), the Christoffel symbols are 


Fog == f= =r sin? p 
1 
= z= di — eg 
Por o n Tag = — sin Y cos y 
1 
Go po _ = 9 p9 _ cosy 
Tra = I'r = r Toy =r ~ sing 
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1.4.4 Volumes and the vector product 


In the three-dimensional Euclidean space a volume element n is 
defined as a three-linear form such as 


= ExExE- hR, (1.102) 
(U1, U2, U3) E n (01, 02,03) € R (1.103) 


3 


When the set of three vectors forms a basis {b1, b2, 63}, for any 
other basis {b} , 65, bs}, the volume element changes with the follow- 
ing expression 


n (b1, b2, 63) = |a}m (b1, ba, b3) (1.104) 


where Ja'i] is the determinant of the endomorphism for basis chang- 
ing already seen in equation (1.18). 

The application 7 can be expressed by a third order skew- 
symmetric tensor Nijk with the following properties. If two of the 
subscripts {i,j,k} equal each other the volume element vanishes. 
Any odd permutation of the subscripts changes the sign of the el- 
ement, any even permutation of the subscripts does not alter the 
volume element. 

For a Cartesian system of coordinates we shall denote the vol- 
ume form by ¢;;, and the above properties become clearer in the 
following scheme 


O when any two of the indices are equal; 
Eijk = gue = 1 when i,j,k is an even permutation of the numbers 1,2,3; 


—1 when i,j,k is an odd permutation of the numbers 1,2,3; 


that means, for example 


112 133 222 
€ €112 = € €133 = € €222 = 0 
123 231 312 
€ €123 = € €231 = € €312 = 1 
132 321 213 
€ €132 = € €321 = € €213 = —l 


In addition, the operator €;;, satisfies the following identity 
€igk€ith = ÔjlÔkh — ÕjhÔkl (1.105) 


Suppose that {0;,0;,0,} is a basis related to a curvilinear co- 
ordinate system, we want to evaluate the volume element in this 
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system through the volume element expressed in the Cartesian axes. 


E On) ðr DIE 

n (dì, j, dk) = & er, dx; €s, dr, a) = 
Ox, Ox Oxe nt 
di On; Ox, Erst = Nijk 


(1.106) 


which, by means of equation (1.104), we find that the coefficient of 
the volume element in the Cartesian coordinates e on the left-hand 
side of the latter equation equals the determinant of the endomor- 
phism of the basis changing, namely 


Oak Ax’ dx. da 
Ox; Ox; dark = [eal det (5) SIOE) 


It is easy to prove that the above determinant is the square 
root of the determinant of the metric tensor related to the generic 
coordinate system fail , i = 1,2,3. So we have 


n (3i, 3j, Ox) = (/I9pgleisk (1.108) 
and in the same way the following contravariant expression can be 
derived o "i 

n (di di, dt) = Vigre (1.109) 


lots 
detto) — det (999). 

The skew-symmetric tensor 7 defines the vector product as 
follows 


where |gpg| = det (gpq) and |g?4| = 


ū x b= ud; x vd; = uwi ð; x dj = uv) ni; (1.110) 
and also 
uxv=u,d' x v;di = ujv,;d’ x di = unin” Ox (1.111) 


We can use the tensor 7 to compute infinitesimal volume, area 
and line elements. Let us begin putting the infinitesimal vector 
along the j-th coordinate curve as follows 


dl; = do! d; (j not summed) (1.112) 
so that the infinitesimal volume is given by 


dV = n (dl, dla, dls) (1.113) 
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hence 
dV = n(0,,d,,d) dalda?dx3 = 
Vge23dr! da dr? = \/gdx'dx*dx? (1.114) 


The above expression for the volume element can also be written 
as 


dv = [dl, x dla] - dl3 (1.115) 
that allows us to attain the same equation expressed in (1.114), 
indeed we have 


[dl x dlə] 5 dls = dx!dr?nz3d? (dx*0s) = 
da! dx? dx? \/ge123d° (ðs) = dal da?dx3 /ge12393 = 
Vgda! da dr* (1.116) 


Infinitesimal area element 


Taken two infinitesimal vectors along two coordinate curves re- 
spectively, the infinitesimal area normal to the vector along the 
third coordinate curve is given by 


dA; = |dl x dl2| = m23]d*|daldx? = 


vg] d3 - d3da!da? = Vgg83dx! dx? (1.117) 
and it is easy to obtain the general expression for any area element 
A; = Vggiidai dz" (1.118) 


where è is not summed and i # j # k. 


Infinitesimal line element 
A generic infinitesimal line element dl? is defined by 
P = |d]? = dl - dl = 
dx" 0; - del d; = dx'da? gi; (1.119) 
whereas, a line element taken along the è — th coordinate curve can 
be represented by the vector 


dl; = dx'd; (i not summed) (1.120) 


and it measures 


g (dl;, dl;) = guds" (1.121) 
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1.5 Covariant differentiation 


In this section we shall briefly introduce some notions concerning 
the derivatives of objects so far discussed, i.e. vectors and tensors. In 
order to differentiate these fields the concept of manifold is required. 
However, in this context it will be restricted to a rough and informal 
description. 

A manifold is an abstract space locally Euclidean so that, for 
each point belonging to the manifold, there is a neighborhood that 
can be described as the Euclidean vector space. When we deal with 
manifolds, the intuitive idea of vectors obtained by simply subtract- 
ing two points in the affine space might no longer be valid. Keep 
in mind, for instance, a curved surface Q € E, i.e. a two dimen- 
sional manifold, and try do define a vector entirely belonging to the 
surface by subtracting two points. It is easy to see that the vector 
cannot belong to the curved surface Q. 

For this reason we need an additional space named tangent space 
TE that allows us to extend the concept of vector spaces so far 
discussed to manifolds. The tangent space is a Euclidean vector 
space consisting of the tangent vectors of the curves through the 
point of the manifold itself. 

In order to use tools for computing volume, area and line el- 
ements, i.e. to define the metric tensor, we shall suppose that we 
always deal with differentiable Reimannian manifolds. For a formal 
mathematical definition see [5]. 

Given a general coordinate system X = tert , i = 1,2,3, let 
a be a vector field a: E > TE and 7 : E > @*TE a k-order 
contravariant tensor, we define the covariant derivative V.,7 of the 
field 7 with respect to the field © as 


Var = ul (ðr + Bee +... + Ti piene, @.... @ Oi, 

(1.122) 

Analogously, for a k-order covariant tensor 7 : E > QT E” the 
covariant derivative becomes 


Var = Optic = nina oT ain)" dad 

(1.123) 

where T'E* is the cotangent space, namely the space that contains 
the dual forms related to the vectors belonging to TE. 

The above expressions are presented only for the sake of com- 

pleteness, while, the covariant derivative of vector fields and second 
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order tensors, will be often used in the mechanics of shell contin- 
uums. In fact, for a second order covariant tensor 7 = r'fd, Q Oy 
the derivative is 


Var = uf (ajr + Thr” + Thir") 04, (1.124) 
while, for covariant tensors T = hed” ® dÝ the derivative becomes 
Vat = uf (3jThk — Tinta — Ty the) d” © d* (1.125) 

and for a mixed tensor T = tha” ® dp the derivative is 
Var = uf (art + ") d'a ô, (1.126) 

Finally, for a vector field we have 
Vin (00° i Tino) a; (1.127) 
and for the dual form 


Vasy (yx = Tion) di (1.128) 


1.5.1 Grad, div, curl and Laplace’s operator 


Gradient. Consider a scalar field f : E — JR, we define the gradi- 
ent of f as the the vector 


grad f = gi Fa, (1.129) 


In a Cartesian coordinate system the above operator simplifies 
in the following expression 


grad f = ore, (1.130) 


Divergence. We define the divergence of a vector field v as the 
following scalar 

div 0 = tr (V0) =v}, + Tijv (1.131) 

In a Cartesian coordinate system the divergence is written as 


divo = tr (Vd) =v, (1.132) 
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Curl. For the sake of simplicity, to define this operator let us de- 


note by V a symbolic operator defined as V = i d°. Now the 


curl of a vector field v can be defined as the following vector 


culo = Y x g’ (7) = 


ə. sed DA | se 
= dat X vj d? =d'x Dai (v;d) = 1 vy OK (1.133) 
where n is the skew-symmetric tensor related to the vector 
product (1.111) and v;; stands for the covariant derivative 


ae h 


|j 


Hence, for a rectangular coordinate system, the curl assumes 
the straightforward expression 

curl = V x v = vja E = (1.134) 
where @ = wyex and wk = vie". Expanding the latter ex- 
pression leads to the following equivalent form 


er € €3 
curl v = det 3, st 3, (1.135) 


VI U2 U3 


Laplace’s operator. We define the Laplace operator of a scalar 
field f the following scalar 


V2 = gii (0:0;f- 40,5) (1.136) 


In a rectangular Cartesian coordinate system the Laplacian 
is written as 
Pf Ef Ef 


yep 
f Ox? j ðr? ôx? 


(1.137) 


A useful remark 


For practical uses of the above expressions concerning the differ- 
ential operators it is necessary to consider a unit system. Attention 
to this aspect must be especially payed when curvilinear coordi- 
nate systems are involved in our computations. For clarity’s sake 
we recall both expressions for a vector 


b= vd; (1.138) 
v = vd (1.139) 
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and we point out that vectors (forms) forming the covariant (con- 
travariant) basis {d;} ({d’}) are not dimensionless. Hence, the vec- 
tor components do not represent a physical quantity, even though 
their geometric properties are correct. So, in order to give vec- 
tor components a physical meaning a normalization of the basis 
is required. To this end we introduce the so called physical basis 
{O<i>} such as 


v= vd; (1.140) 
Next we normalize the covariant basis as follows 
_ È; ð; 
a (i not summed) (1.141) 
|di| Jii 


which, replaced into equation (1.138), allows us to define the phys- 
ical components of 0 as follows 


us = [guv (i not summed) (1.142) 


On the other hand for the dual basis we have 


d<” = (i not summed) 


[gii 
Vei> = viy g” (i not summed) 
As an example, in the following we present the expressions of 


the differential operators discussed in section 1.5.1 for a cylindrical 
coordinate system. 


e Gradient 
Of = oo OF Of = 
— pp CI I zz = 
grad f =g Dp oe + 9 3590 -+ Dad 
_ Of = 1 ðf = Of = 
Fup oP F pd 5,9 08 * DE Og 
Of = 10 
= dI cys -3 Lg + 5 50<z> 


e Divergence 
1 
div y = tr(Vv) = vf, + vh +u%, + o = 
x > 7 p 


= l Cag + vg) + U 
P ? , 
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e Curl 


curl v = n Fu Ok +97 vak + n vsk 
— — 
=0 
n'v pd + n vaak +1 vs28k 
— — 
=0 
1 vð + n? vað + n° v3|3k = 
— 


=0 
TN, |g*| ( (v32 2 vaz) ðı + (vijs = U31) 03 + (va = V1)2) 93) 


which by making use of the cylindrical notation as stated in 
section 1.4.3 (i.e. 1 = p, 2 = J, 3 = z), taking into account 
that vy; = vij due to the symmetry of Christoffel symbols 
in equation (1.98) and considering the physical components, 
allows the above expression to become 


1 Di 
curl v = (Fre = vev») <p> + 
Z 1 £ 
+ (V<p>,z — V<z>,p) Ic9> + (U<4>.0 — V<p>,0) O<z> 


e Laplacian 


Canin 1 of OT % 10f 
_ pæ pV Əz pp Ap 


Note that in the latter expression no normalization has been 
used. 


The divergence theorem 


Consider a generic region V C E bounded by the smooth closed 
surface S. Given a continuously differentiable vector field v € V, we 


have 
faivvav= | onas (1.143) 
V S 


where 7 is the outward pointing unit normal vector of the boundary 


S. 


LECTURES ON SOLID MECHANICS 29 


In components the above theorem becomes 


[+ dV = [vin ds (1.144) 


The divergence theorem holds for tensor fields. For a mixed II- 
order tensor T = a (a* Q On), for example, the theorem states 


[evray = [rm dS (1.145) 


where the k-th covariant component is 
I (mtn + That — Dia) dV = i, Tiny aS (1.146) 
v S 


While for a Il-order contravariant tensor T = 7!" (On ® dk) it 
becomes 


f (i + Dir + rier") dv = f ren, dS (1.147) 
V S 


1.6 Affine space 


Here we shortly introduce the notion of affine space. 

Let E be a n-dimensional vector space. We define the affine 
space associated to E the set of points E equipped with the trans- 
lation +, such as 


+:EXESE:(p,t)+pt+u=p EE (1.148) 
where © = (p'— p) € E represents a free vector, while the pairs 
(p, ©) form applied vectors. 

1.6.1 Free and applied vectors 


This section is restricted to the geometrical interpretation of 
vectors belonging to the Euclidean space and expressed through 
the rectangular coordinate system. So that we have 


gij = JÏ = gi = di; (1.149) 


and 
Nijk = Eijk (1.150) 
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From a geometric point of view an applied vector is represented 
by a line segment AB from point A to point B, where, with respect 
to equation (1.148), A = p and B = p’. If B is moved to the position 
C, then the whole translation from A to C represents the sum of 
the partial translations AB and BC. 

Putting AB = a and BC = b we notice that if they were applied 
in the same point, see figure 1.4, then a practical rule can be used 
to carry out the addition a + b. It consists in moving the vector b, 
in such a way to be kept parallel to itself, into a new position so 
that its starting point coincides with the ending point of a. Thus, 
the line segment from A to the end point of b (in the new position) 
represents the addition a + b. See figure 1.4. This rule is known as 
parallelogram rule because a and b form the sides of a parallelogram 
and a + b is one of the diagonals. 


Figure 1.4: Addition of two applied vectors. 


The subtraction of two vectors applied in the same point can 


be seen as € = @+ (-5) and so it is carried out by means of the 
procedure described for the addition. The vector c = a — b will be 
given by the line joining the starting point of a to the end point of 
—b. See figure 1.5. 

The addition of two applied vectors has the following properties 


1. a@+6=6+4; 

2. (@+b)+e=4 
3. (A+ p)a= àa 
4. A (ua) = (Ap) a; 

5. A (a+b) = àā + àb; 
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Figure 1.5: Subtraction of two applied vectors. 


where a, b, € are applied vectors in E and A, u € R. 
For an applied vector it is possible to define norm, direction, 
sign: 


norm (modulus or length): is the length, measured by a fixed 
— 


unit system, of the line segment AB; 
direction: is the direction of the line passing through A and B; 
. . . . È = 
sign: specifies the sign, i.e. AB = -BA. 


From the preceding discussion about the metric tensor it is 
—> 
known that the length (modulus) of a vector a (= AB) is the square 
root of the scalar product by itself 


| a |= /9 (4,4) = Va-a (1.151) 


Recalling that the metric tensor is a bilinear symmetric positive 
definite form, the following properties can be derived 


1. @-a@=||a|?>0sea40; 
2. a-b 

3. c-(a+b)=c-a+c-b; 

4. A (a-b) = (Aa) -b =4:- (àb); 

The cross product of two applied vectors a, b € E in a Cartesian 


coordinate system is carried out by using the general rule given in 
equation (1.110), so that 


v=axb weEeV (1.152) 


where 
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modulus: || © ||=|| a ||| è || sin@, where @ denotes the angle 
between a and b; 

direction: normal to the plane to which ā and b belong; 

sign: follows the right hand rule. 


Moreover, by virtue of the skew-symmetric tensor €;j;x, the vec- 
tor product vanishes when either one of the two vectors vanishes or 
when the two vectors are parallel. See figure 1.6. 


A 


w 


SI 


-wW 


M 


Figure 1.6: Vector product for Cartesian applied vectors. 


The following properties can be also enunciated 


2. (Aa+ ub) x E= A (ā x 2) + u (b x 2). 


1.7 Surfaces 


Let E be the affine Euclidean space. The submanifold Q C E is 
a surface if dim Q = 2. 

Suppose Q C E is a surface which can be described by an in- 
duced coordinate system of dimension q = m — k, where m is the 
dimension of E and k denotes the number of constraints (codimen- 
sion of Q). Since Q is a surface we have m = 3, k = 1, q = 2. The 
induced coordinate system is given by 


Xİ: Q> Ri: pH «%(p) (1.153) 
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From now on the quantities living on Q will be distinguished by 
the symbol { and the components will be written using superscripts 
and subscripts, running from 1 to 2, in Greek letters. The Latin 
indices will denote components of quantities that are applied on Q 
but lie out, namely belonging to the vector space TgE. 

The unit normal vector is defined as follows 


n:Q—> TQ" so that g(ñ,n)=1. (1.154) 


where g is the metric tensor defined on TE and TO is the orthog- 
onal space. 
Analogously, on the surface Q we can define the the induced 
metric as 
gi:TQxTQ- R 
that in components? becomes 


g' = 9agd® ® dl 


Given two vectorial fields u : Q — TQ and 0: Q — TQ, the 
covariant derivative of v with respect to & can be split as follows 


Vað = Vlo + via (1.155) 

where 
V! :TQxTQ—=TQ (1.156) 
Vt : TQ x TQ > TQ (1.157) 


The application V!! is called second fundamental form of the 
surface. For further details see [5] and [3]. 

We now define the Weingarten* map L as the following endo- 
morphism 


L:=Vn:TQ>TQ:ŭ= Van (1.158) 


3In some books the covariant components of the metric tensor g! are also 
denoted as gir = E, gi2=F, g22 =G. 

‘Julius Weingarten (March 2, 1836 Berlin - June 16, 1910 Freiburg) was a 
German mathematician. 


Source: http://www-history.mcs.st-andrews.ac.uk/. 
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In addition to that, we define the total curvature (Gauss curva- 
ture) K and the mean curvature H of a surface Q as follows 


K :=detL:Q- R (1.159) 
H :=trL: Q > R (1.160) 


Finally, eigenvalues of L are defined principal curvatures. See 


[5]. 
Let L be the second order covariant tensor related to the Wein- 
garten endomorphism L by the metric tensor g', so that 


L:=Vn:TQxXTQ— R: (u,v) = g (L (u), 0) = Van-d (1.161) 


where n = g? (ñ). 
The following differentiation 


0 = Va (g (5, 8)) = g (Vað, A) + g (0, Van) > (1.162) 
g (Vad, n) = —g (U, Vai) (1.163) 


proves that the scalar quantity L (u,v) represents the normal com- 
ponent to the surface Q of the covariant derivative, namely 


Vað = Vlo — L (a, 0) ñ (1.164) 


Dealing with mechanics of shell continuums, equation (1.164) 
will be often used. Hence, in the following we expand its expression 
in components. 

Suppose {ða}, a = 1,2 is a basis related to the induced coordi- 
nate system describing the surface, we have 


Vga = V} 3a — L (ds, Ia) N (1.165) 
and for both right hand terms we have, respectively 
Vh ða = 4" (95) (A (4° (In)) + P20 (In) A (1.166) 


= 6} (72,02) 6, = T2, ð, (1.167) 


L (8g, a) = (L (95) da) = Va; da (1.168) 
= LO. da = Le Jwa = Lpa (1.169) 
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Finally, equation (1.165) in components becomes 
V 0a = Tgad — Lpañ (1.170) 


Note that in the remainder of this book, for the sake of brevity, 
we will use Vg: instead of V5," 

Analogously, for an element of the controvariant basis, recalling 
the general equation for covariant derivatives, and considering the 
above Gauss splitting, we have the following expression 


Ved* = -19,d* — Len (1.171) 


Often, for instance in the case of shell theory, we will deal with 
vector fields that do not belong to the tangent space, so it is useful 
to present an example of derivative of vectors applied in Q but lying 
out of the tangent space. Namely, suppose that v € TgE. We can 
decompose the field v into the tangent and orthogonal component 
as follows 

a = öl + gt (1.172) 


that in components is written as 
T= vd, + eH (1.173) 
Hence, given u € TQ the derivative of è with respect to @ is 
Vað = Val + Vat = Viol — L (a, all) n+Vadt (1.174) 
that in components turns into 
Vad = uf (090% + 19,07 + 0653) Ja +u (vf — Lage”) a 
(1.175) 


In the same way, the dual form v € TOE can be differentiated 
as follows 


Vav = Vaul + Vaut = vio — L(a, vll)n + Vavt (1.176) 
that in components becomes 


Vav = uf (agva — Taty + veL) d7 + uf (veg — LGva) n 
(1.177) 
Examples of surfaces will be provided in appendix A, where, 
within the application of the shell theory, the above results will be 
applied to some well known geometries. 


Chapter 2 
Analysis of strain 


This chapter is devoted to the classical strain theory for deformable 
continuums. In order to offer a comprehensive approach, the first part will 
be treated in curvilinear coordinates, then results in Cartesian coordinates 


will be obtained as a special case. 


2.1 Introduction 


Before introducing the definition of strain it is useful to give 
some preliminary concepts and definitions. 

Let us begin with the definition of body. 

A body C C E consists of a set of particles embedded in the 
three-dimensional Euclidean space. Each particle p € C, i.e. a ma- 
terial point, can be put in one-to-one correspondence with a triplet 
of scalars that univocally determine the position of such a point. 
Namely, for any point p included in the body there exists a coordi- 
nate system X : C C E — IR?. See also the more general expression 
(1.60) on section 1.4. 

From the notions of body and time we can derive the concept 
configuration. Configurations are regions V of the three-dimensional 
Euclidean space E that can be occupied by the body in a particular 
instant. Thus we have 


V = (C,t) = {(p,t) |p e C} (2.1) 


where VY is also called a spacial domain for fixed t. 
It is assumed that: 


e Configurations are open connected sets or domains in the Eu- 
clidean space. 


e On varying of the time t, the configurations of one and the 
same body maintain a continuous one-to-one correspondence 
between different positions of one and the same particle. 


Claudio Borri, Michele Betti, Enzo Marino (a cura di), Lectures on Solid Mechanics, ISBN 978-88-8453- 
854-3 (online), ISBN 978-88-8453-853-6 (print), © 2008 Firenze University Press 
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2.2 Deformation 


Now, beginning with an intuitive statement, we can introduce 
the definition of strain. When the relative position of two points 
included in a continuous media is altered, we say that the body is 
strained. Hence, analysis of strains means to evaluate the change 
of the relative distance between points; this is also called deforma- 
tiont. 


2.3 Strain tensor in general coordinates 


Let V be the region taken by an unstrained state of a body at 
time t, so that 
V= (C,t) (2.2) 
and V’ the configuration of the body in the strained state at instant 
t', that is 
V=(C,t’) (2.3) 
Consider a Cartesian coordinate system equipped with the unit 


normal basis {ē;}, so that for any point p in V and p’ in V’ the 
positional vectors can be written respectively as 


7 = (p—o) = gie; 
F = (p-0)=ye 


We assume that each point in V' is related to its original position 
in V, and vice versa, by the following relations 


SO, 
idea) 
dave) 
In order to avoid penetrations or separations of the material 
particles it is necessary that the transformation of points in VY into 


points in V’ is one-to-one. Namely, to ensure the existence of the 
single-valued inverse of equation (2.6) (or (2.7)) it is sufficient to 


‘We know that in nature all materials are deformable, but sometimes we 
will refer to the abstraction of non-deformable (or rigid) body. This abstraction 
assumes that for every pair of points belonging to the continuum, the relative 
distance remains unvaried throughout the history of the motion. 
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assume that the functions y% and xt are continuous and differen- 
tiable as many times as required and the Jacobian is greater than 
zero”. We write, accordingly 


Consider now a generic curvilinear coordinates system X = {x°} 
so that 


F = 2d; = xd! (2.8) 
where {d;} and {d°} are the covariant and contravariant bases re- 
lated to the curvilinear system and x° = r’, ri = xi. See figure 
2.1. 


Points belonging to the initial configuration V can also be re- 
lated to the curvilinear system of coordinates as follows 


ri = ot (xt, x°, 2°) (2.9) 


where zt are single-valued and differentiable as many times as re- 
quired. 

Moreover, we can use the curvilinear coordinates to describe the 
body in the strained configuration V’, so that 


Ye = ye (x', 27, 2°) (2.10) 


According to section 1.4.2, through the Jacobian matrices, we 
can compute the metric tensors g and g’ associated to the curvilin- 
ear coordinate system for both configurations, respectively. 

For the unstrained configuration the covariant components of 
the metric tensor are 


= & Dei. ‘000. 
Jij = ði È 0; = Oxi P s Əri * = (2.11) 
dx? Oak dx dat 


area (2.12) 


?The Jacobian of the function yi = yi (xå) is the determinant of the matrix 
whose i — th row lists all the first-order partial derivatives of yi. 
*With the exception of singular points, curves, surfaces. 
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Figure 2.1: Unstrained and strained body states. 


while the contravariant components are 


edad = pe = 2.13 
Perea pt (2.13) 

Ox' Ox) pp Ox dai 
= = 2.14 
dah Aak Ort Aah Cnt 

and finally the mixed components are 
; Bos Ox" Ox 

i_ ygi a h, È aai x 
gj =a (0;) = dat or” (2.15) 
_ On Oak p Ort dal (2.16) 


~ deh dai E Axh Axi 
For the strained configuration the covariant components of the 
metric tensor are 


5 Oy OyE 

Jij = dj dI Dai Eh Dai = (2.17) 

Or' Oxi da Ori 
the contravariant components are 

Di | . ð Oxi 

g” =q" g” = n : aye = (2.19) 
e 

Ort Ox! pp OL OLİ (2.20) 


pag, Ay’ Ayh 
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and finally the mixed components are 


dx! h Oyt Til 
Ayre dg ET 
OR OU p Ox Oy? 
~ Oyk Aad E Ayh Aad 


fede (2.21) 


(2.22) 


At the beginning of this chapter we said that the aim of the 
analysis of strain is to evaluate the change of length between two 
points in a continuous medium. We are now mathematically able 
to evaluate this difference 


dl? — dl? (2.23) 


where dl’? = |d} |? and dl? = |dl|? are the arc lengths of the strained 
and unstrained states, respectively. Namely, the vector dl, joining 
the points po and p, during the the transformation, is carried into 
another vector dl’. These two vectors differ in direction and magni- 
tude. See figure 2.2. 

By using equation (1.119) on page 23, we can write the above 
line elements with the help of the metric tensors for both configu- 


rations as 
dl? = gi;jda'da! (2.24) 
dl? = gi dada! (2.25) 
then, the difference 
dl? — dl? = (gi; — gij) da' da? (2.26) 


We now define a symmetric tensor named the strain tensor, as 


1 


so that e 
dl” — di? = 2y;;da' da! (2.28) 


The strain tensor is obtained by subtracting two bilinear forms 
so that it is still a bilinear form. Therefore, given two vectors po 
and go at a fixed time to (let us say the initial unstrained state), 
the strain tensor just measures the difference between the scalar 
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Figure 2.2: Measure of strain. 


product of the vectors p and q at a generic time t (that identifies 
the strained state) and the scalar product at the initial state. 


y:ExE-R (2.29) 
(P,Q) > (0,9) = 9' (P, 3) — 9 (P, 3) (2.30) 
so that 
y (P, a) = pat; (d' ® dI) (On, Oe) = 
= par; = VP’? (2.31) 
Points in V and V’ are univocally determined by the positional 


vectors 7 and 7” respectively. With respect to the generic curvilinear 
coordinate system X we have 


F= xð; (2.32) 
7 = y'ð; (2.33) 
hence, the position p’ relative to p is denoted wu and it is called the 


displacement vector 
u= -T (2.34) 


Considering now that the basis related to the curvilinear coor- 
dinates is given using equation (1.76), we have in the following an 
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equivalent expression 


ð; = or (2.35) 
=} 
dv = z (2.36) 


and by considering the relation (2.34) it becomes 
ð, =r; =F i+ Vit (2.37) 


where the comma denotes the partial derivative and V indicates 
the covariant derivative. 


la 24% 
Vig = z (8+ 3; — 3: - 5) = (2.38) 
1 = Fey 
= 5 ((a + Vit) - (0; + Vja) — ð; -0;) = (2.39) 
at d; -Vju+0;- Vint Viù- Vit (2.40) 
2 


In fact, recalling the general expression (1.122) for this differen- 
tiation, the above equation turns into 


=r, = + (ur + Dîu j Öm (2.41) 


where the Christoffel symbols are referred to the metric tensor re- 
lated to the curvilinear coordinates for the original configuration V 
of the body. 

Finally, using the definition of strain tensor, with some calcu- 
lations we can obtain the expression of the finite strain tensor in 
general coordinates as 


1 . , 
Jij = 9 (viv! + Viu + Viu” V ju") (2.42) 


Expanding the above derivatives the strain tensor assumes the 
following expression 


1 
3 (uis + Uji + Trun) + 


1 
5 (unit + up Thu" + up iF, u S DR upDh e) (2.43) 


Vig = 
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For Cartesian coordinate systems we could repeat exactly the 
above procedure to obtain the strain tensor, but this is equivalent 
to putting zero all the Christoffel symbols in equation (2.43). So 
that for rectangular coordinate systems the strain tensor assumes 
the following expression 


1 
5 (Vij + Uji + UkiUk j) (2.44) 


I= 5 


where we remind the reader again that in the rectangular coor- 
dinates the position of the indices does not make any difference 
because gij = dij. 
2.3.1 Examples of strain in Cartesian coordinates 
Stretching ratio 
Let us define the stretching ratio 6; as follows 
d'—- d dl 
dl dl 


Namely, suppose we have two points in the unstrained state the 


ô = 


1 (2.45) 


difference of which gives a vector dl = dztē;. The corresponding 
vector in the strained state is dl’ = dx"g!.. Therefore the stretching 
ratio 6, gives the relative difference between the length of the vectors 
dl and dl’. 

By means of the definition (2.28) we have 


VWjda'dxi — dl” 
dxkdxk dł 


(2.46) 

then 

Vijda'dxi 
dakdx* 


so that the stretching ratio can be written as follows 


Lai vijda'dai 


Considering a simple extension along one of the x;-axis we have 
dl = é;, the stretching turns into 


6; =4/1 + 27% — 1 (2.49) 


jj +1=1/142 (2.47) 
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Angular dilatation 


Let us consider the vectors dl and ds at a position p in the un- 
strained state which are deformed into vectors dl’ and d3’, respec- 
tively. The difference between the angle amid the deformed vectors 
and the unstrained vectors is called angular dilatation. For the sake 
of simplicity, suppose that dl = è, and ds = ē%. We define the 
angular dilatation the following difference 


T 
a 5 Pio (2.50) 


Tal 
da, 


P dal, T 
Figure 2.3: Angular dilatation. 


See figure 2.3. 
The scalar product of the strained vectors is 


dl’ - ds’ = dl'ds' cos yy (2.51) 


and the modulus of both strained vectors can be written by means 
of the preceding result for the linear dilatation 


dl! = (1 + ôı)dl = 1 + 6; (2.52) 
ds’ = (1 + ĝ2) ds = 1 + by (2.53) 


so that equation (2.51) becomes 


dl’ - ds’ = (1 + 61) (1 + 62) cos yho (2.54) 
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The left hand term of the latter expression can be written with 
the help of the strain tensor, so that by recalling equation (2.27) we 
have 

dl’ - ds’ = (ij + 275) dl'dsì = 2712 (2.55) 


Finally equation (2.54) turns into 
2712 = (1 + 61) (1 + ô2) cos p42 (2.56) 


By virtue of the identity sin w12 = cos £12, the angular dilatation 


becomes 
2712 


(1 + 61) (1 + 62) 


and naturally the above formula can be used to compute also the 
dilatations w23 and w31. 


sin W12 = 


(2.57) 


Area dilatation 


Vectors dl and ds at a position p in the unstrained state de- 
fine an area element dA which is deformed into the area element 
dA’ defined by the strained vectors dl’ and ds’. We define the area 
dilatation ratio the following coefficient 


dA' — dA 


We may suppose for simplicity that dl = è; and ds = é9. See 
figure 2.4. 


Tal 
da', 


P dx h T 


Figure 2.4: Area dilatation. 
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As well known, we have 
dA’ = |d]! x 3'| = dl'ds' sing), (2.59) 


and, recalling equations (2.52) and (2.53), the latter expression be- 
comes 


dA’ = dl'ds' sing, = (1+ 61) (1 + 62) sin yo (2.60) 


finally, through the geometrical relation cos w12 = sin 912 it is easy 
to reach the following expression for the finite area dilatation ratio 


a = (14+ 61) (1 + 62) cos w12 — 1 (2.61) 
that can be alternatively written as 
a= (1 + 61) (1 + 62) V 1— sin? w2 — 1 = 


= y (1+ 1) (1+ 62) - 492, (2.62) 


Volume dilatation 


We define the volume dilatation ratio the coefficient 


dy! — dV 
y= (2.63) 


As in the preceding cases, let us suppose that the initial un- 
strained volume is given by the following unit vectors 


dV = Bi x E2] - 3 = €123 = 1 (2.64) 
Thus, the volume dilatation turns into 
v=dV'-1 (2.65) 


For the strained volume we have 


dV' = [di] x d9] - dl, = (2.66) 
= dl dlydl; = (2.67) 
= (1 +61) (1 + 63) (1 + 63) (2.68) 


Finally, the volume dilatation ratio becomes 


y= (1 + 01) (1 + d3) (1 + d3) — 1 (2.69) 
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2.3.2 Infinitesimal deformations 


In the preceding section we have discussed the theory of finite 
deformations. Now, if all the components of the displacements u and 
the displacement gradient tensor u; j are very small we may neglect 
the squares and the product of these quantities in comparison with 
the first order derivatives themselves. So equation (2.99) becomes 


1 
Eij = 3 (Mig + Ugi) (2.70) 


where e denotes a symmetric second-order tensor named infinitest- 
mal strain tensor. 


Explicit compatibility equations 


Now we want to know if any state of given strain ¢;; yields a 
displacement field uj at every point p € V. To ensure that we have 
found equations (2.70) and to solve the differential equations system 
we discard the components of displacements u; as follows 


269 hk Wi she Ue ihk (2.71) 
2Ehk ij = Uh,kij + Uk, hij (2.72) 
—2Eih jk = —Ui, hjk Una (2.73) 
—2Ejk ih = —Uj, kih + Uk,jih (2.74) 


Summing equations (2.71) to (2.74) yields the necessary condi- 
tion to ensure the existence of the field ù. 


Eij,hk + Ehk,ij — Eih,jk — Ejk,ih = 0 (2.75) 


Infinitesimal stretching ratio 


When we work in the frame of linear deformations, i.e with the 
infinitesimal strain tensor, the stretching ratio is given by the first 
order approximation of the ratio in (2.49), namely 


2 i 
Site EST i lac (2.76) 
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Infinitesimal angular dilatation 


We invoke again the first order approximation of expression 
(2.57), so that, by replacing the finite strain tensor with the in- 
finitesimal strain tensor and by using the latter result for the stretch- 
ing ratio, the angular dilatation assumes the following expression 


a 2€19 da 
V1+2e11v1+ 2629 


With the proper subscripts shifting we can also write the angular 


2€19 (2.77) 


W12 


dilatations w23 and w31. 


Infinitesimal area dilatation 


Recalling equation (2.60), that is 
a = (1 + 61) (1 + 62) cos p42 


the infinitesimal area dilatation is obtained, as usual, by neglecting 
the second order terms, so that 


a ~ di + do = €11 + E22 (2.78) 


Infinitesimal volume dilatation 


From equation (2.69), the first order approximation leads to the 
following expression for the infinitesimal volumetric dilatation ratio 


v ~ 6, + 69+ 63 ~ £11 + E22 + €33 = E4369 (2.79) 


2.3.3 Deformation and rigid body motion 


It is rather intuitive to understand that the motion of a flexible 
body can be made up of rigid translations and rotations as well as 
deformations. To see that from a mathematical point of view, con- 
sider the displacement field & in a point p, as defined in (2.34), being 
defined by the first order approximation from the displacement tig 
on po. 

Uj = uoj + Uj 4X5 
where it is clear that the translational component of the motion is 
wholly yielded by uoj. Consequently the remaining part must store 
the deformation and rigid rotation components. 
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By observing that the gradient of ù may be also written as 
follows 


1 1 
Uji = 5 (Uji + uij) + 5 (Usa — ig) (2.80) 


the displacement field becomes 


1 1 
Uj = Uoj + 5 (Uji + uij) dati + 5 (Uji — wig) dx; = 


= uoj + Ejidxi + wjidx; (2.81) 


where it has been put wji = 3 (uji = u; g): 
So, through the latter expression, the splitting of the displace- 
ment field u appears clear: 


e uo,: pure translation; 
o cj = 3 (uji + uij): pure deformation; 
e wji = 3 (uji — ui j): rigid body rotation. 


In order to give a physical meaning to the operator curl intro- 
duced by equation (1.133) on page 26, we can observe that 


curl ù = EkijUj iek = 
1 1 _ 
= g kij (Uji + uij) + (Se (uji — u3)] Ek = 
—_—___y—_ — 


= EkijWjiCk = WkEk (2.82) 


Using the identity (1.105) it is possible to provet that the skew- 
symmetric component of a tensor is given by 


1 
Wii = g CkigWk (2.83) 


thus the rigid rotation turns into 


1 1 = 
Wjidx; = z Chijkdri = 3° x dl (2.84) 


where © = curl ū and dl = dx,&;. 


4 
EklpWk = EkiperijWji = (b1bpj3 — d15Opi) Wji = Wp). 
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2.4 Shell continuum 


We define a shell-shaped region modeled on a surface Q and 
with thickness 2e as a continuous medium G (e) embedded in the 
Euclidean space E each point of which is determined through a 
coordinate system {x®, €} : G (e) — IR3. Therefore, given p* € G (e) 
it is defined by its position p normally projected on Q - by using 
the surface coordinate system introduced in (1.153) - and by the 
normal coordinate € taken along the unit normal vector n. In fact 
we have 


pi (x° (p) ,€ (p)) (2.85) 


The basis induced by the coordinate system {x®, €} is 104 n}. 

It is worthwhile pointing out that mechanics of shells - by virtue 
of such above statements - is traced back to the theory of surfaces, 
in fact vectors and tensors fields belonging to TgE will always be 
split into the parallel and normal components. 

Note also that the symbol « denotes quantities belonging to the 
shell continuum. 


2.4.1 General assumptions 


The shell theory here introduced is based on the following hy- 
potheses 


Hypothesis 1 The shell is sufficiently thin, so that 


2e ; 

Ta <1 L = min { Rmin, Lmin} (2.86) 
where Rmin and Lmin are the minimum radius and a typical dimen- 
sion of the shell structure, respectively. 


Hypothesis 2 (LINEAR THEORY) Displacements are infinitesimally 
small such that their products can be neglected in the kinematic ez- 
pressions. This assumption allows us to write the equilibrium equa- 
tions in the unstrained shell configuration. 


Hypothesis 3 The material filaments along the coordinate € re- 
main straight throughout the deformation and no length change is 
allowed. Namely, the distance between p* € G(e) and the surface Q 
is unaltered 

E = const. (2.87) 
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Hypothesis 4 (KircHHorr-Love THEORY) The line elements ini- 
tially normal to the shell’s mid-surface remain normal to it during 
the deformation. 


g (dogs Pa) = 0 (2.88) 


where the subscript d is denotes quantities related to deformed con- 
figuration. 

Note that the last hypothesis is nothing but the extension to a 
two-dimensional model of the Bernoulli theory for beams. 


2.4.2 Strain tensor 


A generic point p* € G(e) is determined by the vector 7* referred 
to the global Cartesian axes, so that 


rm=T+en (2.89) 


where £ € (—e,€). See figure 2.5. 

Let us suppose now that a quasi-static motion produces a de- 
formed shell configuration points of which are univocally deter- 
mined by the vector 

Ty = Ta + ana (2.90) 


where £q € (—€, €). 
The displacement field is obtained by subtracting equations 
(2.89) and (2.90), so that 


Ti — T* =T- T +E- N) (2.91) 


where we have made use of hypothesis 3. Equation (2.91) allows us 
to define the positional field as a function of two vector fields 


ö=rq4 -7 VETQE (2.92) 
W=fig-h weTQ (2.93) 

To obtain the strain tensor no more theoretical concepts are 
required. We already know the definition and we just need to com- 


pute the metric tensors associated to the coordinate systems in the 
strained and the original configurations, so we have 


DES ( Yab Ya3 ) 
5 Y3a V33 
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Figure 2.5: Two dimensional sketch of the displacement field for 
Kirchhoff-Love shells. 


where 
1 
Yas = 5 (Iapa — Tao) 2a 
1 
Yas = Ysa = 5 (Gas. — Jas) (2.95) 
1 
733 = 5 (Ma Ra — 1) = 0 (2.96) 
According to equation (1.79) we have 
Iapa = Dha a (2.97) 
and oe 
Jag = dd (2.98) 


where, recalling equation (1.76), we can write 
roa = 5 [8a i — 85-3] = 


Ons "E EVaiia) ; (0g, Az EV ghia) | + 


NIH wpe Ni ple NIE 


da + EV gn) - (Og + EVair)] = 


On: x Og, + da È EV ena ag 03, i EVaña an EV aña > V ghia] 


Ja Og + On: EV ei + Ig Van + E Van: Ven] (2.99) 
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where we realize that the tensor yag can be split in three parts as 
follows 


Yap = Aap + Ewag + E pag (2.100) 
We define the stretching strain tensor as 

Le = a 1 

Qag = 5 [Boa + a — da + dg] = 5 apa — 908) (2.101) 
next, the first bending strain tensor as 
tees z _ _ 
Wag = 5 [Dax . Vga + dg, -Vañd — da * Ven = Og . Vafi] 

(2.102) 

and the second bending strain tensor as 

1 

Paps = 9 [Vana i Vana = Ven Ven (2.103) 


Considering now that the displacements are small enough to be 
negligible the second order terms 


Vav:Vgu=0 
Vad:Vgw=0 


and recalling equations (2.92) and (2.93), the stretching and the 
bending strain tensors become, respectively 


bo : 1 
aap = 5 (daVad daVad) = 5 (vais + Ugla + 2ÉLag) (2.104) 


1, 7 
Wap = 5 (da Va + da + Va) + 


+ ; (Vat: Vgn+Vg0: Van) = 

= ; (ways + Wala + Lye + vi Lao) i 

+ ; G (14Ls + L}Lra)) (2.105) 
Pos = ; (vw) Lye + Wha) (2.106) 


where we have put 
Va0= (4612) + (h)a IM 


v = v +T (2.108) 


a 
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and 
Vath = w by- wLayfi (2.109) 
wa = wh, HWT} (2.110) 

and 
Va Ve = L718, - (wider. = w L gait) = Eva, (2.111) 


Finally, the strain tensor assumes the following form 


1 
Yap = z (vais + UBla + 20 Las) + 


1 
Tzs (ways + Wola + V Lyo + vf Ener} + 


la 
106 (10 Y 
+4 (e (quei) 
1 
+le (wi Lye + wLa) (2.112) 


The stretching strain tensor does not depend on the thickness, in 
fact it describes the deformation of the mid-surface Q. The bending 
strain tensors describe the deformation along the thickness. 

The transversal components of the strain are 


YBa = Ya3 = 2 (Pa ; da) =n: ða) = da ni v’ Lay twa (2.113) 


Kirchhoff-Love strain theory 


If we take into account the Kirchhoff-Love hypothesis, see hy- 
pothesis 4, we have 


Dag a =0 > (A+ 0): (Og + Vað) =0> (2.114) 
Dda + Vað = 0 > Wa = V Lay — V5, (2.115) 


and we observe that the variables reduce just to the field v. Thus, 
the strain tensor turns into 


1 
aag = 5 (Yala + ala + 20ÉLag) (2.116) 


Wag = v” Lyg + vigLya +v” Lya Z ba + ve Li Lag (2.117) 


la 


56 ANALYSIS OF STRAIN 


pag = € (of, LoL) SR E vf Lh) J 


2 ô ô 
yé (v Linh ee aE wf 402) (2.118) 
In the linear theory the second bending strain tensor can be 
neglected because € is very small and its square makes the contri- 
bution of pag insignificant. 
Finally, we have 


733 = Ya3 = Y3a = 0 (2.119) 


Consider now a Cartesian coordinate system where all the Chri- 
stoffel symbols vanish, we immediately realize the well known ex- 
pression of the strain tensor for bending plates 


y si Vag + Uga — 20° (2.120) 
aß 2 a, Ba Neto) ` 


Chapter 3 
Analysis of stress 


This chapter presents the classical stress analysis of a three-dimen- 
sional continuum subjected to both body and surface forces. It begins with 
the notions of stress vector and stress tensor which bring to enunciate 
the famous Cauchy’s theorem, then the static equilibrium equations will 
be derived. 

Next, the graphical representation through Mohr’s circle and the prin- 
cipal directions associated with the state of stress will be also analyzed. 

Curvilinear coordinate systems will be introduced only in the last sec- 
tion, where the analysis of stress for shell continuums will be shortly in- 


troduced. 


3.1 Body and surface forces 


Let V be the configuration of the continuous medium. We sup- 
pose that V is bounded by the closed surface S. Consider a small 
region AV subset of V and a small surface element AS of S. To 
analyze the forces acting on the volume element AY it is necessary 
to account for two types of forces: 


Body forces (or volume forces). These are the forces which are 
proportional to the mass contained in the volume element 


AY. 


Surfaces forces. These are the forces being measured as force per 
unit area of surface AS on which they act. 


A good example of body forces is gravity: pgAV - where p is 
the density of the continuum and g is the gravitational acceleration 
- or inertia. 

Examples of surface forces are: pressure and tension tn (p,7) 
(discussed in depth later on), which two parts of a continuum mu- 
tually exchange. 


Claudio Borri, Michele Betti, Enzo Marino (a cura di), Lectures on Solid Mechanics, ISBN 978-88-8453- 
854-3 (online), ISBN 978-88-8453-853-6 (print), © 2008 Firenze University Press 
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Figure 3.1: Body and surface forces. 


If we imagine that the continuous medium V, equilibrium of 
which we are searching, is a subset of a bigger imaginary continuous 
body, then the tensions exchanged between these two portions can 
be assumed as external force loads. 

In order to write the equations of equilibrium we consider both 
body forces b = b’é; and surface forces fn. See figure 3.1. 

The body forces also produce a resultant moment M = M°*&;, 
where 


i= | (Fx?) dV (3.1) 


3.2 State of stress 


Let V be the configuration of a continuous medium, whose points 
are referred to a rectangular coordinate system 


g: EIR: prg((p-0),&) (3.2) 


where p and o are points of E and {é;} is the unit normal basis of 
Suppose on the body V surface and body forces, e.g. b, Sa Mx, f 
act in such a way to assure the equilibrium state. See figure 3.2. Due 
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to these forces throughout the body internal reactions are activated 
between the material points. 

To understand the stress condition created at generic point p 
within the body V, we suppose to cut the continuous medium by 
means of a generic plane mn, so that two portions Vj and V2 are 
produced. 


Tı 


Figure 3.2: Body V being in an equilibrium state. 


After splitting, the portions of the body on the left and on the 
right side of the section plane 7, lose their equilibrium state. In 
fact, before parting, both V and Vz were in equilibrium due to the 
mutual forces exchanged through the plane. 

Cauchy enunciated the principle that, within a body, the forces 
that an enclosed volume imposes on the remainder of the body must 
be in equilibrium with the forces from the remainder of the body 
itself. 

We denote by AA, the small area surrounding p and by AS, 
and AM, the force and the couple resultants in p stemmed from 
the internal force distribution acting through A An. See figure 3.3. 

Cauchy’s principles implies the following limits 


1. lima 4,0 ee = tn (p, n) = —tn (p, _n) 


2. limAA,-0 aia =0 
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-ta(p,ù) 


AAA 


` ta(p,ñ 
fs, (pù) 


Figure 3.3: Splitting of the continuous media V. 


The vector tn (p, ©) is called the Cauchy stress vector and rep- 
resents the surface force per unit of area acting at point p. The 
second limit assures that the entire state of stress for a fixed point 
p is only defined by the forces, that is the couples are infinitesimal 
in comparison with them. 


It’s important to observe that t, is a linear mapping defined as 
follows 


in:ExE+ E (3.3) 


so that 


in (p) € L(E, E) (3.4) 


and we recognize tn (p) to be an endomorphism which is associated 
to a tensor belonging to E* ® E. 


In the following paragraphs this tensor will be thoroughly ana- 
lyzed. 
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3.2.1 Stress vector components 


Let n be the unit normal vector of the surface AA, located at 
p. We can write the components of stress vector! tn (p, ñ) as follows 


în (pn) = th (Pm) & (3.5) 
so that the normal components of tn (p) can be easily written as 
tn (P, 0) = Ain (P, 2) = th, (p, 7) ni (3.6) 


Let us observe that the stress vector, which represents the entire 
state of stress at p, is completely known if the three coordinate 
components t} (p,m) are known. 


3.2.2 Stress tensor 


Now we want to show that the state of stress at any point of 
the continuum is entirely characterized specifying a linear mapping, 
i.e. endomorphism, represented by the nine quantities called com- 
ponents of stress tensor. 

As usual, p is a point in the medium and tn (p, ©) is the stress 
vector acting on the surface element passing for p with the unit 
normal n. Imagine to have four planar elements, three of which 
are parallel to the coordinate planes, the fourth one is supposed 
passing normal to n, at a very small distance to p. We obtain a 
small tetrahedron. See figure 3.4 

We shall denote by #;, with i = 1, 2,3, the stresses vector? acting 
on the planar surface element orthogonal to the coordinate curves 
xi, namely t; = t; (p, €i). Evidently, every stress vector can be writ- 
ten by its components in the following way 


=e; 4,7 =1,2;3 (3.7) 


where l 
i = €;. ti (3.8) 


The forces acting on the tetrahedron are both surface and body 
forces: 


‘It must be noted that in general the stress vector tn (p, ñ) is not in the 
direction of n. 
?Rigorously t; should be write as te, (p). 
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TL3A 


LI 


Figure 3.4: Stress vectors: the sketch of Cauchy’s theorem. 


- Body forces: bdV 
- Surface forces: —t; (p,—€;) dAj+tn (p, n) dAn, with i = 1, 2,3. 


thus the translational equilibrium of the tetrahedron can be readily 
written as 


— t; (p, —&) dA; + tn (p, 7) dAn + bdV = 0 (3.9) 


that taking into account that dA; = dAnni, indeed we have nê = 
n+ €; = COS (7, €), the above expression turns into 


_ = 1 
—t; (p, —é;) dAnn; + tn (p, n) dAn + 3 P9hdAn =0> (3.10) 


i 1 
—ti (p, Ei) ni + tn (p, 7) + 3pgh =0 (3.11) 


and, for h approaching zero, i.e. the infinitesimal volume surround- 
ing p, the equilibrium becomes 


tn (p, 7) = ti (p, —€i) ni (3.12) 
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Figure 3.5: Stress tensor components. 


Equation (3.12) represents Cauchy?’s theorem which states that 
the stress state tn (p,m) can be completely determined by the stress 
vectors t; (p, —&;), acting on the face with outward unit normal vec- 
tor —€;, where © is considered known. This result also proves that 
we are really dealing with a tensor as introduced by the endomor- 
phism (3.4). 

It will be convenient to use the customary notation, so that 
equation (3.12) may be rewritten in components as follows 


ti (p, n) a; = tl (p, —ēi) me; (3.13) 
from which the stress tensor o is defined as 
t (p,n) = TijNi (3.14) 


The tensor oj; is called the stress tensor, it com- 
pletely defines the state of stress at point p and repre- 


3 Augustin Louis Cauchy (August 21, 1789 - May 23, 1857) was a French 
mathematician. 


Source: http://www-history.mcs.st-andrews.ac.uk/Mathematicians/Cauchy.html. 
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sents the component of the vector t; working in direction 
of x. 


See figure 3.5. 
We can summarize saying that 


o:E>L(E,E), pro(peE*®E (3.15) 
so 
o (p) (7) = tn (p, 2) (3.16) 
that in components, (3.16), becomes 
a (3.17) 


Figure 3.6: Stress tensor components acting on an infinitesimal 
volume element. 


We remind the reader again that the lower and upper indices, 
in this context, can mutually be exchanged; moreover, they can be 
placed both upper and lower. So, generally, we shall also write 


tnj = Ohj"h (3.18) 


3.3 Equations of equilibrium 


3.3.1 Translational equilibrium 


With respect to the body V, bounded by the closed surface S, 
the condition of equilibrium requires that 


[bare [tras =o (3.19) 
V S 
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Making use of (3.16), equation (3.19) becomes 


[rare fo (©) ds =0 (3.20) 
y S 


The divergence theorem can be applied to integral (3.20), so 
that 


[iw f divo (o) dv= | (B+ divo (p) dV=0 (3.21) 


Since the region of integration V is arbitrary, i.e. each part of 
the medium is in equilibrium, integral (3.21) vanishes, thus, at every 
point of V we have 

divo +b=0 (3.22) 


that in components becomes 
Cig + bj =0 (3.23) 


3.3.2 Rotational equilibrium 


As well as the translational equilibrium, we require that the 
moments acting on the body are also in equilibrium, so 


fd a+ | (xt) dS=0 (3.24) 
y 


S 


which in components, by virtue of the skew-symmetric tensor e€, 
becomes 


I (r'b! ekijēr) dV +f (r'é; X TjhNne;) dS =0=> 
V S 

f (r'bierijēr) w+ f (r'ojnnn€Kijek) dS = 0 

y S 


With the aid of the divergence theorem, for the k-th component 
we can write 


esi f (r'b’) dV + €x 5; | (r°0;n) h dY = 0 > 


Ekij i (PD + rinojn +Tiojnn) dV =0 
v 
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Recalling equation (3.23), and that rin = din and since the 
volume V is arbitrary, the rotational equilibrium produces 


EkijTij 7 0 (3.25) 


Therefore equation (3.25) imposes the symmetry of the compo- 
nents of the stress tensor: 


Tij = O ji (3.26) 


The symmetry of the stress tensor can be also seen considering 
the volume element taken in shape of a rectangular parallelepiped, 
with faces parallel to the coordinate planes and with stress vector 
ti acting on the face perpendicular to the x-axis. Denoting the co- 
ordinates {x1, £2, £3} with {x,y,z} - as often happens in literature 
- for the (y, z)-plane the rotational equilibrium becomes 


(0y2dxdy) dz = (c-ydxdy) dz > Oyz = Czy (3.27) 


See figure 3.7. 
If we write the equilibrium for all planes, we obtain again the 
result in (3.26). 


zz 


Figure 3.7: Plane (y, z). Components of the stress tensor acting on 
the volume element. 


3.3.3 Boundary equations 


Let f be the external force acting on the surface Sọ and w the 
displacement field imposed on the remaining portion Su, so that 
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S = Ss U Su. Each point of V lying on the boundary S must satisfy 


the equilibrium and kinematic conditions as follows 


indS = fdS, Vp € Sc (3.28) 
=u, Vp © Sy (3.29) 


&I 


that in components is 


OhjNh = Íi» Vp E Si (3.30) 


3.4 Principal stresses and principal directions 


Let us consider now the sheaf of planes passing through p € V. 
Among the infinite planes there are some for which all the stress 
components vanish except the normal one. These planes are said 
principal planes and their normal directions are said principal di- 
rections. Hence, if n is a principal directions, by definition, we have 
at p 


tn = oñ, o € IR (3.32) 
To find the three principal stresses we impose 
tn = Ojnnne; = oh > CihNhEi = 0Né; (3.33) 


so that 
Cin — CÒ;nnn = 0 > (Cin = odin) ny = 0 (3.34) 


Expression (3.34) is a set of three homogeneous equations in the 
unknown direction n. The solution is nonvanishing if, and only if, 
the determinant of the coefficients matrix is equal to zero; that is 


[oij = odi; | =0 (3.35) 


Solving the determinant above we obtain the cubic equation 
called secular equation in the unknown stresses o 


o — ho? = ho = I =0 (3.36) 


The (3.35), (or (3.36)) has three real roots 07, orr, 0117, which 
are called principal stresses. If o in equation (3.34) is replaced by 
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any one of these eigenvalues, the resulting set of equations may 
be solved for the corresponding direction n. These directions, ñy, 
Arr, Dir are called principal directions. The planes normal to the 
principal directions are termed principal planes of stresses. In other 
words we say that along the principal planes of stresses there is no 
shearing stress. 

Generally there are only three mutually orthogonal principals 
directions. 

The three scalars in (3.36) are invariants as regards to the co- 
ordinate system. They are 


I, = tro 
1 
Io = 2 (Giiojj = Oi; 0%) 


Ís = det (Sij) 


These invariants are physically very important, they in fact al- 
low us to characterize the stress state as follows 


if I3=0 : triaxial state of stress 
if J3=0 and Ig 70 : biaxial stete of stress 
if I= Ñ= 0 and I #0 : axial stete of stress 


Now we want to point out that the principal stresses found solv- 
ing equation (3.35) represent the maximum and minimum stress. To 
see this we make use of the Lagrange multipliers method in order 
to find the extremes of a function of several variables subjected to 
one or more constraints. In this case recalling formulae (1.33) and 
(1.34) we can write the stress tensor in a generic unknown coordi- 
nate system rotated with respect the initial system as follows 


C= a opa (3.37) 
and we also impose the constraint on the unknown matrices a and 


a' such as they are effectively two orthogonal transformations, i.e. 
the condition (1.21). Hence we have 


L (a, à) = a opa -À (ajal) — õi) (3.38) 
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and the stationary conditions are 


OL 

dal = ofa — Val} = (ot — xd) all = (3.39) 
OL i th _ gi 

Dak = apa; — 6; =0 (3.40) 


The first equation, (3.39), yields the following condition 
lof? — AS}| =0 (3.41) 


that is exactly the condition (3.35), hence we can derive that given 
a generic state of stress o;;, the principal stresses 07, Orr, Orry are 
extrem values. 


3.4.1 Normaland tangential components of the stress 
vector 


The last equations of the previous section enable us to know the 
components of the stress vector for every direction we want. Let n 
be the unit normal vector and v the unit tangent vector. It follows 
that the normal and tangent components of the stress vector, o 
and 7, respectively, are readily computed through the usual scalar 
product as follows 


=: (p, n) n= tjn (p, n) Nj = OijNinj (3.42) 
T = tn (pP, ñ) -D = tjn (P, 0) Vj = dg 3.43) 
n 


Figure 3.8: Normal and tangential components of the stress vector. 
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From figure 3.8 on the preceding page it is also clear that the 
square tangent component of the stress vector can be written as 
follows 

rela (3.44) 


3.4.2 Mohr’s circles 


Two dimensional state of stress 


An important graphical interpretation of the above results is due 
to O. Mohrt. Following [10], [9] and [12] let us begin considering the 
above relations (3.42) and (3.43) for a two dimensional problem, so 


that 
| % m 
Gij ( Pa ey. ) (3.45) 


The unit vectors n and P, with respect to figure 3.9, have the 
following components 


e (E) e CE) 


Hence, the normal and tangent components of tn (p, ñ) are 


o = 011 cos” Y + 099 Sin? 02019 sin Y cosy (3.46) 
T = —011 cos Y SIN Y + 099 SÌN Y cos Y + 012 cos? P — 021 sin? p 
(3.47) 


that through some trigonometric manipulations” turn respectively 


‘Christian Otto Mohr October 8, 1835 - October 2, 1918 was a German civil 
engineer. 


Source:http://en.wikipedia.org/wiki/Otto-Mohr. 
5In particular these two identities have been used: 


i) 2sinycosy = sin 29, 


ii) cos? y — sin? y = cos 2y. 
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T24 


Figure 3.9: Normal and tangential components of the stress vector 
in two dimensions. 
into 


1 1 
o= rae cos? pt 3011 (1 — sin? p) + 


1 1 
+ -092 sin? P + =022 (1 — cos? p) + 012 Sin 2Y = 


2 2 
1 1 | 
= (011 + 022) + 3 (011 — 022) cos 2y + 012 sin 2y (3.48) 
T = —0 COS Y sin Y + 099 sin Y COS Y + O12 cos? P — 021 sin? g= 
1 
si; (011 — 022) sin 2p + 019 cos 2p (3.49) 


Next, by squaring both terms of the latter equations and sum- 
ming term by term, the variable 2y disappears, hence 


1 7 1 È 
G — 5 (011 + om) + 7 = (5 (011 — 0) + o? (3.50) 


If we represent the above equation in a two dimensional Carte- 
sian system with o and 7 as abscissa and ordinate, respectively, we 
realize it represents the equation of a circle in the form (x — ro) + 
(y — yc)? = R2 where 


1 
tes (011 + 022) (3.51) 
yo =0 (3.52) 


are the coordinates of the center and 


2 
R= NE (011 - on)) + 02, (3.53) 
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is the radius of the circle. This circle is known as Mohr’s circle and it 
represents all the possible states of stress in p. Namely, there exists 
a one-to-one relationship between each state of stress tn (p, ñ), i.e. 
o and 7, and points belonging to the circle. To show that, let us 
assume y is the angle between the xj-axis and the stress vector 
tn, as figure 3.10 depicts. To find the correspondence between the 
stress state and the circle let us observe that y defines a principal 
direction so, by definition 7 = 0, and we have 

tan 2y = MiA (3.54) 

O11 — 022 

and with the aid of picture 3.10 in the above equation we recognize 
that P,P* = 2042 and P2P* = 011—022. Consequently the following 
expressions hold 


1 
r cos 2y = 5 (011 — 022) (3:55) 
rsin 2y = 012 (3.56) 
(3.57) 


that substituted into (3.48) and (3.49) and making use of some 
trigonometric identities® yield, respectively 


1 
Co (011 + 022) + rcos2(y— £) (3.58) 
T=rsin2(Y- p) (3.59) 


Thus, given a generic plane oriented as y equations (3.58) and 
(3.59) are a parametric representation of a circle and so a one-to- 
one relationship between the state of stress and the Mohr’s circle is 
established. See figure 3.10. 

We define P* = (011, —012) as the pole of the circle. The line 
passing through P* having inclination y with respect to the vertical 


SIn particular these identities have been used: 


i) cos2ycos 2p = 3 (cos 2 (y + p) + cos 2 (y — ¥)), 
ii) sin2ysin 2% = 3 (cos 2 (y — £) — cos 2 (y + ¢)), 


iii) sin 2% cos 2% = 3 (sin 2 (y + ọ) + sin 2 (y — ¢)). 
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line joining P) and P* intersects the circle in P, and the angle 
PCS, is right 2( — 7), so the coordinates of the point PL, in 
the (o,7)—plane, are just those expressed by equations (3.58) and 
(3.59). 

Thus, we have proved that, given a stress vector orientated as 
y, once Mohr’s circle is known, the normal and tangent components 
of a stress vector can be graphically found provided the inclination 
vy in known. 

On the other hand, if the normal and tangent stresses are known, 
Mohr’s circle enables us to find directly the principal direction. In 
fact point S has coordinates o = OC + R and 7 = 0, so that the 
line P*S) defines the angle y that fixes the principal direction. See 
figure 3.10. 


Th Oj 1-922 


Figure 3.10: Normal and tangential components of the stress vector 
for in two dimensions. 


Two other relevant features on Mohr’s circle are those for which 
the tangent component of tn (p, ñ) is maximum. These directions 
can be found through the same procedure. Indeed the lines P*T; 
and P*T represent the directions along which the stress vector has 


maximum shear component. See figure 3.10 and 3.11. Analytically 
these maximum and minimum values are 


1 2 
Tmax = Vou — 022) + 403, (3.60) 


1 
Tain = -5y (0u — 0o22) + 40, (3.61) 
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and the directions can be computed putting zero the following 
derivative 


dr T T 
— =0=>cos2(7-09)=->7=9+- 3.62 
ay SO, (3.62) 


Let us summarize now the key items to draw and use the Mohr’s 
circle when a plane state of stress o;;, with i, j = 1,2 is known with 
respect to a generic system {x}, £2}. See figure 3.11. 


1. Compute the radius R and the abscissa of the center C of the 
circle, equations (3.51) and (3.53); 


2. Identify the pole P*; 
3. Identify the principal direction drawing a line from P* to both 
Sı and S2. The inclination of the latter defines the principal 


directions; 


4. Compute the principal stresses 07 and 07; at the extreme 
points Sı and So, respectively; 


5. Compute the maximum and minimum shear stresses Tmin and 


Tmax- 


ZI 
> 


Figure 3.11: Graphical determination of principal directions. 
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Three dimensional state of stress 


Consider again the state of stress in p referenced to the principal 
axes and let the principal stresses be ordered according to or > 
Cir > ozzy. Assume the three principal stresses are known, so that, 
in accordance with equations (3.42) and (3.44), we write 

o = omi F orn? F onmi 
2 2 2 
T +o? = (orm) + (07,n2) + (o771N3) 


and being n. n = n? + n? +n? = 1, by solving for the directions ni, 
we obtain 


2 7? + (o — orr) (o — arrr) 
nı = 


_ T? +(0o-— orr) (o — or) 
> fs) (3.64) 
2 
“a “i (3.65) 


(arrr — 01) (orr — O11) 


In the above equations o7, Cr, Ozzy are known; o and T are 
functions of ni. 

In order to interpret these equations graphically we note that in 
equation (3.63) or — orr > 0 and 07 — Grrr > 0, and n? is positive. 
Therefore 


(0 — orr) (0 — orr) +77 > 0 (3.66) 
When the equality sign holds, this equation may be rewritten 
as 
1 a 1 A 
aa (orr + ou) | +72 = $ (orr — ar)| (3.67) 


which is the equation of a circle in the (0, 7)-plane, where we assume 
o as abscissa and 7 as ordinate. The circle in figure 3.12 is termed 
Cı and has the center in 3 (or + arrr) on the o axis, and radius 
3 (Orr — ort): 

Examining equation (3.64) we observe that orr — crrr > 0 and 
ort — 07 < 0, so we have 


(o — orr) lo — 07) +r? <0 (3.68) 
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Orr Orrr 
2 


Figure 3.12: Mohr’s circles. 


The boundary of the area of this equation, i.e. in the case of 
equality sign, defines a circle as before in the (0, 7)-plane 


2 


o ; (or + ou)| i +r = È (or — ou) | (3.69) 


named C3. See figure 3.12. 
The same procedure allows us to obtain from equation (3.65) 
the circle C3, indeed, we have the following condition 


(o — or) (0-071) +7° >0 (3.70) 


that at the boundary yields 


o 3 ; Gi arr) See È Ce arr) i (3.71) 


Finally, from inequalities (3.66), (3.68), (3.70), it follows that 
admissible values of o and 7 lie in the shaded region of figure 3.12 
bounded by the circles Cy, Co, C3. The value Tmax and Gmax can be 
readily provided from figure 3.12, so that 


(or — orr) (3.72) 


(or + 0111) (3.73) 


and as a consequence, the surface elements supporting these stresses 
are found replacing the above values into equations (3.63) to (3.65). 
For further details the reader is referred to [1]. 
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3.5 Stress quadric of Cauchy 


Consider an elements of area dA with a normal vector n. As 
previously stated, the stress vector tn (p, ñ) can be decomposed into 
a normal component o and a tangential component 7. 

Let us introduce now a local system of axes €; with origin in p 
equipped with the unit normal basis {é;}. See figure 3.13. 


Figure 3.13: Stress quadratic of Cauchy. 


Let # be the vector taken along © joining p with a generic point 
p', namely (p' — p) = rn. This vector can be equivalently expressed 
by the following expressions 


SI 


=rħñ (3.74) 
The first equation provides the j-th component of 7 as follows 
&=T-€= rn; (3.76) 


that, replaced into the expression of the normal component of the 
stress vector, yields the following relation 


rio = 0ijté; (3.77) 


We recognize equation (3.77) as a quadric form’. 


"We remind that any quadric form F can be expressed as 
F (0) = MijVvivj (3.78) 


where ọ = (v1, v2, v3)” is a vector expressed with respect to the chosen basis, 
and M;; is a certain symmetric matrix that depends on F and on the basis. 
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So we restrict the coordinates of é; by requiring the end point 
p' of F to lie on the quadric surface 


F (£1, €2, €3) = oij€iéy = +k? (3.79) 


where k is an arbitrary real constant and where the sign is chosen 
in such a way to make the surface real. As a result we have 


k2 


ae 


oO=a 


(3.80) 


Since r? is a positive quantity, k? will be taken with the posi- 


tive sign whenever the normal component o is a tension and with 
negative sign when it represents compression. 

Next, deriving equation (3.79) and by using equation (3.76), we 
obtain 


OF : 
BE a rtn (p, 7) (3.81) 


which allows us to realize that the quadric form (3.79) has some 
properties of a potential function, indeed the partial derivative of 
F with respect the i-th coordinate gives, except for the constant r, 
the force component (i.e. the component of the stress vector) right 
in the i-th direction. 

Furthermore we observe that the above derivatives, equation 
(3.81), denote the direction of the normal © to the plane tangent to 
the quadric surface (3.79) at point p’, so that the right-hand term 
of equation (3.81) just establishes the stress vector tn (p, n) is also 
normal to this tangent plane. 

The above results have been directly taken from [1], to which 
the reader is referred for any further detail. 


3.6 Stress-deviator and spherical components of the 
stress tensor 


Every state of stress oj; may be decomposed into a spherical 
portion and into a portion sij known as stress—deviator by the fol- 
lowing equation 

Oy = OM oi; + Sij (3.82) 


where om = Loi is the arithmetic mean of the normal stress, i.e. 
spherical stress component (or hydrostatic stress). Equation (3.82) 
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may be solved for sij 
Sij = Cig — OMbij (3.83) 


where the latter components are termed stress—deviations. 
Namely, 


011 OM 012 013 
Sij = 021 022 — OM 023 
031 032 033-0m 


It is possible to prove that both the stress tensor o and the devi- 
ator tensor s have the same principal directions. The characteristic 
equation for the deviator is 


s + Jos + J3 =0 (3.84) 


where the deviator invariants are 


1 
Jo = Tg iSi 
J3 = det Sij 


3.7 Stress in shell continuums 


3.7.1 Shifters 


Before reasoning upon the stress state characterizing a shell con- 
tinuum it is worth introducing some geometrical relations linking 
points belonging to the mid-surface Q with corresponding points 
belonging to the shell thought as a three-dimensional continuum. 

Therefore, let us recall the relation already met to compute the 
components of the metric tensor Ge as see equation (2.99) on page 
53, between the basis in p* € G(e) and the basis in p projection of 
p* on Q along the normal coordinate curve €. So we have 


ð% = ða + ELÉ Ög (3.85) 
w= (3.86) 


which in a short notation assumes the following form 


dx = SPO, (3.87) 
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Hence, with respect to the basis associated to the coordinate 
system {x®, E} the tensor S has the following components 


1+€L1 €L2 0 
S? = ELL 1+€53 0 
0 0 1 


Therefore, the superficial part of S can be expressed by the 
following tensor product 


st = d 9 or (3.88) 
so that 
St (ðs) = (d7 @ 4) (dg) = 95 (3.89) 
In the same way we define F' as follows 
Fi =0,@d"7 (3.90) 
so that 
Fİ (2”) = (& & d7) (2”) = (3.91) 


Tensors SÌ and FÍ are called shifter tensors. 


3.7.2 Contraction of surface forces 


Consider now a curve c : IR — Q representing the intersection 
of the surface Q normal to Q which splits the shell continuum G(e) 
into two portions. 

Let D € TQ be the unit normal vector applied in p outward 
pointing from c and let 1 € TQ be the unit vector tangent to c 
applied in the same point. Then the three unit vectors {p, ln} 
form a local basis in p. A similar triplet of vectors can be defined in 
p* as {pf IE, n}. Note that the symbol « denotes as usual quantities 
belonging to the shell thickness. See figure 3.14. 

In order to ensure the equilibrium condition, the portion of the 
shell included by Qe must exert on the remaining part of the contin- 
uum a tension such as for each point p* is entirely described by the 
stress vector t*. Moreover the stress vector t* can be equivalently 
expressed by Cauchy stress tensor as follows 


È (pò, v*) — g* (p*) vas (3.92) 
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Figure 3.14: Local bases in G(e) and in Q. 


where o* is the contravariant form of the stress tensor defined in 
p*. For the sake of brevity hereafter o* (p*) will be denoted simply 
by o. 

Now our goal is to establish a relation between the stress state 
distributed along the surface Qe and the stress state along the 
boundary of the mid-surface of the shell. This can be done by 
means of a reduction, i.e. a contraction, of the stress per unit area 
to a stress per unit line. 

Therefore, let us define two vector fields n and m such as 


f slo id f P(p*, v*)dA* (3.93) 
f m(p, v)di = f (pt — p) x Fp", v") jd A" (3.94) 


e 


Equalities (3.93) and (3.94) guarantee that the stress system n 
and m is statically equivalent to the stress system #* along the fiber 
€ passing through p. 

The oriented elemental area in equations (3.93) and (3.94) with 
respect to the local basis {0n} is given by the following vec- 
torial product 


v*dA* = dil* x dén (3.95) 
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and since dll* = dl<d% 


*, equation (3.95) can be equivalently ex- 


pressed as follows 


vdA* = dl° ds x den = nigdi*déd' = egg /g*di*déd'’ (3.96) 


where g* = det (sts). 


Moreover, back to the mid-surface we notice it is possible to 


write 
dll x n = vdl (3.97) 
which in the coordinate system {x®, €} becomes 
dl? Oy. XT = Magli d' = €aga/gal™d? (3.98) 


where g = det (gag). 
Equation (3.92) and (3.96) allow us to rewrite equations (3.93) 
and (3.94) as follows 


I n(p, v)dl = f Teas y/g* dl” dtd. (3.99) 
Cc Qe 

m(p,v)dl = | (p* — p) x G€ggr/grdl°déd*? (3.100) 
C Qe 


Next, by virtue of the shifter FÌ, the latter equations become 


f n(p,v)dl = f Teagr/ grdl° dé (dy © d*7) d? (3.101) 

C Qe 

f m(p,v)dl = f En X Teagr/grdldé ( 8 d*7) d? (3.102) 
c Qe 


which, taking into account equations (3.97) and (3.98), become 


| Lewd i i [fo (5, © d*”) vdlde (3.103) 
flo v)dl = ff ea x [Xo (0, @d*) vdldE (3.104) 


and finally 


+e = 
n(p.v) = f g0 (0,049) vag (3.105) 


ni 
m(pv)=nx | ggo(d, ad") vdg (8.106) 


=e 
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where we have put g = \/g*/g 
Both integrands in (3.105) and (3.106) can be further simplified 


just substituting o = od OF Q 3;) and v = vyd® as follows 


n(p,v) = cs go" Jjas) Va (3.107) 


TE 


m(p,v) =ñ x (J go dja Va (3.108) 


and using once again equations (3.85) and (3.86) they assume the 
following form 


+e +e “4g 
n(p,v) = ( go" dé + f go"eus r7) OyVa+ 


TE —€ 


A ( f a go"“dt av, (3.109) 


+e +e = 
m(p,v) =n x ( f go Ed + I goeder) ) Va (3.110) 


where we can finally define two tensors N and M 


N = N” (0,8 Og) + N% (dq ® n) (3.111) 
M = M” (ða & dg) (3.112) 
respectively as 
te te 
N° = i go dé + ni go Ede LY (3.113) 
N% = | gor dé (3.114) 
and 
te te 
Me? = I goh EdE + / goede Lf (3.115) 
such as 
n(p,v) = Nv = Nv, 0g + NE voit (3.116) 


m(p,v) =" x Mv =ñ x M vaðg (3.117) 
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Two fields n and m are called surface stress vector and surface 
couple vector respectively; while the fields N and M are termed 
surface stress tensor and surface couple tensor. 

From the above results it is immediate to notice that the surface 
stress vector n belongs to TgE, consequently it can be split into a 
superficial part and an orthogonal part as follows 


n = nll + nt (3.118) 

where 
nl = NP vaðg (3.119) 
nt = N“y ii (3.120) 


while the surface couple vector m belongs to TQ so that 
m = ml (3.121) 


As the last remak we point out that the coefficient g involved in 
the integration of Cauchy stress tensor along the thickness depends 
only on the geometrical features of the mid-surface Q, in fact it is 
easy to prove the following expression 


g = det (st) = 1+€H+&K (3.122) 


where H and K are the mean curvature and the total curvature of 
the surface Q defined in equations (1.160) and (1.159). 


3.7.3 Body forces and load density 


Suppose the the curve c : IR — Q is closed in such a way as to 
capture a surface portion Q’ C Q bounded by 0Q = c. Assuming c 
to be a directrix, that is a curve through which a line generating a 
given ruled surface always passes, the generatrices directed along n 
define a cylinder Ge(e) C G(e) with thickness 2e and also bounded 
by the surface QU QEU Q_e. 

We assume that the volume forces acting at every point belong- 
ing to the cylinder G,(¢) and the load density acting at every point 
on the upper and lower surfaces Q° and Q_, can be integrated along 
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the thickness to yield a new force system defined on the mid-surface 
Q' as follows 


: Q' > Tg E (3.123) 
:Q-TQ (3.124) 


IO QI 


where g = qf Og + qñ represents the load vector and 5 = © x sP Og 
represents the load—moment vector. 
See [13] for details. 


3.7.4 Eulero’s equations 


The equilibrium equations for the mid surface portion Q’ can 
be written as follows 


1 n(p,»)al+ f qdQ' =0 (3.125) 
aQ! Q' 


J 


f (m(p,v) +7 x n(p,v)) dl +f (Fx g+3)dQ'=0 (3.126) 
dQ! 
which yield 
Nvdl + | qdQ’ = 0 (3.127) 
dQ! ? 
i. (nx Mv+7rx Nv)dl +f (Fx q@+s)dQ’=0 (3.128) 
OQ! , 


Making use of the divergence theorem enounced in equation 
(1.145) on page 29, and due to the arbitrariness of Q’, the above 
equations become 

divN+g=0 (3.129) 
div(ñ x MO, +7 x N°") +F x G+5=0 (3.130) 


Equations (3.129) and (3.130) can be written in components as 
follows 


VIN + LENE +g? =0 (3.131) 
VoN% + Loy N+ =0 (3.132) 
Vi Me — NSE + 58 —0 (3.133) 
eg (L3M® — N°9) =0 (3.134) 
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where equations (3.131) assure the translational equilibrium in the 
tangent plane, while (3.132) represents the translational equilibrium 
along the normal direction. Next, two equations in (3.133) impose 
the rotational equilibrium about the surface axes, respectively. Fi- 
nally, the last equilibrium condition (3.134) gives the symmetry to 
the tensor LSM? = NW. 


PROOF 

Here we want to show all steps we made to pass from the equilib- 
rium equations (3.129) and (3.130) to the corresponding expressions in 
components (3.131) to (3.134). 

Let us start form equation (3.129). We invoke the definition of di- 
vergence for second order contravariant tensors already used in equation 
(1.147), so we have 


(div N)” = NO + Ta NIA+TENG = 
= NYP + NE + Te NYP +72 NE LTE NO LTE NO = 
= NOP + NOS + TR, N° + DE, NE+ 
ATENO + TE NOS + DE N + TE NO 


Now we just need to separate the tangential and normal components 
as follows 


y NE = No8 Lpo NYE 4 pe NY 4 po nae 
(divN)? = NO? + Ta, NIP 478 N° +08, N (3.135) 


(divN)* = NYE + TI NE + TÉ NE + DE NO (3.136) 


By virtue of the the identity (Van)? = LÊ = Dhe equation (3.135) 
becomes 
(divN)® = VEN + LENS (3.137) 
where we have just collected the surface divergence® terms into 


VIN = NYP +Te NW! +r NO (3.138) 


Equation (3.137) proves the in-plane translational equilibrium ex- 
pressed in (3.131). 

Concerning equation (3.136), the translational equilibrium along the 
normal direction is readily proved remembering both Piy = Lay and? 


VaN” = NÆ +T NE + TE NG (3.139) 


8In literature the divergence of the surface tensor N°f is often denoted by 
Noe, 
|a 


°In literature the divergence VÎ,N°£ is often denoted by Ne 
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Hence we obtain 
(divN)§ = Va N° + Loy N (3.140) 


which finally proves equation (3.132) 
In order to prove equations (3.133) and (3.134), first we simplify equa- 
tion (3.130) by taking into account equation (3.129). So it becomes 


divn x Mð, +7 x div (M°"d),) + dive x N°*9,+5=0 (3.141) 


We can split the divergence of the tensor M@" in accordance with the 
results in (3.137) and (3.140), thus we have 


Van x Med, +7 x (VI M°? + LEM) Og+ 
+ x (VE M% + LoyM™) i+ Fax NWO, +5=0 (3.142) 
which after further algebra becomes 
Lð, x M%6,, + Lð, x Mn +n x (VIMP + LEM) g+ 
+0, x N° 6, + bq x NE ñy +7 x dg = 0 (3.143) 


Collecting the normal and tangential terms we obtain the following 
three scalar equations 


Ms (LIM® + N™) =0 (3.144) 


and 
nx (ViMe? — NPS + 8°) dg =0 (3.145) 


which finally proves the rotational equilibrium (3.133) about the surface 
axes. 


Usually a new variable is introduced to make easier possible 
further calculations; in fact we define the pseudo-stress tensor 
the symmetric tensor 


NS? = NOP — LEM" (3.146) 


It is straightforward to notice that N = N only when either a 
membrane stress state holds or for flat shells, namely when Wein- 
garten’s tensor is identically zero. 
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3.7.5 Membrane state of stress 


In this last section we introduce an hypothesis on the state of 
the stress that enables us to derive a closed form solution for several 
shell geometries without invoking the constitutive law. Examples of 
these closed form solutions will be provided in appendix A. 

A shell continuum is subjected to a membrane stress state when 
both the following condition hold 


N% =0 (3.147) 
Me? =0 (3.148) 


Hence, the equilibrium equations become 


VaN +g? =0 (3.149) 
Loy N + 4° =0 (3.150) 
Nag NÊ =0 (3.151) 


where equation (3.149) represents the translational equilibrium along 
the tangent plane; equation (3.150) represents the equilibrium along 
n and finally equation (3.151) states the rotational equilibrium 
about n and establishes the symmetry of N. 


Chapter 4 
Equations of elasticity 


Chapters 2 and 3 of these notes do not specifically concern with the 
elastic media, in fact they can be understood for a generic continuum and 
studied independently. In this section we shall combine the previous results 
in order to to investigate the response of elastic bodies under the action 
of forces. 

A body is called elastic if it has the property of recovering its original 
shape when the forces which produce the deformations are removed. This 
property can be characterized mathematically by certain relationships con- 
necting force and displacement, that are also called constitutive laws. In 
particular we will analyze the linear constitutive law as a generalization 
of the Hooke’s law. 


4.1 The material law 


It was Robert Hooke! who in 1676 gave the first rough law 
of proportionality between forces and displacements for an elastic 
body. In order to understand the key features of elasticity, let us 
consider a thin rod with an initial cross section Ag, which is sub- 
jected to a variable tensile force F. We suppose that the stress is dis- 
tributed uniformly over the area Ao and the initial cross-sectional 
area stays constant. The stress is obtained by dividing the force at 
any stage by the area Ao. So, o = F/Ao. The relationship between 
F and the axial strain e is plotted in figure 4.1 on the next page. 

Figure 4.1 shows that until the point P the relationship o — € 


‘Robert Hooke (July 18, 1635 Freshwater (Isle of Wight) - March 3, 1703 
London) was an English scientist. 


b 
Ft 


Source: http://turnbull.mcs.st-and.ac.uk/history/Biographies/Hooke.html. 
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B 


9 € 
Figure 4.1: Hooke’s law. 


is nearly a straight line with the following equation 
o= Ee (4.1) 


where the constant of proportionality E is known as modulus of 
elasticity or Young’s modulus. 

The greatest stress that can be applied to the rod without pro- 
ducing a permanent deformation is called elastic limit of the mate- 
rial. When the force F is increased beyond this limit the material 
goes in the elastic-plastic field. Namely, firstly the material reaches 
the yield-point Y at which the rod suddenly stretches, then the ma- 
terial reaches the ultimate stress at U where it offers the maximum 
stress. If the elongation increases again both the cross sectional area 
Ao and the stress decrease until the rod breaks at B. 

From now on we shall study only the elastic range. 


4.1.1 Generalized Hooke’s law 


Here we want to extend the results of Hooke’s law to a multidi- 
mensional state of stress and strain. So, in accordance with equation 
(4.1), let us write a linear relation 


Tij = CijhkEhk i j,h,k = 1,2,3 (4.2) 


The coefficients Cijnk are independent from the position of the 
reference point in the continuous medium, in other words we require 
the homogeneity of the body, that means uniformity in structure 
and composition. It can also be shown that the elastic constants 
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Cijnk are 81 components of a fourth order tensor which is termed 
elasticity tensor. 

Since the stress tensor oj; is symmetric, an interchange of the 
first two indices in (4.2) does not alter its meaning. In addition to 
that, the symmetry of the strain tensor ensures also the symmetry 
of the last two indices, so that 


Cijhk = Cjihk (4.3) 
Cijhk = Cijkh 


That means that the 34 components of C reduce to 36 indepen- 
dent constants. Let us show the expansion of a generic component 
of the stress tensor, that is 


011 =C1111€11 + C1112€12 + C1113€13+ 
C1121€21 + C1122€22 + C1123€23+ (4.5) 
C1131€31 + C1132€32 + C1133€33 


Equations (4.3) and (4.4) allow (4.5) to be rewritten as follows 


o1 =C1111£11 + l1122€22 + C1133€33+ 
2C1112€12 + 201113€13 + 2C1123€23 


Thus, the whole elastic matrix can be written as 


011 Cini Cree C1133 2C i112 2C1123 201131 €11 
022 C2222 C2233 2C2212 2C2223 2C2231 E22 
033 |_ C3333 2C3312 2C3323 2C3331 €33 
o | 2Ci212 2C1223 2C 231 €12 
023 sym. 2C2323  2C2331 €23 
031 203131 €31 


which, making use of the symmetry relationships expressed in (4.3) 
and (4.4), simplifies as follows 


011 Ci €12 C13 C14 C15 C16 Ell 
022 C22 C23 C24 C25 C26 E22 
033 _ C33 C34 C35 C36 €33 
012 C44 C45 C46 E12 
023 sym. C55 C56 E23 


031 C66 E31 
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Later on, see equation (6.10), we will also introduce another 
symmetry condition that has been assumed in the above. Namely, 
the condition 


Cijnk = Chkij (4.6) 


that further reduces the independent elastic constant from 36 to 
21. So, the latter material equation represents the constitutive law 
for an anisotropic elastic material. However, most of the engineer- 
ing materials have some symmetry properties which allow further 
reductions of the elastic constants. 

The highest degree of symmetry leads to the so called isotropic 
material. We define an isotropic material an elastic continuum which 
has the same response in any direction, so that the elastic tensor is 
not influenced by any rotation of the references axes. 

Let the elastic tensor be represented by Cijnk with respect to 
the cartesian coordinate {x*} whose basis is B = {&;}. With respect 
to a rotated system {x} with basis B’ = {&{} the elasticity tensor 
is Cisne: By the definition of isotropic material, we expect that the 
elasticity tensor does not change. In order to show this, let us recall 
the transformation relations (1.36) on chapter 1. Here we are dealing 
with a Cartesian coordinate system, hence it does not matter if the 
indices are all subscripts. So, we have 


1 USA 
Cijnk = Uila jmCimnoGohAnk 
keel t 1 
si ditAjrmAhodknClmno (4.7) 
but to ensure the immunity against the rotation of the reference 


system, we impose 


Ci;nk — Cirino (4.8) 


that is only satisfied if the elasticity tensor assumes the following 
form 
Cimno = Adimbno + 101nOmo F Kôloðmn (4.9) 


where A, u, k are elastic constants?. 


?This can be proved by replacing equation (4.9) into (4.7), as follows 


if 1 / / I 
Ciink = ditljm@hn0%kn A imino + LOindmo + KdtoÎmn) = 
= 4! 1 II TIGER DIR: DIRO | ERA 5 
=, dimAjmAhodko on Ha;nAj04hnAko ne Kaj0AdjnAhnko o 


bij Onk + udinòjk + kd;kÒ;h 


that is exactly the expression (4.9). Note that we have used the identity aj,,a,, = 
Ôpq provided by equations (1.21) and (1.24) on page 7. 
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In equations (4.3) and (4.4) we have already noticed the sym- 
metry of C in relation to the two front and two back indices, let us 
show now that one more reduction is possible 


Cijhk = AdijÒnk + udinòjk + KÒikÒjh (4.10) 
Cijkh = AdijÒkh + udikòjh + KOihOjk (4.11) 


where, subtracting term by term and considering the symmetry of 
the unit tensor 6;;, equations (4.10) and (4.11) lead to the only 
possible condition 


H (Oindjk — dikdjn) + K (OikOjn — Oindjr) = 0 > 
H (dindjk — OikOjh) — K (Oindjk — dikdjn) = 0 > 
(u — K) (indik — Sikdjn) = 0 (4.12) 
which is only true if (u — x) = 0. So, the relationship between & 


and u further reduces the number of elastic constants to 2. Namely, 
we have 


Cijnk = AdijOnk + u (Oindjk + dikdjn) (4.13) 


The Hooke’s law becomes 


Cij = Cijnkenk = Adi; OnkEnk + H (dindji + dikò;n) Ehk = 


AĝijEhh + 2HEij (4.14) 


where we have used hkEhk = Ehh = tTEhk- 

Equation (4.14) is the generalized form of Hooke’s law, valid 
for homogeneous, isotropic, elastic bodies. A and p are called Lamé 
constants? , 


*Gabriel Lamé (July 22, 1795 Tours - May 1, 1870 Paris) was a French 
mathematician and engineer. 


od 


Source: http://turnbull.mcs.st-and.ac.uk/history/Mathematicians/Lame.html. 
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The trace of the stress tensor is readily computed by contracting 
the indices, so that 


Cii = 3AEnh + 2uE > (4.15) 

Ti = (2u + 3A) Enh > (4.16) 
Tii 

hh (2u + 3A) (4.17) 


where we can put troj; = 0; = X and tre;; = ey = O. 
The above expression (4.17) is useful if we solve (4.14) for ej. 
In fact, we have 


1 
Eij = Qu = pe (4.18) 


and in observance of (4.17) we obtain 


1 A 


= HA Ôi DI 4.19 
Bu 2u (8A + 2p) ” ey 


Now, let us consider an axial state of stress. The stress tensor is 
011 0 0 
Tij = 0 0 0 
0 00 


form (4.19) we have 


€11 = 1 1 À = 
T Gram) 17 


A p 
— i ce elet 4.20 
16142) 0!! (4.20) 

À 
= = SS i 4.21 
E22 E33 24 (BA ce 21) O11 ( ) 
(4.22) 


Let us define Poisson’s ratio v as follows 


pace i (4.23) 
€22 E33 2(4+A) 
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À u=G E 7 
di i I a = OD) 
A, V n Sa I) g 
pb) = - š Eu 
Hik re 5 2u(1+v) - 
Ev | TE = ; - 


Table 4.1: Relationships between the main elastic constants. 


According to Hooke’s law in the original form, see equation (4.1), 
we can see that 


1 A= 4 _ w(3A+ 2p) 


= => F= 
E w(3rA+ 2p) A= u 


(4.24) 


So, we have proved that Lamé constants can be replaced by 
E and v which lead to writing the alternative expressions of the 
constitutive law 


iù 1+yv) 04; — Udi; 4.25 
ISE j j 
E v 
Ox 2 (es + ue) (4.26) 


Table 4.1 shows the relationships between elastic constants. 


4.2 The linear elastic problem 


In this section we are going to sum up equations and unknown 
quantities which define the classical linear elastic problem. Then 
we will estimate the distribution of stresses and strain as well as 
displacements at all points of the body when certain boundary con- 
ditions are given. Let us balance the unknowns and the equations, 
we have fifteen unknowns (6 stress components + 6 strain compo- 
nents + 3 displacement components) for all points in the continu- 
ous and just fifteen equations (6 equilibrium + 6 compatibility + 3 
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boundary conditions). So, for a given linear elastic body V we have 


4.27 
4.28 


4.29 


) 
) 
) 
4.30) 


( 
( 
( 
( 


In order to solve the linear elastic problem we start from the 
known quantities (4.27) to (4.30), and through the following avail- 
able equations 


- compatibility equations 


Eij = ; (ug +u;i) onV (4.31) 

- equilibrium equations 
Cijj+b =0 onV (4.32) 

- constitutive laws 

Tij = i (e + ide) on V (4.33) 

- boundary conditions 
Tijnj = f on So (4.34) 
ui = ûi on Su (4.35) 


we will formulate two boundary-value problems. 


4.2.1 Boundary value problem in terms of stresses 


This first boundary value problem can be stated as follows: 


Determine the distribution of stresses and displace- 
ments in the interior of an elastic body in equilibrium 
when the body forces are prescribed and the distribution 


of the forces acting on the surface of the body is known". 


4Sokolnikoff [1]. 
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Following the above formulation, the procedure for solving the 
problem would suggest writing the available equations entirely in 
terms of stress. To this aim let us start from equation (2.75) 


Eijhk + Ehk,ij — Eih,jk — Ejk,;ih = 0 (4.36) 


and consider the constitutive law (4.25), so that 


l4v 
E. (Cijhk Ong Chak Oe = 
V 
=F Onde + ÔhkOnn,ij — Om Omak — ÎjkOnn,ih) (4.37) 


Equation (4.37) represents a set 34 = 81 equations since all the 
four indices è, j, h, k run from 1 to 3. Not all of these equations are 
independent, indeed the system (4.37) contains only 6 independent 
equations. A first reduction of equations is due to the contraction 
h=k that yields 


Oij,kk + Okk,ij — Tik,jk — Tjk,ik = 
V 
ey (dijOnn,kk + dkkOnn,ij — dikOnn.jk — djhOnnsik) (4.38) 
that, by denoting Y = tro;j = o; and oj; ~~ = Woy, becomes 
V 


2 2 
E (dij V E + Vj) (4.39) 


2 
V‘ oij +Dij — Tik,jk — Tjkik = 
By virtue of the equilibrium equations (4.32), the above expres- 
sion can be rewritten as follows 


1 V 
2 2 


which is a set of 6 independent equations. 

Next, in order to express V?¥ as a function of the body force 
b, we put h = i and k = j in equation (4.37), so that, after a bit of 
algebra, we have 


v 
Cijig = VE — CEF n Vi 
1-V9 


= >: 4.41 
1+ oe ( ) 
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and finally, by invoking the derivative of the equilibrium equation 
that provides the relationships bj; = 0ij;ij;, we get 
l+»v 


To pei (4.42) 


2 
Y= 
V = 


Now, going back to equation (4.40) and making use of the latter 
result, it is not a difficult task to obtain the following expression 


v pe 
Tag A = @ + bjii + 1 pbgdiv i) (4.43) 


V7 oj + 
Equations (4.43) were derived by Michell® in 1900 and by Bel- 
trami in the 1892 for the special case when the body forces are 
absent. Nevertheless, it is common to refer to equation (4.43) as 
Beltrami-Michell equations. 
In case of missing or constant volume forces equation (4.43) 
assumes the straightforward form 


1 
Voi; + fay? = 0 (4.44) 


4.2.2 Boundary value problem in terms of displace- 
ments 


The second boundary value problem can be stated as follows: 


Determine the distribution of stresses and displace- 
ments in the interior of an elastic body in equilibrium 


>John Henry Michell (October 26, 1863 - February 3, 1940) was an Australian 
mathematician. 


Source:http://en.wikipedia.org/wiki. 
ĉEugenio Beltrami (November 16, 1835 Cremona - February 18, 1900 Rome) 
was an Italian mathematician. 


$ 


Source: http://www-groups.dcs.st-and.ac.uk/ history/Biographies/Beltrami.html. 
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when the body forces are prescribed and the displace- 


ments of the points on the surface are prescribed func- 


tions’. 


By replacing the constitutive law in the form of (4.14) into equi- 
librium equation, we obtain 


(AdijEkk) ; + 24€ ij,j +b; =0 (4.45) 


that is 
VEkk,i + 2HEijj + di = 0 (4.46) 


and in accordance with the compatibility equations we have 


Uk ki + 4 (Ui jj + Ujij) +b; =0 (4.47) 
Nup,ki + UV? Us + Hukk + di = 0 (4.48) 
(A+ H) Ukki + uV?u; +6; =0 (4.49) 


that in the vectorial form reads 
(A+ u) grad diva + uV?a+b=0 (4.50) 


Equation (4.49) (or equivalently equation (4.50)) is called Lamé- 
Navier equation and together with the boundary conditions ex- 
pressed by equation (4.35) define the boundary problem inn terms 
of displacements. 

Once the first boundary value problem has been solved, i.e. when 
the displacements are known, the state of strain and hence the 
state of stress can be found though equations (4.31) and (4.33), 
respectively. 

Further attention should be focused on the case when body 
forces do not occur or they are constant. First, consider the di- 
vergence of equation (4.49) 


(A+ H) gii + UV? wii + dii = 0 (4.51) 
that yields 


AV Ub + 24V? Ukk + bii = (A + 2p) V2ukk + bi = 0 (4.52) 


"Sokolnikoff [1]. 
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which, under the hypothesis of b; = const., so that bi; = 0, gives 
V7 un, = VO =0 (4.53) 


where we have set © = tre;j = Ei. 
Moreover, recalling (4.17) it is also proved that 


Von: = 0 (4.54) 


We can finally say that if the volume forces are constant, the 
boundary linear elastic problem in terms of displacements turns 
into a general boundary values problem of a biharmonic differential 
equation. 


4.3 Constitutive equation for shell continuums 


The Kirchhoff-Love hypothesis and the inextensibility of mate- 
rial fibers along n allows one to consider the shear stress components 
NS® unrelated to strains, so that the constitutive problem can be 
solved through the plane stress model. Thus, components N§° are 
found only by means of the equilibrium equations. The analytical 
derivation of the constitutive equations is beyond the scope of this 
book, so we will just present the final equations that will be used 
in the appendix A in order to solve some case studies. However, 
readers can find thorough discussions in [3] and [16]. 

Suppose a membrane state of stress, the constitutive equations 
are the following 


N08 = DH May, (4.55) 
M8 = BHM Ey, (4.56) 
where 
l-v 2v 
He eG a go) n 


The fourth-order tensor H®# has the following symmetries 


HOO — ppodu — pobuà — pPruas (4.58) 
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Finally, coefficients D and B are the in-plane and the bending 
stiffness, respectively, defined as 


DE < (4.59) 


A (4.60) 


Chapter 5 
Principle of Virtual Work 


This chapter is entirely devoted to the Principle of the Virtual Work. 
In particular the relations between equilibrium, compatibility conditions, 
and virtual work will be highlighted. 


5.1 Virtual work 


Virtual work can be defined as the work done on a deformable 
continuum by all the forces acting on it when the body is subjected 
to a small hypothetical displacement field - unrelated to the forces - 
which is consistent with the constraints present. The latter is named 
virtual displacement and is denoted by an asterisk. 


Suo | H 


vA S 


I y N 


I 
Figure 5.1: Forces and constraints acting on the continuous. 
As figure 5.1 shows, let us suppose to split the surface S into two 


separated boundary surfaces, in such a way that surface forces are 


Claudio Borri, Michele Betti, Enzo Marino (a cura di), Lectures on Solid Mechanics, ISBN 978-88-8453- 
854-3 (online), ISBN 978-88-8453-853-6 (print), © 2008 Firenze University Press 
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prescribed on S, and a boundary displacement field © is prescribed 
over the remaining boundary surface denoted by Su. Namely, the 
entire surface results as the sum S = Su U So, where Su = S}, U S! 
and So = S} U SX. 

Consider now a virtual displacement field už which yields the 
deformation 


1 
Ey = 3 (CHE + ui ;) (5.1) 


In order to calculate the virtual work done by the volume and 
surface forces we write 


w= | futds, + | butdV (5.2) 
Se y 


and we recall also the equilibrium equations discussed in chapter 3 


Ciji + bj = 0, VpeV (5.3) 
cijnj = fi, Vp € Sc (5.4) 
Cij = Oji Vpeyv (5.5) 


By replacing both equations (5.3) and (5.4) into equation (5.2) 
we have 


w=f aisnjuf dS, — | onipuido (5.6) 
So v 


Next, consider the following identity 


f oi dv = I ed: if crutzdV (5.7) 
y i y y 


that through some simple algebra yields the following equations 


f Comu”) a = | omnuzave 
V V 


1 1 
+5 I Thi (uj), + Up i) dV + 5 | Chi (uz n — Ups) dV = 
v v 
eH mmm 
=0 


= f omnia 4 | onsi + 
V V 


i, Fug = I (oniut) n dV — f mew (5.8) 
y y y 


LECTURES ON SOLID MECHANICS 105 


where we have split the virtual displacement gradient u% xa Ei + 


wij, See equation (2.81), and used the fact that the product of a sym- 
metric tensor 0;j by a skew-symmetric tensor wij = 3 (ut, = už) 
always vanishes. Moreover, the divergence theorem allows us to 
write 


f (Chi) h dV = Oni; NAAS (5.9) 
V SoUSu 


so finally equation (5.8), provided that u* = 0 on Su, becomes 


w* a. cijnjuj dS — aniuimnas+ | oniehidY (5.10) 
So So Vv 


=0 


As a result we have proved the following expression, also known 
as principle of virtual work, holds 


fuzdse + | bai = | ez; (5.11) 
So V V 


In equation (5.11) we shall define the left-hand side group of 
terms as external work, and the right-hand side one as internal 
work, the reason why the following alternative names are often used 


w*= L= f fiutdS, + f but dV 
o y 

w* = L; = f CijEi;V 
y 


To obtain the above results we started from the equilibrium and 
compatibility conditions and it is interesting to notice that we have 
never used any constitutive laws. So the PVW can be applied to all 
continuous material with the only limitation of small displacements. 

Let us take a look to the physical meaning of the internal work. 
Consider an infinitesimal volume element dV, shown in figure 5.2. 

In a two dimensional case the work done by forces acting on dV 
for each deformation Ei; can be seen in sketches 5.3, where the axial 
virtual dilatation along the x direction and the angular dilatation 
in (x1, £2)-plane are depicted. 

Thus, the infinitesimal virtual work done by 022dr1dx3 is 

dL; 


* * 
iza m O22dL1dL3E32dL2 = 0726994) 
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da! 


Figure 5.2: Elemental volume element. 


o La o Tz 
Endr? 
~ 
21 
T Tı X 
i da da! 
da? da? 
(a) Axial dilatation. (b) Angular dilatation. 


Figure 5.3: Virtual deformation. 


In the same way from figure 5.3(b) the work done by the force 
09dx1dx3 is 


dL; 


* * * 
o 091dxdx393dx2 = 2091€9dxdx3dx3 = 2091 €23 dV 


Computing the infinitesimal work dL* done for each deformation 
and integrating on the entire volume V we obtain 


L= f cijEijdV (5.12) 
V 
The principle of the virtual work for rigid bodies can be readily 
derived from the general expression (5.11) where, since ej; = 0, 
yields 


L*=0 (5.13) 


e 
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5.1.1 A simple example 


As an application of the virtual work principle let us consider a 
simple system of bars! 

All the forces are summed in a resultant force F applied to the 
end point of two bars, as shown in figure 5.4. The external force 
provides only an axial state of stress in the bars. 


F 150kN 

E | 206 x 10° N/mm? 
A 15cm? 

l 2m 

li 4.47m 

lo 4m 

a 63.435° 


By means of the PVW we want to find the magnitude of the 
real displacement at point A along x; and xa-directions under the 
effect of F. To this end we consider first a unit explorer load in 
x-direction (to compute uj, then a unit explorer load in the z2 
direction (to compute us). 

First of all we solve the equilibrium problem, so that 


real case a case b case 


N, 2F/sina 2/sina 0 
Na —2F 2 1 


Consider case a in which we shall compute uj. PVW reads as 
follows 


Lé=luj= | Mejdl+ | Noe5dl = L; 
ly lo 
where 
sr D 
E = = 
1” EA EAsina 
se 
2 EA EA 


1 Although no notions on mechanics of beams or frame structures have been 
introduced so far, the intuitive meaning the reader can assign to some quantities, 
like the unit axial deformation £, could be enough to get the main idea of this 
example. Moreover, after reading chapter 8, the reader will be able to have a 
more comprehensive view of this application. 
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Figure 5.4: Example. 


So we have 


For case b we write 


Peirs i Medar 
l2 


where 
«_ Nx 1 
6 = — = - = 
27. EA EA 
So we have 
Si 2F 
lus = Fa” 


Finally, with respect to the cartesian positive direction, using 
the above expressions we easily obtain 


uj = 18.62mm 


Uy = —3.88mm 
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5.2 PVW, Compatibility conditions, Equilibrium 


In the previous section we showed that the equilibrium and the 
compatibility conditions lead to PVW. Now in this section we shall 
prove that two of three conditions are enough to obtain the third 
one. 


PVW + Compatibility > Equilibrium. Let us start from the 
PVW and the compatibility equations 


f totas + | ugar =f cea» 
S V y 


1 
Ey = z (wig + Uju) 


By splitting the displacement gradient the PVW can be rewrit- 
ten as follows 


fsutas+ | ugar =f oyiw- f oywyav 
S v V v 


Now, making use once again of the identity (5.7), the latter 
equation becomes 


J fiutdS + ni biutdV = 
S y 


fc, dv- | gia — | oyw > 
fot = f Coun- fias- f oway 


where the divergence theorem has been used. 


Due to the arbitrariness of the displacement field, if com- 
patible, the latter equation is only satisfied if each argument 
vanishes, therefore 


Gij j + bi = 0, VpevV 
oijnj— fi=0,  Vpe So 
CijWij = 0> 


Cij = Oji VpevV 
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PVW + Equilibrium > Compatibility. To prove this statement 
consider the expression of the PVW (5.11) 


fugds + | tnt = | cet; 
So V V 


and the equilibrium equations 


Ciji + bj = 0, Vpev 
Cijnj = fi, Vp ES, 


that replaced into the equation of the PVW lead to 


1 oimsuîdS, — | ongutay = | onea 
So V y 


and by using again the divergence theorem the above equation 
becomes 


[ (osi); - | ops = f oyeyav 
v ? v ¥ 


and by expanding the first integral on the left-hand side we 
obtain 


forio = f oscar > 


1 1 
[ews (uj; + ui ;) dV + [ess (uj — už ;) dV = 


=0ijWij=0 
| GijEi; dV 
v 


that proves the compatibility condition of the virtual displace- 
ment field i 
Ey = 5 (Wij tia) 


Chapter 6 
Energy principles and variational methods 


This chapter deals with some of the most important results concern- 
ing energy principles in elasticity. Firstly, we will start introducing the 
strain energy and how it is related to the work done by the external forces 
(Clapeyron’s theorem), secondly we will introduce two important theorems 
which make use of the strain energy: uniqueness of the solution for the 
elastic boundary-value problem and the theorem of reciprocity. 

Finally, the equilibrium condition will be interpreted as the stationary 
condition of the potential energy and accordingly some energetic theorems 


will be enounced. 


6.1 The strain-energy function and Hooke’s law 


Suppose the body V lies in a natural state at time t = 0. Under 
the effect of surface and body forces, the continuum has at time 
t the strained configuration. With respect to the usual cartesian 
system each point of V is found by x; + uj (x;,t). Where {x;} de- 
notes the coordinates in the unstrained configuration at t = 0. The 
displacement fields may be derived as follows 


Ou; 
ot 


dt = ù;dt 


During the deformation the work done by all the forces acting 
on the body is denoted by W. The rate of W is given by 


dw 
ETA f fiùidS + l dii; dV (6.1) 
By making use of equilibrium boundary equation (3.3.3) in chap- 


Claudio Borri, Michele Betti, Enzo Marino (a cura di), Lectures on Solid Mechanics, ISBN 978-88-8453- 
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ter 3 and the divergence theorem, equation (6.1) becomes 
dW ; : 
E = [ruwa f Costa), dV = 
= f biù;dV + ‘| Cij,jttdV + si oj; ts jdV = 
V V V 
= f biùidY + f Oij jùidY + | Gij (Éi; + Wiz) dV = 
V V y 
= f bid + [oj jd + | Gij (Éi; + wij) dV 
y V y 


Due to the equilibrium condition, and remembering that c;jwij = 
0, the latter equation may be written as follows 


Ww = Foi = |o iw (6.2) 
Now let us suppose that there exists a function ¢ = @ (€;;) such 
that 
Da = (6.3) 
so, (6.2) becomes 
w = 1 = “i dV = a fos) dV (6.4) 


Let us define strain energy the following integral 


= | è (eij) dV (6.5) 


where @ is said volume density of strain energy or elastic potential. 

Hence, since we are considering the instant t when the body lies 
in an equilibrium configuration, so that the kinetic energy vanishes, 
equation (6.4) states that the work W done by the external forces 
in altering the configuration of the natural state to the equilibrium 
state at the instant t is equal to the strain energy ®. Therefore, the 
latter can be considered as the energy stored in the deformable body 
when it is brought from an initial natural state to the equilibrium 
state. 


LECTURES ON SOLID MECHANICS 113 


We assume now that the strain energy density function ¢ can 
be expanded in a Mc Laurin series 


oe =0+ (3) +3( ca ) cent (6.6) 


de; OEijOEhk 


where we discard all terms of order 3 or higher. The constant @ (0) 
is the energy density associated with the initial stress state while 


( oq ) 0 
de; 
is the initial stress state. 


Now recalling the Hooke’s generalized law (4.2), and taking into 


account equation (6.3), we obtain 


Osi; 


TR = Cink (6.7) 
and 3 22 
Si 0 (6.8) 
OEnk OE; KOEI; 
so from equations (6.7) and (6.8) we have that 
ao 
(A 6.9 
JRR depk0eij ( ) 


and due to Schwartz’s theorem the symmetry of the elasticity tensor 
has also been proved 
Cijhk = Chkij (6.10) 
This important result can be substituted into equation (6.6), 
where, by assuming that both the energy density at initial state of 
stress and prestresses vanish, we obtain 


1 
$ (eij) = z CishkEijEhk (6.11) 
hence i 
= | è (€;) dV = J CijhkEijEhk (6.12) 
V V 


In the case of a simple axial state of stress the density of strain 
energy is 


1 1 
d= gOunenen = e 


therefore the dashed area in figure 6.1 represents the density of 
strain energy for an axial state of stress. 
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€14 


Ø 


l I 
O11 


Figure 6.1: Density of strain energy in the case of axial state of 
stress. 


Theorem 1 (CLAPEYRON’S THEOREM) If a body is in equilibrium 
under a given system of body forces b; and surface forces f;, then 
the strain energy ® is equal to one-half the work done by the external 
forces (of the equilibrium state) acting through the displacements ui 
from the initial state to the equilibrium state. 


To prove Clapeyron’s theorem let us recall the PVW expres- 
sion. See equation (5.11) in chapter 5. To do this it is necessary 
to require the equilibrium state of the body and the consistency of 
displacement and strain fields, thus 


fiuidSo +f bju;dV = f oijeijdV (6.13) 
So v V 
where replacing the expression of the strain energy (6.12) we have 


28 = fiuidSo +f biuidY (6.14) 
So V 

On the right-hand side of equation (6.14) we recognize what we 

have defined as external forces work. Therefore we have proved that 


1 
b= SL (6.15) 


As done in equation (6.3) we suppose now the existence of the 
conjugate strain energy density 


d* = $* (cij) (6.16) 
such as 


ap _ 


00%; 


Eij (6.17) 
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so we define 
ğ* = | o (sij) dV (6.18) 
y 


as the conjugate strain energy. 
Through equations similar to the strain energy case, we can 
write 


* La 
p* (cij) = 5 CishkTijOhk (6.19) 


where Ci;hk = Oi and €;j = CishkThk- We can also prove that 


š 1 
è (Eij) = & (cij) = getij (6.20) 
in fact, we have 
* RS 1 
$“ (oij) = 5 CishkTijOhk = gfig0ij = è (Eiz) (6.21) 


Often it is useful to know the whole strain energy of a deformable 
body without knowing the internal state of stress. Figure 6.2 shows 
the deformation of a beam under a concentrated external load. 


u 


a] 


Em» 


F Pe Le 


Figure 6.2: The simplest application of Clapeyron’s theorem. 
We can compute the strain energy easily as 


1 


1We suppose the self weight vanishes. 
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6.1.1 Superposition principle 


We set the elastic equilibrium boundary-value problem as fol- 


lows 
divo+b=0  VpeV (6.22) 
e=symVu Vpev (6.23) 

E v 

= t .24 
o (+e) VpeVvV (6.24) 

and 
on= f Vp E So .25) 
u=u YpESu (6.26) 


Superposition principle is a general tool that can be applied to 
many physical linear systems. It states that if a number of inde- 
pendent influences act on the system, then the resultant influence is 
the sum of the individual influences acting separately. 

Namely, let us suppose to have two systems of body and surface 
forces 


(6, FO} (6.27) 
pa, JO} (6.28) 


Every force system is related to the following strain and stress 
state, respectively 


fa e, 300} (6.29) 
fa, 2), gd} (6.30) 


Formally the the principle of influence superposition allows us 
to state that for all A1, A2 € JR, given the following forces system 


{Ayo + dob?) , AF) + of} (6.31) 
then the following set of displacement, strain and stress 
{GY + Agu), Axe) + Age, Apo + Ago} (6.32) 


is the solution of the equilibrium boundary-value problem. 
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To prove that it is enough to substitute the above fields in the 
elastic problem equations. 

However, there are some special cases for which the superpo- 
sition principle does not hold. Indeed, let us consider the strain 


energy ® and put £;;j = a + ee We obtain 
p (či) = a CijhkEijEhkdV = = (6.33) 
1 
= ni Cisne ( ssi +e) (e +e) dV = (6.34) 


= 5 f Comet? OVF 5 | Cime Pea 


1 
+3 F Cie e dV + 5 T Cipe Dell) dV (6.35) 
Finally, due to the symmetry of the elasticity tensor we have 
© » - © va ee ef) = (6.36) 
1 
= © (cP) +e ey Ce Get Day (6.37) 


The last term in equation (6.37) represents the coupling con- 
tribution to the strain energy which disproves the superposition 
principle for the strain energy. 


6.1.2 Uniqueness of the solution 


The solution of the boundary—value problems formulated in sec- 
tions 4.2.1 and 4.2.2 is unique. To show that let us assume, by ab- 
surd, that it is possible to obtain two solutions for the boundary— 
value problem expressed by equations (6.22) to (6.26) 


fa), 0), o} (6.38) 
fa), A o) (6.39) 
Equation (6.22) allows us to put 


o = 03 +b =0 (6.40) 


while equation (6.25) allows to write 


ony = = oe n; = fi (6.41) 


118 ENERGY PRINCIPLES AND VARIATIONAL METHODS 


Now, by virtue of the superposition principle, it is clear that the 
following function 


represent a special solution fulfilling the equilibrium equation pro- 
vided that b; = 0 and fi = 0. In fact we have 


1 2 
Ong = (0) — of) , =0 (6.42) 


Cijjnj = (0) — 05) nj =0 (6.43) 
Thus, for this special solution Clapeyron’s theorem (6.14) writes 
as follows 


26 =0 (6.44) 


that is 
f È (eij) dV = f CijhkEijEhkdV =0 (6.45) 
V v 


But since @ is a positive definite quadratic form, the above in- 
tegral can vanish only when ¢;; = 0 and so the following identity 
has been proved 


e) = (6.46) 


Therefore, if the components of the strain tensor for two solu- 
tions must be identical, then, by means of the constitutive law, it 
follows that the components of the stress tensor must be identical 
as well 


Finally we want to remark that the equality e) = ¢(2) does 
not exclude rigid body motion vo. In fact ul = Ñ? + uo satisfies the 
uniqueness of deformation since it does not produce any deforma- 
tion, i.e eo = 0. But we shall suppose the constraints along S,, are 
able to inhibit all rigid displacements, so that 
a) = a (6.47) 


i 
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6.1.3 Theorem of reciprocity 


Now we introduce a general reciprocal expression relating the 
equilibrium states of a body under different loads. To do this let 
us consider two equilibrium states of an elastic body one of which 
subjected to displacement field & due to the body and surface forces 
b and f, respectively; the other equilibrium state is characterized 
by the displacement field & due to the body and surface forces b' 
and f’, respectively. The work that would be done by the forces b; 
and f; if they acted through the displacements u; can be written as 


follows 
f tas | ootav = | onea (6.48) 
S y V 


and in the same way, the work that would be done by the forces b; 
and f’ if they acted through the displacements u; can be written as 


follows 
i fiujdS + | buid Y = i, 0j;EijdV (6.49) 
S v V 
Through the symmetry of the tensor of elasticity we notice that 


Gij = CijhkEhk 


Tij = CijhkEhk 
hence, 
fore = (6.50) 
y 
= | CismeeisV= | Cruiseisehed = (6.51) 
V V 
= f ced. (6.52) 
y 


Equations (6.50) and (6.52) prove that (6.48) and (6.49) are 
identical, so 


[todas [outer = | fas + | tanav (6.53) 
S V S V 


Equation (6.53) can be enunciated in the following theorem 


Theorem 2 (BETTIS THEOREM) If an elastic body is subjected to 
two systems of body and surface forces, then the work that would be 
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done by the first system b; and f; acting through the displacements 
ul, due to the second system of forces is equal to the work that would 
be done by the second system of forces bi and f! acting through the 
displacements u; due to the first system of forces. 


Theorem of Betti is also termed theorem of reciprocity. 


6.2 Variational methods 


In this section we present an alternative approach aimed at find- 
ing the state of stress in a continuous elastic body. This method 
basis itself on some minimum principles that characterize the equi- 
librium state of bodies. Namely we shall see that it is possible to 
construct some integrals relating the work done by the forces act- 
ing throughout the deformation and to show that these integrals 
have their minimum values when the distribution of stress in the 
body corresponds to the equilibrium states. So searching for equi- 
librium state is reduced to certain standard problems of calculus of 
variations. 

This method is strongly used in computational mechanics to 
solve the equilibrium problem in the finite elements method. 


6.2.1 Potential energy 


Let us start by introducing the functional U called potential 
energy of deformation. We shall show that this potential attains an 
absolute minimum value when the displacements of the body V are 
those of the equilibrium configuration. As usual, see also figures 5.1 
on page 103 and figure 6.3, let b; be the body forces and f; the 
surface forces prescribed on Sc. 

We suppose that over the remaining part of S, i.e. S,, the dis- 
placements û; are known. We denote the displacements which sat- 
isfy the equilibrium configuration as u; and consider an arbitrary 
small? displacements du; only if consistent with respect to the com- 
patibility conditions imposed over Su. To the field du is also re- 
quested to belong to class C? and to be zero over Su. We shall term 
du as virtual displacement field?. 


?Compatible with the hypothesis of linear elasticity, de;; = 3 (Sui, j + du;,i). 
3Note that in chapter 5 the virtual displacement 5u was denoted by &°. 
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S| |F 
Ci 
ù 


Tı 


Figure 6.3: Forces and displacements acting on the body V lying in 
the equilibrium state. 


The elastic problem is given as follows 


divo +b=0  VpeV (6.54) 
€ = sym Vù VpEV (6.55) 
n= L (e + i Zire) VpevV (6.56) 
and 
uw=u WES, (6.57) 


We define the potential energy of deformation U as follows 


uU (Eijs ui) = (eij) — Y (ui) (6.58) 
where 
2 (e) = / o(eyyay (6.59) 
V 
LU (ui) = f iuw + [ fiuidS (6.60) 


The potential energy U is the sum of the strain energy ® and 
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the conservative loads potential Y. The first variation of U is 
ôU =U (u + ôu, £ + de) -U(U,E) (6.61) 
= © (e + de) — ŲĮ (ui + ĝui) — P (€) +Y (ui) = (6.62) 


= ğ (e + ôe) — | bi (u; + du;) dV— i fi (ui + du;) dS+ 
V Ss 
= ® (€) + / bjugdV + | fiuidS (6.63) 
V Ss 
In equation (6.37) we have just seen that 
p (e + de) = (€) + (de) +f CijhnkEnkÒEij (6.64) 
V 
so neglecting ® (de;;) as a second order infinitesimal, we obtain 


U= | ayseyav—f réuar— | féuas (6.65) 
V V S 


We recognize in the right-hand side terms of the above expres- 
sion the Principle of Virtual Work (5.11). So it is easy to notice 
that the stationary point for potential energy U corresponds to the 
equilibrium condition. We know, in fact, that the PVW and the 
compatibility conditions of displacements du; result in the same 
thing and that is equilibrium state. 


Theorem 3 (STATIONARY VALUE OF POTENTIAL ENERGY) The total 
potential energy U of an elastic body has a stationary value in the 
class of the geometrically permissible displacements for the true dis- 
placements which correspond to the state of equilibrium. 


It is also possible to prove the following stronger theorem. 


Theorem 4 (MINIMUM POTENTIAL ENERGY) Of all displacements sa- 
tisfying the given boundary conditions those which satisfy the equi- 
librium conditions make the potential energy an absolute minimum. 


6.2.2 Complementary energy 


Now we proceed to prove another important minimum theorem. 
As usual, see also figure 5.1 on page 103, let V be a body in an 
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equilibrium state under the volume forces b; and surface forces f; 
prescribed on Sc. 
The elastic problem is given as follows 


divo +b=0  VpeV (6.66) 
e=symVu VpevV (6.67) 
E v 
‘di pray (e + E ue) VpevV (6.68) 
and 
on=f Wes, (6.69) 
(6.70) 


We term U* conjugate potential energy of deformation or com- 
plementary energy and we shall show that this potential reaches an 
absolute minimum value when the displacements of the body V are 
those of the equilibrium configuration. We define the complementary 
energy of deformation U* as follows 


U" (Tij, fi) = B* (oij) — V(f) (6.71) 

where 
®* (aij) = fe (aij) dV (6.72) 
W* (fi) = f fiuidS (6.73) 


We suppose that over the remaining part S„ the displacements 
a; are known. We denote the displacements which satisfy the equi- 
librium configuration as u; and consider an arbitrary small varia- 
tions of o and f: éo and df, respectively. It is required that the 
fields o + do and f + ôf assures the equilibrium state, so that, by 
virtue of equations (6.22) and (6.25), we easily obtain 


O04; = 0 in Y (6.74) 
do;;ni = 0 on Sc (6.75) 


The first variation of U* is 
su* =U*(f+6f,0+ 60) - U*(f,0) .76) 
=" (o +60) — D* (o) — W* (F) + W* (f+5f) (6.77) 
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which in components becomes 


du* (Sij) = ğ* (Sij) +f CishkThkÒ0;j (6.78) 
y 
—* (Sij) = i. ô fiuidS (6.79) 
S 
finally, 
6u* (Sij) = 1 d0;jEij a l ô fiuidS (6.80) 
V S 


where we recognize on the right-hand side of the above expression 
the Principle of Virtual Work (5.11). So it is easy to notice that 
the stationary point for complementary energy U/* corresponds to 
the equilibrium condition. We know, in fact, that the PVW and 
the compatibility conditions of displacements du; lead to the same 
result: the equilibrium state. 


Theorem 5 (STATIONARY VALUE OF COMPLEMENTARY ENERGY) The 
total complementary energy U* of an elastic body has a stationary 
value in the class of the statically permissible state of stress for the 
true state of stress corresponding to the equilibrium. 


It is also possible to prove the following stronger theorem. 


Theorem 6 (MINIMUM COMPLEMENTARY ENERGY) The complemen- 
tary energy U* has an absolute minimum when the stress tensor cij 
is that of the equilibrium state and fulfills the conditions (6.74) and 
(6.75). 


6.2.3 Theorems of Castigliano 


Results discussed in the previous sections give us the means to 
find some other important results which go by the name of theorems 
of Castigliano. Therefore, let us suppose that a body V is subjected 
only to concentrated loads F} as figure 6.4 shows. 

The potential energy and the complementary energy are, respec- 
tively 


© (eij) LR ui (6.81) 


U* = &* (0;) DI Fru; (6.82) 
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Figure 6.4: Concentrated loads acting on the body V. 


Rearranging the scalar product on the right-hand side of equa- 
tions (6.81) and (6.82) we obtain 
U = ® (eij) — Fhux (6.83) 
u* = d* (Sij) = Fru, (6.84) 
where uz is the component of the displacement vector at the point 
of application of F* in the direction of this force. Now we suppose 
that the potential energy U results only from the displacements ug 
and the complementary energy U* results only from the external 
concentred loads FF, therefore the energies above can be written as 
follows 
U (ui; cay Un) = PU Pu (6.85) 
CPP) (6.86) 
where the first variations are 
OU = ® (u1, Uk + Sup, un) — Ð (u1, Uk, Un) — Fou, = 


2 SO UE si — FE6u = (Se > r") Sup (6.87) 
k 


D* (FF y* (FF 
=" a ) 5p Fup = (FE - SFE (6.88) 
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and by means of the principles of stationary value, equations (6.87) 
and (6.88) vanish, and discarding the trivial solutions ĝu! = 0 and 
6FF = 0, we obtain 


dd, 


OP _ pk (6.90) 
Our 


where F* is a generic concentred load and uy is its corresponding 
displacement. 
Equation (6.90) may also be enounced in the following form 


Theorem 7 (CASTIGLIANO’S THEOREM) If an elastic body is sub- 
jected to concentred loads and supported in such way that each rigid 
body motion is inhibited, then the displacement component ux of the 
point of application of F? towards its direction, is obtained from the 
partial derivative of the complementary energy with respect to the 
particular force. 


Menabrea’s theorem is a particular case of Castigliano’s theo- 
rem. In fact, if F* is a reaction due to a constraint which does not 
allow any displacement, then equation (6.89) will become 


om 2h (6.91) 


Equation (6.90) represents the second Castigliano’s theorem, 
which states 


Theorem 8 (II CASTIGLIANO’S THEOREM) If an elastic body is sub- 
jected to concentred displacements uz and supported in such a way 
that each rigid body motion is inhibited, then the component of a 
concentrated load acting in the direction of such displacement is ob- 
tained by the partial derivative of the potential energy with respect 
to the particular displacement component. 


Illustrative example 


Let us consider the simple system showed in figure 6.5. We have 
a rigid body which is supported by means of two elastic devices. 

The entire elasticity of the system is concentrated in A and B, 
where we have a rotational stiffness with spring modulus km and a 
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Figure 6.5: Example of Castigliano’s theorems. 


translational stiffness with spring modulus k;, respectively. So, the 
constitutive relationships are 


M = kn (6.92) 
R= kyl (6.93) 


and we may also set the following geometric relation 
1 
= 210 
meg 
According to equations (6.5) and (6.11), through the above con- 


stitutive relationships, we can compute the strain energy ® as fol- 
lows 
1 1 1 1 k 
® = MV + RIV = km + Ekl? = 20° (2 +k) (6.94 
gD gra Tg ci ves Ea 


Now, by virtue of second theorem of Castigliano, we can find 
the external load 


so that 


Chapter 7 
Strength of materials 


This section is intended to give only an overview about some selected 
criteria aimed at determining whether the state of stress characterizing 
an elastic continuum is secure compared with conventional material limits 
derived from experimental tests. 

For a thorough investigation on the strength of materials the reader 


is recommended to referred to [7], [9], [10], [12]. 


7.1 Introduction 


Usually the mechanical properties of materials are investigated 
by simple experimental tests, for example the tensile test and com- 
pression test offer the two failure stresses oh and of which allow 
the evaluation of the riskiness of a combined state of stress result- 
ing from the solution of the linear-elastic problem. 

When the real state of stress is simply tension or compression, 
we can directly compare the results with the experimental values 0), 
and og and easily evaluate if the working stress is lower or higher 
than the yielding point (or rupture point for fragile materials). 

The problem becomes more complicated when the actual state 
of stress is combined, hence various theories have been developed 
in order to find laws which, from the behavior of the materials in 
simple compression or tension, predict the condition of failure under 
any kind of combined stresses. The following sketch is rappresenta- 
tive of this statement 


analytical How to compare? experimental 
{ij} == {00:70 
Ci 
It is known that the whole state of stress is defined by three 


principal stresses, so according to equation (3.36) in chapter 3, we 
assume that I3 # 0. Furthermore, suppose the eigenvalue problem 


Claudio Borri, Michele Betti, Enzo Marino (a cura di), Lectures on Solid Mechanics, ISBN 978-88-8453- 
854-3 (online), ISBN 978-88-8453-853-6 (print), © 2008 Firenze University Press 
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(3.35) yields the three principal stresses ordered as follows 
Or > OI > OIN 


where tensions are taken positive and compressions negative. 
The law we are looking for assumes the following general form 


f (01, 011, 0111) = constant (7.1) 


where f (o7,077, 0777) is a comparable quantity and the constant 
can be found applying the criterion f to the simple state of stress 


f (0,0, 00). 
7.2 Maximum stress theory 


This theory, also called Rankine’s criterion, assumes the maxi- 
mum stress as the criterion for the material failure. Accordingly we 
write 


Cia = |o1| (7.2) 
ol = lorrr] 7.3 


where ojqg is the ideal stress, that is the comparable stress. 

In case of ductile materials the theory assumes that the yielding 
starts when the maximum stress becomes equal to the yield point 
stress of the material in simple tension or when the minimum stress 
becomes equal to the yield point stress of the material in simple 
compression. Namely, the failure conditions are 


lor] = 99 
lorr = o0 7.5 


Whereas, the safety side is ensured by the following conditions 


lor| < 00 (7.6) 
lorrr < ob (7.7) 


This theory is not comprehensive and presents some limits. Con- 
sider a specimen under simple tension, sliding occurs along the plane 
where the stress does not attain the maximum value. Moreover, con- 
sider an homogeneous isotropic material weak in simple compres- 
sion, it can sustain very large hydrostatic pressure without yielding. 
That proves the magnitude of the maximum tensile or compressive 
stress alone does not define the yielding condition. 
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7.3 Maximum strain theory 


In this theory, historically attributed to Grashof, it is assumed 
that the yielding of a ductile material starts when either the maxi- 
mum strain, i.e. elongation, equals the strain 0j/E at which yield- 
ing occurs in the simple tension or when the minimum strain, i.e. 
compressive strain, equals the strain of /E at which the yielding 
occurs in simple compression. Therefore, by recalling constitutive 
equations (4.25) and (4.26) on page 95, the criterion is stated as 
follows 


[Emin] = 50 7.9 
that is 
Li tit (7.10) 
Emax = E OI V\OII OTIT)) = E 2 
pa Ade (7.11) 
Emin = E OIII_V\O0ItT0II}))= E : 


keeping the order or < orr < orrr. The ideal stresses to be com- 
pared to the experimental values of and of are 


Cig = 01 — V (orr + orr) (7.12) 
Og = On — Y (or + orr) (7.13) 


The failure conditions are 


lor — v (orr +0111) | = 00 (7.14) 
lorrr — v (or + orr) | = 95 (7.15) 


where o) and of are, as said before, the failure point stresses in 
tension and compression, respectively. The safety side is ensured by 
the following conditions 


lor —v (orr + orr) | < oo (7.16) 


lorrr — v (or + orr) | < 0G 7.17) 
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7.4 Beltrami’s theory 


For the first time in 1885 Beltrami proposed an energetic ap- 
proach to make comparable a combined state of stress with a simple 
tension or compression state of stress. Indeed, the quantity of strain 
energy stored per unit of volume of the material has been assumed 
as a basis to define the stresses at which the yielding starts. 

The elastic potential or the volume density of strain energy, as 
defined in chapter 6 (equation (6.5) on page 112) for the principal 
state of stress and strain assumes the following form 


$ 


= Sf (of +07, + oirr — 2v (01011 + oror + 0110111)) (7.18) 

Beltrami stated that the failure of a body under a combined 
state of stress occurs when its volume density of strain energy equals 
the elastic potential at the yielding point for a simple tension, so 
that 


o=— (7.19) 
where it is assumed that the material has the same behavior both 
in compression and tension, i.e. of = 0) = Co. 


Therefore, taking into account the previous expression for ¢, we 
obtain 


oF + oF, + ofp, — W (01011 + 010111 + 0110111) = 95 (7.20) 


Finally, the linear elastic behavior of the material is ensured 
when the comparable stress giq, given as follows 


Cid = +/03 + ofr + of — 2v (01011 + ron + orror) (7.21) 
is lower than the yielding point stress. So we have 


[Sial < 00 (7.22) 


7.5 Von Mises’ criterion 


Following Beltrami’s approach, in 1913 R. von Mises proposed 
a new method to evaluate the failure state for materials. In this 
theory not the whole elastic potential is assumed as responsible 
for yielding but only the potential energy ®p due to the deviator 
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stresses. Namely, Mises’ criterion assumes negligible the hydrostatic 
state of stress to evaluate when the yielding starts. 

Recalling equation (3.82) on page 78, any state of stress and 
strain can be split as follows 


Cij = ombij + Sij (7.23) 
Eij = embij + ej; (7.24) 
where om and em are the spherical state of stress and strain re- 
spectively, while s;; and e;; are the deviator stresses and strains, 


respectively. Hence, the total potential energy can be written as 
follows 


1 1 
® (eij) = 50E = 3 (OmMbij + sij) (Emdi; + eij) = 
1 3 
38ijCij + 30MEM (7.25) 
2 2 
Let us now define 3 
Pu = gMEM (7.26) 


as the spherical energy, that is the potential energy associated to 
the volumetric variation of the body and 


1 
Pp = g Siei (7.27) 


as the deviator energy, that is the potential energy associated to 
the shape variation of the body. 
As assumed in Beltrami’s criterion, the material is supposed to 
have the same behavior both in simple tension and compression. 
Now, evaluating the deviator energy depending only on the 
stresses, we have 
1 Sij 1 
alia  4G II 
The limit value for yielding is obtained by applying the criterion 
to the simple tension oo which can also be split as follows 


Bp = (7.28) 


00 0 00/3. 0 0 
00 O |= 0 00/3 0 + 
0 0 00 0 0 00/3 
—00/3 0 0 
0 -00/3 0 (7.29) 


0 0 200/3 
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consequently, the limit value for the deviator energy results 


12 
dI = —50) 7.30 
D = 146300 ( ) 
Finally, imposing the criterion (7.1), that here becomes 
1 1 2 
0 2 
the ideal stress assumes the following expression 
3 
Tid = \f 3Si98% (7.32) 


The above expression can also be given through the principal 


stresses as follows 
3 
Cid = g i ij = 
32 2 2 2 2 diva 
5 (si, + 539 + 83 + 251, + 2573 + 2853) = 
3/2 2 2 \_ 
2 (si + 57, + sîr) = 


NÈ (Cor — om} + (orr — om)? + (or — om?) = 


2 2 2 
Vof +07, + 0%, — 0101 onor — 010111 (7.33) 


7.6 Criteria comparison 


To highlight out the differences among the four methods above 
discussed let us consider the graphical interpretation for each of 
them. For the sake of simplicity consider the case where crr = 
0. Thus, the whole state of stress is given by the two principal 
stresses {0,, 0,1} not necessarily sorted as equation (7.1) shows. 
Moreover, we shall assume the material has the same behavior in 
simple tension and compression. 


7.6.1 Maximum stress 


By the preceding assumptions the comparable stress in this case 
turns into 
Cia = max{|o7|,|077|} (7.34) 
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and the rupture limits are 


Or = 100 (7.35) 
OIL = E00 7.36) 


that on a Cartesian plane depict a square domain. See figure 7.1. 


A Or 
2 BZ 1 
6 0, A 
— 00 To fon, 
3 B?— o 4 


Figure 7.1: Rupture domain for the maximum stress criterion. 


The lines in the figure represent the values of oy and 077 at which 
yielding starts. The lengths OA and OB represent the yield points 
in simple tension along the directions of a7 and 077, respectively. 
In the same way A’ and B’ represent the yielding points for simple 
compression. Moreover, the four conditions above ensure that any 
point within the square 1234 define an elastic configuration. Hence, 
we can define lines 1234 as rupture (or yielding, for ductile material) 
boundaries. 


7.6.2 Maximum strain 
Here the ideal stress assumes the form 
Cia = max{|o7 — vorr), lorr — vor|,|-v(or+or)|} (7.37) 


and the yielding limits are 


or vor] = £00 (7.38) 
OJI — VOI = x09 (7.39) 
—v (or + orr) = +00 (7.40) 
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— Colv 


Figure 7.2: Rupture domain for the maximum strain criterion. 


that on a Cartesian plane form a domain as shown in figure 7.2. 

Note that the drawing in figure 7.2 has only an illustrative in- 
tention, in fact, in the reality, points 5— 6 and 8—9 are much closer 
each other when v approaches values around 0.3. 

Figure 7.2 also shows that if two principal stresses are equal 
and opposite in sign, the maximum strain theory indicates that the 
yielding starts at a lower value than the maximum stress theory 
would indicate, see points 4 and 10, for instance. 

On the other hand, since a tension in one direction reduces the 
strain in the perpendicular direction, two equal tension can have 
higher values at yielding than the maximum stress theory. 


7.6.3 Beltrami’s criterion 


Here the ideal stress assumes the form 


Cid = Joi + o? — 2Qvororr (7.41) 
and the yielding boundary is described by the ellipse 
o? + 02; — 2vororI = og (7.42) 
that intersects the orz and oy; axes at points 


or = +00 (7.43) 
OII = 00 (7.44) 
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A Or 


Figure 7.3: Elastic domain for Beltrami’s criterion. 


The intersections with the bisector orr = 07, passing through 
the first and third quadrants, occurs at points A and B, see figure 
7.4, which have coordinates respectively 


00 
dale ee 7.45 
arc 2(1—v) CER, 


while the intersections of the ellipse with the bisector orr = —07, 
passing through the second and fourth quadrants, occurs at C and 
D having the following coordinates 


DEE. shat ee OE 
Gey. pg ey 
PR EI phi ca (7.47) 

V2(1+v) V2(1+v) 


respectively. See figure 7.3. 


7.6.4 Von Mises’ criterion 


From equation (7.33) we derive the expression of the ideal stress 
when a two-dimensional state of stress occurs 


Cid = lo} +0}, — oror (7.48) 


hence, the yielding limit is found when Gia = 09 so that 


o0 = VÈ +03, - oon (7.49) 
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The latter represents in the (07, 017)-plane the following ellipse 
05 = 07 + OF; — 1011 (7.50) 


having as the major axis the bisector of the first and third quadrants 
and as minor axis the bisector of the second and fourth quadrant. 


401 


Figure 7.4: Elastic domain for the Mises’ criterion. 


The intersections with bisector line orr = 07 occurs in points E 
and F, see figure 7.4, which have coordinates 


o]= 011 = too (7.51) 


while the intersection points with the bisector orr = —07 are G and 
H which have coordinates 


90 90 
op=—-— of =+ 7.52 
I V3 II V3 ( ) 

00 00 
op=+— -—or=—> 7.53 
Ù V3 A V3 l ) 


respectively. See figure 7.4. 


7.6.5 Comparison 


In order to compare all the above criteria let us consider a two- 
dimensional state of stress defined as follows 


Tij = ( de, dp ) (7.54) 
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It is straightforward to represent this stress state through Mohr’s 
circle where we find 


xc=0 
yc=0 
R= 012 


and the coordinates of the pole in the (a, 7)-plane are P* = (0, — 012). 
For further detail on Mohr’s circles see section 3.4.2. 

Moreover, the principal stresses located at points S4 and S2, see 
figure 3.10, are respectively 


OI = —OJI = O12 


Thus, if we define the following general expression for the yield- 
ing point 
012 = 00 
then for each criterion we can readily evaluate the coefficient a 


considering the equations above discussed. For the sake of clarity 
the relevant results are collected in table 7.1. 


CRITERION a a (v = 0.3) 
Maximum Stress 1 1 
Maximum Strain iw 0.77 

Beltrami’s asta 0.62 

2(1+v) 
: ’ 1 
Mises Z 0.58 


Table 7.1: Criteria comparison. 


PART II 


Theory of elastic beams 


Claudio Borri, Michele Betti, Enzo Marino (a cura di), Lectures on Solid Mechanics, ISBN 978-88-8453- 
854-3 (online), ISBN 978-88-8453-853-6 (print), © 2008 Firenze University Press 


Chapter 8 


Saint-Venant’s problem 


This section is devoted to the theory of the beam. We will first in- 
troduce some necessary adjustments of the three-dimensional theory of 
elasticity in order to provide an ad hoc mathematical model for prismatic 
structural elements. Then the mechanical behavior of the beam will be 


analyzed considering separately four fundamental cases. 


8.1 Statement of the problem 


The solution of the general boundary—value problem presented 
in section 4.2 often presents some mathematical difficulties because 
of the complicated form of the boundary conditions. Frequently it 
is necessary to introduce some simplifications in order to ensure 
solutions for technological applications of the theory of elasticity, 
so that the mathematical solutions of the problem represents only 
an approximation to the actual situation. 


In order to simplify the boundary conditions let us assume the 
following principle on which the theory of beams is founded. 


If some distribution of forces acting on a portion of a 
body is replaced by a different distribution of forces act- 
ing on the same portion of the body, then the effects of 
two different distributions on the parts of the body suffi- 
ciently far from the region of application of the forces are 
essentially the same, provided that the two distributions 
of forces are statically equivalent. 


where “statically” equivalent means that two distributions of forces 
have the same resultant force and the same resultant moment. This 


Claudio Borri, Michele Betti, Enzo Marino (a cura di), Lectures on Solid Mechanics, ISBN 978-88-8453- 
854-3 (online), ISBN 978-88-8453-853-6 (print), © 2008 Firenze University Press 


144 SAINT-VENANT'S PROBLEM 


principle was proposed in 1885 by J. C. B. de Saint Venant!. 


First of all we declare the fundamental hypothesis which define 
a method of solution. 


Shape of solid. We define a beam as a particular body bounded 
by acylindrical surface called the lateral surface and by a pair 
of planes normal to the lateral surface called the bases of the 
cylinder. We shall also suppose the cross section is constant 
and the beam’s length is much larger than the cross section’s 
linear dimension: l > r, see figure 8.1. The Cartesian coordi- 
nate system is positioned having the x3-axis taken along the 
length of beam and parallel to the lateral surface. This axis 
usually coincides with the axis of the beam passing through 
the centers of gravity of the bases. The cylinder is assumed to 
be of length l so that one of its bases belongs to the (x1, £2)- 
plane and the other is taken at x3 = l. 


Loads. It is supposed that the lateral surface of the cylinder is 
load free and that the loads act only on its bases x3 = 0 and 
x3 = l. Moreover, the forces at the ends assure the equilibrium 
condition of the cylinder. It is also supposed that the body 
forces b are zero. 


Constraints. According to the previous point, we shall suppose 
that the cylinder is unconstrained, and the forces acting on 
the bases fulfill the global equilibrium equations. However, to 
assure that no rigid displacement is allowed, at least one point 
of the beam, say G, must be fixed. So we will assume that to 
inhibit any translation, at x1 = £2 = x3 = 0, the following 
constraints hold 


ui = Ug = u3 = 0 (8.1) 


'Adhémar Jean Claude Barré de Saint-Venant (August 23, 1797 Seine-et— 
Marne - January 6, 1886 Seine-et-Marne) was a French engineer. 


Source: http://www-history.mcs.st-andrews.ac.uk/PictDisplay/Saint-Venant.html. 
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Figure 8.1: Prototype of beam. 


and to avoid any rotation we require that 
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Material. Let us assume an homogeneous isotropic linear elastic 


material. 


State of stress. We shall assume the stress vector normal to the 
lateral surface is zezo. So, for such a unit vector n = n°; = 


1 


Onn = Cijn n’ = oynin! + oan n? + 2019n!n? = 0 (8 


only if 
(n°)? + (n?) =1 


nle, + n%é normal to the lateral surface, in accordance with 
the location of the coordinate system, we have 


5) 


(8.6) 


Hence, equations (8.5) and (8.6) induce? a particular form of 


the stress tensor 


0 0 013 
Tij = 0 0 023 
031 032 933 


It can be immediately proved putting first n! = 1 and n? = 0, then n! = 0 


and n? = 1. 
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The equilibrium problem of an elastic beam with a free lateral 
surface subjected to loads only on the bases can be formulated as 
follows: determine the three components of stress 013,023,033 and 
the displacements u; that satisfy equations 


Tij j =0 on V 


By virtue of stresses and loads assumptions, the above problem 
becomes 


013,3 = 0 
0233 = 0 on Y 


031,1 + 032,2 + 033,3 = 0 


Next, in accordance with the hypotheses, the boundary equa- 
tions are 


TijNj = fi on £3 = 0, T3 = l (8.10) 


that can be expanded as follows 


T3 = 0, n= (0,0, —1) 


013Nn3 = —013 = fi (8.11) 

033n3 = —023 = fa (8.12) 

0333 = —033 = fa (8.13) 
T3 = n= (0, 0, 1) 

013N3 = 013 = fi (8.14) 

0933 = 023 = fo (8.15) 


033N3 = 033 = fs (8.16) 
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/ ) 


(a) Basis x3 = 0 


(b) Basis 73 = l 


Figure 8.2: Unit normal vectors on the bases of the cylinder. 


Constitutive equations (4.25) become 
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1 
€11 = ~ 7033 (8.17) 
1 
E22 = ~ 7033 (8.18) 
1 
€33 = Zz (lity) 033 — V033) = È (8.19) 
€12 = 0 (8.20) 
l+v 
gn E la (8.21) 
l+v 
€23 = TE log (8.22) 


Now, making use of Beltrami-Michell’s equations (4.44), it is 
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possible to write 


033,11 =0 (8.23) 
033,12 =0 (8.24) 
033,22 =0 (8.25) 
=0 8.26 
033,11 + 033,22 + 033,33 + I+ 7,733,38 ( ) 
023,11 + 023,22 + 023,33 + 033,23 = 0 (8.27) 
l+v 
1 
— =0 8.28 
013,11 + 013,22 + 013,33 + 1% 773318 ( ) 


and by means of equations (8.7) and (8.8), the above expressions 
turn into 


033,33 + 7 103383 = 0 = 033,33 = 0 (8.29) 
= (J 8.30 

023,11 + 023,22 + 1% 77 33:28 ( ) 

= 0 8.31 

013,11 + 013,22 + T+ 773318 ( ) 


Equations (8.23), (8.25) and (8.29) suggest that the compo- 
nent 033 must vary linearly with 21,272,273, while equation (8.24) 
imposes that it cannot contain the product x1x2. Hence, 033 = 
033 (£1, £2, £3) must assume the following form 


033 = a + bzı + cx2 — (d + ex1 + fx2)x3 (8.32) 


See also [9] and [11]. 


8.1.1 External and internal forces 


In this section we point out that in most of the practical circum- 
stances we know the resultant force F and the resultant moment M 
acting on the ends of a beam rather than the real external surface 
force distribution f. Indeed, if we accept the Saint Venant’s prin- 
ciple, we could not care about the nature of the stress distribution 
which produces the resultants F and M, just because we are inter- 
ested in portions of beam sufficiently far from the ends where the 
influences of boundary stress distribution is not decisive. 
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Let us fix our attention on the bases of the cylinder. Firstly, 
on the basis 73 = 0, denoted by Ao, we assume a distribution of 


surface forces given by f°, next, on the basis z3 = l, denoted by Aj, 
we assume a distribution of surface forces fl. 


For the basis 73 = l the resultant external force are related to 
the external surface forces as follows 


l= | fidi (8.33) 
Ai 

T= | fidA (8.34) 
Ai 

N= | fidA (8.35) 
Ai 


and, of course, we can do the same for the basis x3 = 0 


n= | fedA (8.36) 

Ao 

F? = 1 PdA (8.37) 
Ao 

Ñ? = fOdA (8.38) 
Ao 


where 7} and 109 lying in the plane of the basis, are responsible for 
bending and sharing of the beam, while N, taken in xg-direction, 
is responsible for tension or compression. See figure 8.3. 


The couple M, analogously, may be split into a component M 
along the x3—axis which provides the twisting for the beam and the 
components Mı and My which are responsible for bending. Thus, 
for both bases we have 


M = | fixodA (8.39) 
Ai 


z f flaidA (8.40) 
Al 


=f (- fizz + fin) dA (8.41) 
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and 
wp = | foadA (8.42) 
Ao 
MI = -f feaidA (8.43) 
Ao 
M} = / (- fz + fn) dA (8.44) 
Ao 


where the moments have been computed with respect to the centers 
of gravity for 4; and Ao, respectively. See figure 8.3. 


Figure 8.3: Equilibrated components of force and couple resultants 
acting on the ends of the beam. 


As already mentioned, the readers should be aware that in many 
actual problems we shall know just the components of resultants F 


and M rather than the actual distribution of surface forces f, so, 
usually, we will solve an inverse problem. To this end in the following 
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we want to find how the resultant forces are transmitted inside the 
beam. 

From the end z3 = l, towards the opposite basis x3 = 0, the 
external actions propagate in such a way that, for any cross-section 
taken at the distance x3 from the origin of the coordinate system, 
the components of resultant forces and moments are, respectively 


N+ (x3) = N' (8.45) 
Ti (#3) = Ti (8.46) 
T3 (23) = T5 (8.47) 
and 
My (x3) = Mj — 7} (l — 23) (8.48) 
My (a3) = M3 + Tj (1 — 23) (8.49) 
Mł (z3) = My (8.50) 


where the sign convention has been established by the following 
rule: on the “positive side” of a generic cross section whenever the 
forces orientation and the coordinate axes are concordant, then the 
positive sign is assigned; concerning the components of the couples, 
the positive sign is ascribed whenever the moments are concordant 
with the following rotations: £2 —> £3, £3 > £1, £1 > T2. 

Here the reader will also realize that by virtue of the equilibrium 
condition of the cylinder, taking into account the external actions 
from the opposite side, here crudely named “negative side”, for the 
same generic cross section considered above, the local equilibrium 
condition must be assured. In this view, the sign convention as- 
sumed when the “negative side” of a generic cross section is taken 
into account, is: the external forces have positive sign if they are 
discordant with the coordinate axes and, concerning the moments, 
they will be assumed with a positive sign if they induce the follow- 
ing rotations: £3 — £2, £1 > £3, £2 — T1. In the light of this sign 
convention, we notice that equations (8.45) to (8.50) can be equiva- 
lently written, for the same generic cross section at x3, considering 
the external actions on the “negative part”, that means 


NT (x3) = N° 
T; (z3) =f? 
Tz (x3) = Î} 
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and 
My (a3) = M? + Ty 23 
My (z3) = M3 — Tag 
Mg (x3) = M3 
Naturally, for each cross section the external actions on the 


right-hand side (“positive part”) must equal those on the left-hand 
side (“negative part”)?, indeed we can write 


N (#3) = N* (23) 
TI (#3) = T? (x3) 
Ty (£3) = Ty (x3) 
and 
My (z3) = M} (z3) > M? + T$x3 = Mj — TS (1 — 23) 
My (a3) = My (x3) > M3 — Tag = M} + Tj (l — zs) 
Mz (a3) = Mf (a3) > M$ = M; 


When x3 = 0 the first set of the above equations remains unal- 
tered, while the second, i.e. the bending forces, becomes 


a es m> Di 


In order to pass from M; (0) to ML the signs must be corrected 
in accordance with the SILA sign convention, that is 


T= —-T; (0) =-Tj 
T} = —T (0) = -T} 
and 
MÌ = -Mī (0) =-M{+Tl 
My = —M; (0) = -M, -Til 
M3 = -M; (0) = —M3 


3For this reason hereafter we will omit the primes (-)* or (-)7. 
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Now, imagine to cut the beam at the cross section located at 73. 
The set of forces {T1, T2, N, M1, M2, M3} can be seen as external 
forces, thus, for this section - that is assumed being seen from the 
right-hand side - the well known boundary condition oj; = fi holds. 
This justifies the following expressions 


li 031dA (8.51) 
A 
B= | oszd A (8.52) 
A 
Nel o33dA (8.53) 
A 
and 
m= | 03322dA (8.54) 
A 
Mo = -f 033r1dA (8.55) 
A 
M= | (-031%2 + 0321) dA (8.56) 
A 


8.2 Four fundamental cases 


Now we are ready to present the four fundamental cases that, 
by virtue of the superposition principle, allow us to entirely solve 
the problem of elastic beams. They are 


e Extension of a beam by axial force applied at the ends. 


e Bending of a beam by couples whose moments lie in the plane 
of its bases. 


e Torsion of a beam by a couples whose moment is normal to 
its basis. 


e Flexure of a beam by transverse forces applied at one end of 
the cylinder, while on the other end it is acting an opposite 
transverse force and a couple in such a way the equilibrium 
condition is fulfilled. 
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8.3 Beam under axial force 


8.3.1 State of stress 


Consider the basis x3 = l, and assume that the N° is the only 
nonzero external force 


N' 40 (8.57) 
È = È = M! = M = Ml = 0 (8.58) 


therefore, in accordance with equations (8.45) to (8.50), for a generic 
cross section at x3, the external forces are 


N=N'! (8.59) 
Ti = Ty = Mı = Mp = M3=0 (8.60) 


Making use of equations (8.51) to (8.53) and the boundary equa- 
tions (8.14) to (8.16), the surface forces are 


fs #0 (8.61) 
fi = fo=0 (8.62) 
So, considering expression (8.32) for the normal stress we can 
write 
i I. Liana delie «a 
A 


and recalling equation (8.9), where 031 = fi =0 and 039 = fo = 0, 
the above equation turns into 


N= i) (a+ ba, + cz2) dA = (8.64) 
A 
=a dA+5 f nidA+e | xodA (8.65) 
A A A 
=0 =0 
so we have proved that 
N 
= — 8.66 
a=3 (8.66) 


Equations (8.39), (8.40), through condition (8.58), represent a 
linear system whose solution yields 


b=c=0 
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hence, the expression of the normal stress is finally 


N 
Se 8.67 
33= 5 (8.67) 
and the stress tensor assumes the following form 
0 0 0 
Tij = 0 0 0 
0 0 N/A 


8.3.2 State of strain 


The deformation of a beam under an axial force results directly 
from the constitutive laws (8.17) to (8.22), so that 


al (8.68) 
eq. = -V —— : 
> EA 
a (8.69) 
E22 = —v—— ‘ 
# EA 
N 
eee 8.70 
E33 EA ( ) 
and the strain tensor is 
—vN/EA 0 0 
Eij = 0 —vN/EA 0 
0 0 N/EA 


8.3.3 Displacement field 


In order to compute the displacement field uj = u; (xj) of the 
cylindrical body we have to solve the system of differential equations 
obtained from compatibility relationships. Hence, the equations are 


N 
Wi = €1 = -VEA u12 + u21 = 262 = 0 
u22 = ER = Vey and U13+U31 = 26133 = 0 
N z L 
u33 = €33 = EA u23 + U32 = 2623 = 0 


The integration of the first group of equations gives 


u = -vgt +a (x2, T3) 
uz -vpt + b (z1, £3) 
ug. = PT3 + (21,22) 
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where a, 6,y are unknown functions. Replacing the above expres- 
sions in the second group of differential equations, and assuming 
inhibited any rigid body motion, some calculations drive us to write 
the following solution 


N 
ui = ce gi (8.71) 
N 
See 72 
ug av Ge (8.72) 
N 


Through the above results we are able to know the strained 
shape of a beam subjected to an axial force. Let p be a point inside 
the beam, so that p = (x1, 72,73), we notice that the displacement 
u3 of p does not depend on the position of p in the cross section area, 
so we have constant displacements along x3—axis for each cross sec- 
tion. Let p' = (x, 24, 24) be the position of p after the deformation, 
so that 


£i = £1 + (8.74) 
xh = £3 + ug (8.76) 


we can define unit axial strain as 


(u3 + duz) — u3 _ ON 
dis E33: = EA (8.77) 


that integrated along the entire length of the beam gives the axial 
elongation 


Al= [endes = E331 (8.78) 
l 


To the elongation Al corresponds a cross-section contraction p 
that can be measured by the vector sum of u and us as follows 


P = Ue] + U22 (8.79) 


The magnitude of p is 


vN 
p = Vu? + ug? = BAY L1? + T2? (8.80) 
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di N 


La Ty 


Figure 8.4: Strained state of a beam subjected to an axial force. 


Equation (8.80) shows that the displacement of a point p € A 
depends solely on its distance from the center of gravity of A. We 
say also that each point p moves towards G along radial direction. 
In fact, by considering a polar coordinate system {r, 9} with the 
origin in the centroid G, the related basis is obtained by making 
use of the well known transformations, see page 20 in chapter 1, 


= Ox Ox ee = 

= a + J7 = sin VE, + cos VEz (8.81) 
F Ox Ox 

dy = T é, + T — è, = r cos VE — r sin VEz (8.82) 


hence, to ensure the radial direction, it is enough to prove that 
Ov - p = 0. Indeed we have 


ðv. p= u sind — uz cos = z3 cos -2 + tang (8.83) 
S T2 


but since it is also known that x1/r2 = tan ð, (8.83) always van- 
ishes. See figure 8.5. 
8.3.4 Strain energy 


By virtue of Clapeyron’s theorem, see equation (6.15) on page 
114, the strain energy is 


da 1 N21 
= A= SN [= -—_ 8.84 
2 PES EA on) 
or, in the same way, by using me PLV, we have 


1 1 N°1 
D = z fornendo = 54 [Frate “JEA (8.85) 
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Figure 8.5: Strained state of a beam subjected to an axial force: 
radial contraction. 


8.4 Beam under terminal couples 


8.4.1 Introductive sketch 


Before going ahead it is useful to consider an idealized model of 
beam made up of long filaments parallel to the axis of the cylinder. 
By virtue of Saint Venant’s hypotheses we know that the stresses 
are zero on the lateral surfaces, in the direction perpendicular to 
filaments’ length, and act only on the ends of the filaments. 

Let us consider now a beam subjected to a pair of equilibrated 
couples at the ends. See figure 8.6. Because of the couples, the lower 
longitudinal filaments, those towards the center of curvature, will be 
contracted and the extrados, the upper portion of filaments, will be 
extended. We shall assume that the central line, i.e. the line passing 
through the centroid of all cross sections, is unaltered in length and 
the cross sections lie always in a plane normal to the cental line. 


Figure 8.6: Beam under terminal couples. 


The upper longitudinal filaments, initially parallel to the x3— 
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axis, under the effects of the couples, at a distance f from the 
central line, turn from dso into ds 


ds = (R+ f) dd (8.86) 


where R is the radius of curvature of the central line. Now we are 
interested in evaluating the elongation of a generic material fibre, 
so we define the extension e as 


o ds ds _ (R+ f) dd Rdd _ f 


dso Rdd R 


This linear elongation may be thought to be produced by a 
longitudinal stress 033 which in accordance with equation (8.19) is 


(8.87) 


given by 
=—f (8.88) 
033 = 
3R 


These intuitive considerations allow us to understand that the 
normal stress produced by the flexure of a beam is not constant 
along the cross section and is proportional to the distance from the 
central line. 


8.4.2 State of stress 


In this case we assume that the external forces acting at the 
end of a beam are couples whose moments lie in the plane of the 
cylinder. So that 


M #0, M #0 (8.89) 
N=Tl=Ti= M= 0 (8.90) 


therefore, in accordance with equations (8.45) to (8.50), for a generic 
cross section located at x3 the external forces propagate as follows 


Mı = M! (8.91) 
Mə = M? (8.92) 
N=T, =T = M3 =0 (8.93) 


The boundary conditions (8.14), (8.15), (8.16) on page 146 in 
this case are 

f3 #0 (8.94) 

fi = fo =0 (8.95) 
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Hence, we can write 


f 33dA =0 (8.96) 
A 
Mı = f 033r9dA (8.97) 
A 
Mo = -f 033r1dA (8.98) 
A 
which become 
I (a+ br, + cx2) dA = 0 (8.99) 
A 
Mı = f (a + ba, + cz2£2) dA (8.100) 
A 
Mə = -f (a + ba, + cx2) cd A (8.101) 
A 


By integrating equation (8.99) we find a = 0, therefore equations 
(8.100) and (8.101) become 


Mı = f bzizodA +f cto dA (8.102) 
A A 
Mə = -f bora f cer 1dA (8.103) 
A A 
that can be written in matrix form as follows 
J12 Ji b =! Mi 
Jo J12 (6 T — Mə 


In order to simplify the solution of the above linear system, we 
may rewrite the system with respect to the axes of inertia (£, 7), so 


it becomes 
0 Je b\ | Me 
h 0 ce) \ —M, 


that resolved for b, c, gives 


Me 
ano 8.104 
sd ( ) 
M, 
bos (8.105) 
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and finally, the stress induced by the couples is 
My Me 

=——- > 8.106 

033 I, E+ n ( ) 

By means of equation (8.106) we are able to find the line of zero 

stress below of which all the fibers are in tension and above of which 

they are in compression. This line is termed neutral azis. We shall 

denote it by n — n. See figure 8.7. In the (€,7)—plane, the equation 


of the neutral axis is provided as follows 


My Me 
033 noT Je" ( ) 
Mn (2) 
n= (È) e 8.108 
Me \ Pn l ) 


where we have set Jẹ = peA and Jy = pr At. Now if 7 is the angle 
between the axes (€,n — n), we can put 


(8.109) 


Figure 8.7: Projection of the couples and rotation axis. 


da (2) (8.110) 


tando Pn 


‘Note that pg and p, are the radius of gyration taken along n and € respec- 
tively. They satisfy the equation &/p7 + n° | pe =a; 
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Furthermore, let w be the angle between the axes (€,s — s), see 
figure 8.7 on the preceding page, we notice that 5 +ô = w so we 
can write 


cotw = cot (5 +6) = — tanò (8.111) 
while equation (8.110) becomes 
pe\? 
tan ytanw = — (=) (8.112) 
Pn 


The latter equation allows us to find the neutral axis given any 
bending pair Mı and Mə. In fact, we can summarize the procedure 
to find the neutral axis in the following items 


- find the axis s — s, orthogonal to the bending vector M = 
Mie + M9€3; 


- find the principal axes of inertia, i.e. rotate the system {x1, £2} 
into {€,7} and compute the radii of gyration; 


- w is now known; 


- y, i.e. the position of the neutral axis, can be computed by 
equation (8.112). 


The axis (f — f) orthogonal to the neutral axis n — n is called 
flexural axis and the angle it forms with the axis s — s gives a 
measure of the bending deviation. 


8.4.3 State of strain 


The state of strain results directly from constitutive laws (8.17) 
to (8.22). Hence, chosen the system of principal axes (£, n, 73), with 
respect to the origin G, we can write 


V 
Ege = — 733 (8.113) 
Em = -5033 (8.114) 
T33 
=v 8.115 
€33 E ( ) 


A simpler solution is obtained by introducing a new coordinate 
system {x, y, z} where, as showed on figure 8.8, x is the neutral axis, 
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y is the flexural axis, and z = x3. The advantage of this coordinate 
system is that the normal stress can be written as 


Ozz = MY (8.116) 


and the constant m is readily computed by 


M, = I 033YdA = m | yY dA = mJy (8.117) 
A A 
hence 
Ma (8.118) 
ma — 3 
Jz 


Finally, equation (8.106) assumes the following monomial ex- 
pression 


(8.119) 


Figure 8.8: Neutral axis and flexural axis. 


Equation (8.119) is also known as Navier formula. See also [9] 
and |11] for a proof based on geometric considerations. 
The state of strain associated with the following state of stress 


0 0 0 
Tij = 0 0 0 (8.120) 
0 0 (Mz/Jz)y 


assumes the form 


Eij = 0 —vky 0 (8.121) 
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where 
k= —— (8.122) 


x has an important geometric interpretation which will be discussed 
in the next paragraph. 

Now, recalling equations (2.78) and (2.79) on page 49, it is pos- 
sible to prove that 


AA= i) (cos + Eyy) dA = -ovr f ydA=0 (8.123) 
A A - 
AV = f (Ex + Eyy + Ezz) dV = [af Ky (1—2v)dA 
y l A 
(1 — 2v) if ydA =0 (8.124) 
= 0 


So we notice that throughout the deformation the initial volume 
and the area of every cross-section is unaltered. 


8.4.4 Displacement field 


As the foregoing case, the field of displacement u can be found 
by integrating the system of differential equations (2.70) on page 
48. Let us rename the displacement components with respect to the 


{x,y,z} coordinate system, as u1 = u, ug = v, uz = w. Hence, the 
compatibility equations are 


Ue = Exe = —VRKY Uy tUr = 2exy = 0 
Uy = Eyy = —VKY and UztWa = 2p, = 0 
Wz = Ex, = KY UzgtWy = 2e = 0 


By integrating the first group of the above differential equations 
we find 


u = —vKyxt+a(y,z) 
2 

v = —vK + B(z,z,) 

w= kyz+97(2,y,) 


that replaced into the second group of differential equations yield 


—vKzt+aly,z),+8(2,z), = 0 
a(y,z),+7(2,y) x 
B(z,z),+Kz+7(z,y), = 0 


Il 
o 
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The first of the above equations can be derived with respect to 
y and then to x; the second with respect to z and then to x; finally 
the third with respect to z, then to y. So we obtain 


where equations (8.125) and (8.127) yield 
a (y,z) = A+ Ay + A"z + A"yz (8.131) 


equations (8.126) and (8.129) yield 


2 2 


b (x,z)= UR — KS + B + B'x + B"z + BY xz (8.132) 


while equations (8.128) and (8.130) yield 
y(z,y) =C +C'zr + C"y+ Czy (8.133) 


where A, A’, A”, A”, B, B', B", B", C, C", C”, C" are unknown con- 
stants. 

Equations (8.131), (8.132), (8.133) can be replaced into the sec- 
ond group of the initial differential equations to obtain the following 
system 

A' + DB! ye (A” + B"') z 0 
A" + C' + (A” + Cc") y 0 
CO" + B" + (B" + 6") x = 0 


from which we can derive the solution 


A! + DB! = 0 A! + Bu = 0 
A! + C! ta 0 and A! SE Cl” = 0 
C" + B" = 0 Bl + GU = 0 


where the left-hand group imposes the conditions that A’ = —B’ 
A” =—C' and C" = 
A! = B" = GU! = 0. 


bi 


—B", while the right-hand group assures that 
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In the light of the latter results and making use of the boundary 
conditions which inhibit any rigid body motion, see equations (8.1) 
and (8.3), we finally obtain the components of the displacement 
field related to a beam under terminal couples 


U = —VKyx (8.134) 
v= 3; (v (y? — 2?) + 2?) (8.135) 
w = KYz (8.136) 


Equations (8.134), (8.135) and (8.136) show that the filaments 
lying in the neutral plane, i.e. y = 0, do not suffer any extension. 
The longitudinal material fibres on the side of y > 0 are extended, 
whereas the filaments on the side y < 0 are contracted. 

Now we are able to know the strained shape of the beam. Let 
p and p’ be the positions of a point within the beam before and 
after the deformation, respectively. Hence, the coordinates of such 
positions are p = (x1, £2, £3) and p' = (£1, 74,74). By virtue of the 
above displacement field, we can relate the initial coordinate to the 
strained one as follows 


x =x+u=%x-VKIY (8.137) 
K 

y=y+v=y 3 (2? — v (a? —y’)) (8.138) 

z =z+w= z+ kyz (8.139) 

Focusing on the central line x = y = 0 the above equations 

become 

x =0 (8.140) 

y = i (8.141) 

"At (8.142) 


where we notice that the points on the central line, after the defor- 
mation, go into the points 


yf = -x2 (8.143) 


that describes a parabola whose radius of curvature is given by the 
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following formula 


1 Eg] 
= dz!2 
R” E (8.144) 
which can be approximated by 
1 dy 
— S| 8.145 
R Per ( ) 
3 dy' 2 
assuming (45) to be small. 
Thus, equation (8.145) leads us to write 
1 |Ma| 
== |r| = 8.146 
TEE- (8.146) 


where the constant EJ, is termed modulus of flexural rigidity. 
Points belonging to the central line, i.e. x = y = 0, are subjected 
to the following displacements 


u=0 (8.147) 
v= 3? (8.148) 
w=0 (8.149) 


where in this case v is called the elastic curve (or deflection line) 
and describes the plane curve, i.e. a parabola, that the center line 
assumes when the beam is subjected to pure bending. 


8.4.5 Strain energy 


We have three equivalent tools to compute the strain energy. 


Work done by the external forces. By virtue of the Clapey- 
ron’s theorem, see equation (6.15) on page 114, the strain 
energy is 


1 
® = Meda (8.150) 


where Mi is the only external force and gl, is the rotation in 
the (y, z)-plane at the point of application of the couple, i.e. 
at the end z = l, with a = y = 0. See figure 8.6 on page 158. 
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The rotation È. is the first derivative of the displacement v 
and the sign is assumed through the right-hand rule, so that 


l 
= vv. = kl = l 8.151 
hence, 
ee 
1 (i) I 
De ~~_4_ .152 
2 EJ (9192) 


As figure 8.6 shows, both Mi and Yy are negative with respect 
to the Cartesian axes. 


Work done by the internal stresses. By recalling the equation 
6.5 on chapter 6, we can set 


1 
b= 5 J ousia (8.153) 
2 Jy 


that in the current application becomes 


1 fM l M,\? 
o=-_]/ — ue riti 2 
5 [spun sa (32) ra 
l oan o M21 
=> (=) Jo dA ET. (8.154) 


where we have just used the tensors (8.120) and (8.121). 


Work done by the internal forces. Note that here the sign con- 
vention is taken in accordance with that assumed in section 
8.1.1, so we have for a generic cross section 


d® = M, (z) dp, = Mz (z) kdz (8.155) 


hence, by integrating along the entire beam we find the fol- 
lowing expression 


1 f! 1 M21 
®=- | M = 0 1 
5 f (z) kdz DEI, (8.156) 


that due to the well known relation for forces transmitted 
along the beam, it is nothing but equation (8.152). 


LECTURES ON SOLID MECHANICS 169 


In equation (8.155) we just used the relationship between the 
curvature and the rotation of the section, in fact 


dex 


dz 
that stems form 


_dv(2) 


Pa = dz ld = KZ 


8.5 Beam under torsional couples 


8.5.1 Circular bar 


To capture the basic ideas on the torsional problem let us start 
from the simple case of a circular bar with one end fixed in the 
plane (21,2). At the other end, i.e. x3 = l, suppose there to be 
applied a torsional couple lying around the x3—axis. 


Figure 8.9: Circular bar under torsional couples. 


Under the hypothesis that all cross-sections parallel to the plane 
(x1, £2) remain plane, we can intuitively assume that the magnitude 
of the rotation in a generic section perpendicular to the x3-axis 
depends proportionally solely on the distance from the fixed end, 
see figure 8.9. Such as 

Ù = kaz (8.157) 


where k is the twist rotation per unit of length, i.e. the relative an- 
gular displacement of a pair of cross-sections that are unit distance 
apart. 

Consider now a generic cross-section, as figure 8.10 shows. 

The hypothesis of plane sections means that u3 (p) =0 Vp of 
the bar, moreover the circular shape ensures that a generic point 
p lying on a cross-section can just rotate keeping unaltered the 
distance r from the origin. 
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M3 
LA 


yT 


Figure 8.10: Rotation of a point p lying on a generic cross-section 
of the circular beam. 


Hence, with respect to figure 8.10, we can write 


u = — (r cos 8 — r cos (V + 8)) (8.158) 
uz = rsin (8 + ð) — rsin 8 (8.159) 


that by means of equation (8.157) and considering that 


zı =r cos (8.160) 
vg = rsin ZB (8.161) 


becomes 


ui = r cos (kz3 + 8) — z1 = 


= r (cos ka3 cos 3 — sinkax3 sin 3) — z1 (8.162) 
ug = r sin (8 + kz3) — x2 
= r (sin kz cos 8 + cos kx3 sin 8) — z2 (8.163) 


Next, under the assumption that J is small such that 


sin kx3 ~ kz3 


coskx3 œ 1 
we finally get the expressions of the displacements 


uy = —kx3x2 (8.164) 
u2 = kx3x1 (8.165) 


Now, making use of compatibility and constitutive equations 
governing the linear static problem, see section 4.2, we can obtain 
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the strain and the stress tensors as follows 


0 0 —kxo 


1 
Eij = z 0 0 kzı (8.166) 
—kx2 kai 0 
E 0 0 —kxo 
Tij ya a AY, 0 kai (8.167) 
2 (1 tr v) —kxo kay 0 


In order to verify whether the state of stress stemming from 
the assumption on the displacement field is consistent with the hy- 
potheses on which the Saint Venant’s model is founded, we want 
to be sure that on the lateral surface the external forces vanish. To 
verify this consider the unit normal vector n = cos Be, + sin Geo, 
and hence the stress state is 


1 
031n + o32N2 = -3% (x2 cos 3 — xı sin p) = 0 (8.168) 
JESN 
=0 
The above is the proof the solution is right. 
On the other hand, the boundary condition at x3 = l, where 
n = €3, requires that 
053N3 = fj (8.169) 
thus, f3 = 0 and by virtue of equations (8.39) and (8.40), the only 


non vanishing component which produces the above state of stress 
is 


M3 = f (foo = fixa) dA = I (02301 — 01312) dA 
A A 


kE 2 2 kE 
= —_—- dA = Jo = kd, 8.170 
where u = EST see table 4.1 on page 95, and Jo is the polar 


moment of inertia for a circular cross-section. 

Usually in practical applications the problem presents an inverse 
formulation, namely, the unknown is the state of the stress and the 
given datum is the external couple M3, so we can easily derive 

E M M3 

r2= 
21+ v) po Ip 

E Ms M 
023 = z1 = LI 8.172 

2(14+ Vv) uJo Jo ( ) 


013 = T2 (8.171) 
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8.5.2 Cylindrical bar 


The hypothesis of plane cross-sections, i.e. u3 (71,2) = 0 is 
solely allowed for circular cylinder. It can be proved, in fact, that a 
generic-shaped cross-section under a torsional couple warps. This 
can be seen, for example, looking at equation (8.168) which would 
not be satisfied if the unit normal vector were not given with respect 
to a circular cylinder. 

Therefore, to remove the plane sections hypothesis we shall as- 
sume the following displacement field 


Uy = —kx3%x2 (8.173) 
u2 = kx3x1 (8.174) 
uz = kg (£1, £2) (8.175) 


where is an unknown function that must be determined in order 
to satisfy all the required conditions. 
The strain and stress tensors become 


k 0 0 PI — T2 
Eiji = 3 0 0 p 2 + Tı (8.176) 
2 
p1ı— T2 LP2+21 0 
0 0 PI — T2 
Tij = ku 0 0 p 2 + Tı (8.177) 


PITT LP2+21 0 


The equilibrium condition on V leads to 


013,3 = 0 (8.178) 
023,3 = 0 (8.179) 
031,1 + 032,2 = 0 (8.180) 
hence 
013 = 013 (£1, £2) (8.181) 
023 = 023 (£1, £2) (8.182) 
P11 + ¥,22 = 0 (8.183) 


Moreover, the boundary condition on the lateral surface imposes 


031M + 032n2 = 0 (8.184) 
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and considering the tensor in (8.177), the above condition becomes 


Pini — Tan] + g,2N2+ TIN? = 0> (8.185) 
Vo -n = ron, — tn (8.186) 


So, let A and 0A be the cross-section of the beam and its bound- 
ary, respectively, the torsional problem of a cylindrical beam can be 
stated as follows 


{ V?p=0 VpEA 


Vo- n= rn) — rino VpedA G18) 


where the first equation stems from (8.183) and the boundary con- 
dition is provided by (8.186)?. 

The function y = ọ (z1, 22) is named the torsion function. 

The problem (8.187) is known as Neumann’s problem and con- 
sists in determining a function which is harmonic in a given region 
and whose normal derivative is prescribed on the boundary of the 
region. 

Here we will not give the entire analytical solution for Neu- 
mann’s problem, but we shall just give some general statements. 
The whole problem is solved in [1]. 


Stress function 


Since ọ (z1, £2) is harmonic on A it is possible to construct 
the analytic function y + iw of complex variable zı + îx2, where 
w (x1, 22) is the conjugate harmonic function linked to ~ (#1, x2) 
through the following Cauchy—Riemann equations 


G1 =Y2 (8.188) 
p2 = -Y1 (8.189) 


The theoretical background of the above statements is beyond 
the scope of this book, anyhow the reader can find a comprehensive 
formulation of the torsion problem in [1], [2] and [6]. 


5Note that Vy is the gradient of the scalar field y and the scalar product 
with the unit vector n gives its normal derivative as follows 
dp 


SE = Vp.n=gradp.n= e n= pam 
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Suppose that the boundary of the cross-section A is described 
by a curve c : IR > IR? such as st (z1 (s) , x2 (s)). We can find 
the tangent vector t to the curve as t = &;&;, where ti = di and 


i = 1,2. Hence, the conditions n-n = 1 and ñ -t = 0 allow us to 


compute the component of n as 


t2 


POSE eee 8.190 
JEF E 
pre (8.191) 


VET + #3 
Replacing the above expressions into the second equation of 


(8.187) and by expanding the gradient of the torsion function we 
have 


(0121 + 0,282) ñ= T2 vere: +2} vere: (8.192) 
that is 
ty iy 
pam PO ara Vata 
pits — Poti = x2%2+ 2121 (8.193) 


and now, making use of equations (8.188) and (8.189), the latter 
becomes 


W2h2 + dadi = Leto +1121 (8.194) 
then d ti 
A a T (xj + 23) (8.195) 
so that we finally obtain the expression of the function w as 
1 
p (£1, £2) = 5 (aj + x2) + const. (8.196) 


The arbitrary integration constant does not affect the final so- 
lution in terms of stresses and deformations, in fact two different 
constants will lead to two solutions which differ from one another 
only by a rigid motion. 

Moreover, from Cauchy-Riemann equations (8.188) and (8.189) 
it follows 


22 = P12 (8.197) 
Wil = 921 (8.198) 
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and by summing member by member it becomes 
V24=0 (8.199) 


Hence, Neumann’s problem (8.187) turns into the following Dirich- 
let’s problem 


{ V2 (21,72) = 0 VpeA 
Y 


T1, £2) = 5 (x? + 3) VpedA 200) 


The function Y% = Y% (x1, 22) is named the stress function. 
Suppose we are able to solve the problem (8.200). The stress 
tensor then becomes 


0 0 P2 — £2 
Oj = ku 0 0 — 1 + 2X1 (8.201) 
ba -t2 Vita 0 


thus, by virtue of equations (8.50), (8.56) and (8.41) and by con- 
sidering the boundary conditions (8.14) and (8.15), we can set 


M3 = ku f (— (da — z2) x2 + (1 + £1) 21) dA 
= ku f (xi + 23) = (dai + Y 2x2)) dA 


= ku (4. E [eae + 202) dA) (8.202) 


that solved for the elastic constant gives 


k= M3 
~ (Jo — fg (bara + 4,202) dA) 


Finally the state of stress in a generic cross section x3, given an 
external twisting action, is 


(8.203) 


3 

= BI 8.204 

031 DIA (Y2 — x2) ( ) 
M3 

= —— — 8.205 

032 a (Yı = 21) ( ) 


where we have defined 


= Jo — fa (rei + 4202) dA 


8.206 
7 (8.206) 
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8.5.3 State of strain 


Making use of equation (4.25) and table 4.1 we can write the 
strain tensor as follows 


1 0 0 013 
Eij 2G 0 0 023 = 
031 032 0 
M 0 0 (Y2 — x2) 
ST 0 0 (in) | = 
ad (a — x2) (#1 — 21) 0 
0 0 (pa — £2) 
M: ; 
= T. 0 0 (p2 +21) (8.207) 
evo \ (p1= 22) (p2 t2) 0 


8.5.4 Displacement field 


At the beginning of section 8.5.2 we already introduced the dis- 
placement field associated to the torsion problem. However, we im- 
plicitly assumed that the rotation axis for the point p € A is coin- 
cident with the axis x3. It is possible to show that this restriction 
does not affect the validity of the results in terms of stress. 

Anyhow, if we repeated an integration procedure similar to that 
we made to compute the displacement fields in the axial force and 
pure bending cases, we would find that here the displacement field 
assumes the following general form 


uy = —kz3 (£2 — £5) (8.208) 
uz = kzz (zı — xf) (8.209) 
uz = kg° (x1, £2) (8.210) 


where c = (x$, £$) is the point about which the rotation occurs and 
° is the torsion function relative to the rotation point c. 
The state of stress consistent with the above displacement com- 


ponents is 
0 0 Py r2 +5 
f= 5 0 0 Po +21 — ri 
pı —t2+% Potti- x 0 


(8.211) 
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and consequently it becomes 


0 0 py — t2+ 25 
Tij = Gk 0 0 Ph +%1- Zi 
gi T2: +13 Yo +r- dî 0 
(8.212) 


By using the above state of stress instead of that in (8.177) the 
equilibrium equation (8.180) assures the following condition 


V? =0 (8.213) 
Moreover, the boundary condition on the lateral surface implies 
(p5 — x2 + 5) ni + (ph +21 — rf) n =0 (8.214) 
which can be also written as 
(p5 + x3) ni + (ph — ef) no = Tani — cino (8.215) 
which is equivalent to the following expression 


(yo + 5x1 — TY L2) 1 nı + (yo + #521 — T{L2) » ng = Lon, — Linz 

(8.216) 

From the latter it is straightforward to realize that the new tor- 

sion function ~° = y°+ x52, — xf x2 must satisfy the same condition 

on OA that p must satisfy. In addition to that, condition (8.213) 

guarantees that the Laplacian of ° vanishes. Thus, Neumann’s 
problem assumes the following form 


l V? =0 VpeA 


VÈS- n= ron) — rino VpedA (8.217) 


Due to the uniqueness of Neumann’s problem the two torsion 
functions y and $° can only differ each other by a constat value, so 
that 


d=p+t (8.218) 


from which 
p? = p — x581 + LÍT +t (8.219) 


If we use the torsion function in (8.219) to compute the state of 
strain in (8.211), we will immediately find that the state of strain 


178 SAINT-VENANT'S PROBLEM 


remains unaltered compared to that obtained by using the torsion 
function y. That means that the kinematics related to y and ° 
only differ one another by a rigid body motion which does not alter 
the results in term of stress. 

By virtue of these remarks we can state that the hypothesis of 
assuming the torsional axis coincident with x3-axis was reasonable 
and acceptable. 

As the last point of this section we want to find the position of 
the point c called the center of twist obtained by the intersection 
of the axis of twist, that is the axis parallel to the generators of a 
cylinder undergoing torsion - located so that the displacement of 
any point on the axis is not affected by any rotation, and a generic 
cross section. 

From equations (8.208) to (8.210) it is possible to put zero the 
mean value of the displacement ug and the mean value of the rota- 
tions of a given point p € A by putting 


l p° (x1, £2) dA = 0 (8.220) 
A 
J p° (£1, £2) 2d A = 0 (8.221) 
A 
if p° (01,02) 21dA = 0 (8.222) 
A 
which, making use of equation (8.219), give 
t=0 (8.223) 
f e (£1, £2) rad A + zÍ Jz = 0 (8.224) 
A 
f y (z1, £2) zt1d A — x5Jy = 0 (8.225) 
A 
and finally 
1 
= -=f Y (£1, £2) ad A (8.226) 
Je JA 
1 
T= | p (£1, £2) 41dA (8.227) 
Jy JA 


8.5.5 Strain energy 


We can equivalently use three tools to compute the strain energy 
associated to a beam undergoing torsion: 
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Work done by the external forces. By virtue of the Clapey- 
ron’s theorem, see equation (6.15) on page 114, the strain 
energy is 


1 
= 5 M30! (8.228) 


where M} is the only external force and W is the rotation in 
the (x1,2)-plane at the point of application of the torque, 
i.e. at the end x3 = l. 


As already stated, the twist rotation 0 per unit length is given 
in accordance with equation (8.157) as follows 


ge g 


= = .22 
ra (8.229) 


thus, the whole twist rotation all along the beam is readily 
given by the integral 


Mi Mi 
d'= faa = | ga =$] 8.230 
l i 1 GoJo iù GoJo ( ) 
and finally 

a 2 

1Mi1 
= .231 
2 GoJo (A 


Work done by the internal stresses. By recalling equation 6.5 
on chapter 6, we can set 


1 
D = 5 [ousia (8.232) 
2 Jy 


which in this specific case becomes 


1 
= F (031€31 + 032€32) dA = 
y 


1 
= za | a+ dA = 


N 


M3 2 2 
a? - dV = 
20027? | (a x2) + (p2 + x1) ) 
Mil 2 2 
IGorJ? | (pi ite + 9 + 0201) dA + 0 


(8.233) 
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The above integral can be rewritten taking into account the 
following identities 


gi — vate = (p (91 — 22)) 1 — 011 (8.234) 
Pa + p20 =(P (9,2 + 21)) 3 — 0022 (8.235) 
hence 
i. (på — ite + p% +9201) dA = (8.236) 
| (vy — x2)), + (p (2.2 + zija) dA+ 
= f pV? pd A (8.237) 
A 


Now we can realize that the last integral contains the first 
condition of Neumann’s problem so that on the domain A it 
vanishes. Moreover, by using the divergence theorem for the 
first integral at the second member we obtain 


(pi — p122 + p% + 9201) dA = (8.238) 
A 
f, ((p (p1 — x2)) ni + (p (p,2 + £1))n2)ds = (8.239) 
f p(Vp-n— x12n1 + x1n2) ds (8.240) 
OA 


The latter integral includes the boundary condition of Neu- 
mann’s problem that is identically zero. 


Finally we have proved that 


/ (Pi — p12 + 94 + 9,201) dA=0 (8.241) 
A 


and so the strain energy can be expressed as follows 
M31 
= — > 5 oJ, 8.242 
where since M3 = M3 we obtain 
_ Msi 
B 2GoJo 


that is the same result we found through the previous method. 


(8.243) 
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Work done by the internal forces. Note that here the sign con- 
vention is taken in accordance with that assumed in section 
8.1.1, so for a generic cross section we have 


1 
db = 5 M3d0 (8.244) 


and, by integrating along the entire beam, we obtain 


1 1 M2 
© = 5 | Mabars = sl 
l 


8.245 
2 GoJo ( ) 


8.5.6 Torsion of tubular beams: Bredt’s theory 


To solve the problem of tubular beams under torsional couples 
we can make use of an approximate theory that requires just the 
equilibrium equations. Consider a generic domain A and a closed 
curve c within that domain. See figure 8.11. 


Figure 8.11: The sub-domain A, bounded by the curve c. 


Let 7 be the tangential stress vector lying in the domain A so 
that 
T = 031€1 + 032€2 (8.246) 


The equilibrium condition (8.9) leads to the following alterna- 
tive expression 
Tii = 0 i=1,2 (8.247) 


that is nothing more than div 7 = 0. 
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Now let A. C A be the area included by the curve c, invoking 
the divergence theorem (1.143), we have 


f divrade= f(7-n)de=0 i=1,2 (8.248) 


Equation (8.248) proves that given a generic region A, the stress 
flux through its boundary c always vanishes. 
Consider now two curves cj and co as shown in figure 8.12. 


Figure 8.12: Stress flux within a small region included by two closed 
curves and two generic transversal sections. 


Due to the result in (8.248), no stress flux passes through the 
arches a — c and b — d, so for the closed area abcd the flux balance 
is given as follows 


-f 7 taads + f T- teads = 0 > 
Sab Scd 


= | Tab (8) ds + ni Ted (8) ds = 0 (8.249) 


Sab Scd 


where ft is the unit vector normal to the transversal sections a — b 
and c — d, respectively, while Sab and sp are the the lengths of the 
transverse sections, i.e. the thickness of the tubular section. 
Assuming that s is sufficiently thin, we can consider the average 
value 7”° instead of T = 7 (s), so the above integral can turn into 


Tab Sab = ed Sed (8.250) 
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where 
Tab = — | T-tards (8.251) 


Finally, since the sections Sab and Spe have been arbitrarily cho- 
sen along the all tubular section, equation (8.250) provides the fol- 
lowing result 

T” s = constant (8.252) 


and moreover, if we suppose that the flux lines are parallel to the 
midline c, i.e 7 = |7| = 7™, we have that the tangential resultant 
for unit length is given as follows 


dF, = Tsdc (8.253) 


Now suppose a generic equilibrium direction is fixed by the angle 
a as showed in figure 8.13, 


Figure 8.13: Stress resultants. 


the resultant force acting on the cross-section has the following 
expression 


girs cos adc = Ts f cos adc = 0 (8.254) 
c c 


Accordingly, the in-plane rotational equilibrium is satisfied by 
imposing 


M3 = freh (s) de = rs $h (s) dc = 2Ts Aec (8.255) 


where o is a generic point with respect to which we compute the 
moments, A, is the area included by the midline c. 
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Hence, the inverse problem is solved as follows 


M3 


= 8.256 
T 25Ac ( ) 
8.6 Bending and shear 
8.6.1 External forces 
The external forces acting on the base x3 = l are 
T 40, T40 (8.257) 
Ñ! = M = M = Mi-o (8.258) 
hence, the equilibrium condition imposes 
f! a fldA #0 (8.259) 
A 
T} sii fidA 40 (8.260) 
A 
N! | fidA =0 (8.261) 
A 
so that 
ff #0 (8.262) 
fi £0 (8.263) 
fA=0 (8.264) 
and the following condition has to be satisfied 
Mi = f (- fira Es far) dA=0 (8.265) 
A 


On the other hand, the rigid body equilibrium requires that on 
the base x3 = 0 the external forces acting are 


Pero, 0 (8.266) 
M? =, MS =—-Tl (8.267) 
N°=N'=0 (8.268) 

M? = -M =O (8.269) 
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Consequently the boundary conditions are 


Î #0 (8.270) 
AFA (8.271) 
Î$#0 (8.272) 


8.6.2 State of normal stress 


Let us start from a generic form of 033 as stated in equation 
(8.32). The condition 033 = 0 at x3 = l allows us to reduce the six 
unknowns to three, in fact we have 


033,c3=1 = a + bx + cxg + (d+ ex, + fz2)l = 0 (8.273) 
033|e3-0 = a + br + cry = f? (8.274) 


and the following expression could be a solution 
033 = (a + Bay + yx2) (l — #3) (8.275) 


where a, 3, y are the unknown constants. 
To compute the above constants we shall impose the equilibrium 
condition on the base x3 = 0, so we have 


NO = -f 0c33dA = -f (a+ Bx, + yz2)ldA = 0 (8.276) 
A A 


M? = -f 033r2dA = -f (a + Bx1 + yz2) lzr2d A = Til 
A A 
(8.277) 


MI = -f 033r1dA = -f (a + Bx, + 7x2) led A = Til 
A A 
(8.278) 


where we have made use of equation (8.275). 
Equations (8.276), (8.277), (8.278) represent a linear system in 
the unknowns a, 8, y, that is 


alA = 0 
-bh2 -7h = Ë (8.279) 
-bh -y2 = —T! 
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In the same way as made for the pure flexure case we may write 
the above system with respect to the axes of inertia (€,7), so it 


becomes 
0 =J yì T 
Je 0 BI \ = 


£ 
SERES 8.280 
B ( ) 
T 
n 8.281 
y J, ( ) 


and accordingly the stress produced by the couples is 


O 
033 = P (l — x3) (8.282) 


Taking into account that the bending moment due to the exter- 
nal forces Tz and T} propagates along the beam as 


Me=-Ti(1- 23) (8.283) 
M, = Tf (l — 23) (8.284) 


then, the state of stress normal to a generic cross-section assumes 
the following expression 


Me Mn 
033 = ae (8.285) 
Je In 


where we want to remak that the above expression is similar to 
equation (8.106), but here Mg and M, are not constant, they depend 
on the cross-section position, i.e. x3. 

A simpler solution is obtained by introducing the coordinate 
system (x,y,z) where x is the neutral axis n — n, y is the flexural 
axis f — f 

My T! (1-2) 


= 8.286 
y z ” ( ) 
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Figure 8.14: Shear regions. 


8.6.3 State of tangential stress: Jourawski’s theory 


This section is devoted to an approximate shear theory widely 
applied in practical cases, it is named Jourawski’s shear theory. 

Consider a generic cross-section of the beam and suppose to 
split the area A in two regions A; and Ag. See figure 8.14. 

We call I the line that divides the section and r the line normal 
to l. So it is possible to define a local coordinate system assuming 
{l, r} as Cartesian axes. Hence, | and 7 form a two-dimensional basis 
for the system. Let us consider now a three-dimensional portion of 
the solid included by two surfaces normal to the x3 axis, at x3 and 
x3 + dx3, respectively, and the plane 77. See figure 8.15. 

Let 73 be the tangential stress vector lying in the domain A, so 
that 

T = 031€] + 032€2 (8.287) 

The stress flux 73, passing through the line l is given by the 


scalar product 73-7, so that the equilibrium condition of the portion 


V is 
= / 033dA, + i. (033 + 033,3) dA, — f Tr3dl = 0 (8.288) 
Al Al 


lab 
then 
| 033,3dA1 = I Tradl (8.289) 
Al lab 
By using the result in (8.286), the latter equation becomes 
TI 
| ydA, = / Tradl (8.290) 
Jy Al lab 
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Figure 8.15: Beam splitting. 


Invoking the mean value theorem, we can set the mean stress 
T;3 along the chord lap as follows 


T3 = 7—- | Tradl (8.291) 


hence we can state that 


T, 
— ydAi = labT,3 > 
x J Ai 
Ty S41 
m yh le 
= 8.292 
Tr3 Jalat ( ) 


where Ty = fi is the shear force acting along the flexural axis; Six 
is the static moment of the area A; with respect to the neutral axis; 
J, is the entire cross-section moment of inertia with respect to the 
neutral axis; lap is the length of the chord. 

Equation (8.292) allows us to compute the mean value of the 
shear stress acting in the r direction normal to a generic chord 
which splits the section in two portions. Jourawski’s theory does not 
depend on the chord position, it is just required that it separates 
the cross-section in two parts. Moreover, the chord lap can be a 
polygonal line and in the case of tubular section it can cut the 
section more than once. 

The practical application of Jouwraski’s theory is allowed when 
the chord length is sufficiently small, so under this condition we 
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can approximate the mean stress with the actual one without loss 
of accuracy 
T3 ~ Tr3 (8.293) 


Furthermore, equation (8.5) on page 145 ensures the no-stress 
condition along the normal to the lateral surface of the cylinder. So 
in points a and b, see figure 8.14, the stress 7,3 must be tangent 
to the boundary lines of the cross-section. Considering the above 
condition on the smallness of lab, if the boundaries and the chord 
are orthogonal we can write that 7 = 73,7, i.e. there are no other 
components of the shear stress vector except the one along the r- 
axis. More details will be given later on symmetrical sections. 

As a concluding remark we want to show that the shear stress 
73, does not depend on which portion of the section we choose. In 
fact, if we consider the flux towards the area Ag we have 


moreover, we know that 
Sig + Sox = 0 => Siz = — Sor (8.295) 


so the shear stress equals 


zm = TyS2x 
ee diglab 
nd TS 
ge e vaa (8.296) 


a Ielab Jalan 


8.6.4 Tangential stress for symmetrical cross—sections 


Consider now a symmetrical cross-section under a shear force 
passing along the axis of symmetry that coincides with the flexural 
axis. See figure 8.16. 

Suppose that the section width is sufficiently small to consider 
valid Jourawski’s theory, then the shear stress along the chord lab 
is given by the following expression 


Ds 
Jzlab 


On the left side of figure 8.16 is showed the distribution of the 
static moment related to the portion A, computed with respect to 


Tay = (8.297) 


190 SAINT-VENANT'S PROBLEM 


8 


Figure 8.16: Symmetrical cross-section. 


the x-axis, namely 


Sin (1) = I, vil I "aiid (8.298) 
1 y 


Often in the practical application it is required to compute the 
maximum shear stress, so it easy to observe that since 


Ty = Tay (y) (8.299) 


then the maximum value is found by imposing the following condi- 
tion 


T3y,y = 0 


which implies 


(+) Ria (8.300) 
y 


IJ), tdy Pay 


Moreover, equation (8.298) assures that dS), = lydy, therefore 
the latter equation becomes 


y- £ =0 (8.301) 
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that is 
y= © — (8.302) 


This result implies that along the x-axis, i.e. when y = 0, the 
chord / has either an extreme value or is constant and that is suffi- 
cient to assure the condition expressed by equation (8.300). In fact 
we have 


dl 
Ce ify= 
7 0, ify=0 (8.303) 


therefore ni 
yP lax 
(T3y)max = dele 


(8.304) 


where lo denotes the length of the chord at y = 0. 

Now we can split the moment of inertia into the sum J; = Jiz + 
Joy = Sxh1+ Sho, where hı and hg are the distances between the 
centers of area C and C2 of the two portions separated by means 
of the chord lo and the center of the whole section G, respectively. 
Hence, if hy + ho = ho, we can write 


Je = Sxho (8.305) 


which leads us to write 


(T3y) max = T3y (0) = —— (8.306) 


| 


Figure 8.17: Maximum shear stress for symmetrical cross-section. 
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As showed in figure 8.17 the centers Cı and Cə coincide with 
the points of applications of the resultant forces produced by the 
normal state of stress 033 = 033 (y). 

As a concluding remark of this section we want to show that by 
using the equilibrium equation in (8.9), which has not been used 
yet, it is possible to know the distribution of the stress along the 
x-axis. To this end let us write equation (8.9) with respect to the 
neutral and symmetry axes, respectively 


T3x,0 + T3y,y + 933,3 = 0 (8.307) 
and the derivative with respect of x allows us to write 
T3r,ca = 0 (8.308) 


because equations (8.297) and (8.286) tell us that the first two terms 
of equation (8.307) vanish, so that 


Tg, = ax +0 (8.309) 


where a and #8 are two integration constants that must be found by 
means of the boundary conditions. The stress boundary conditions 
are known due to equation (8.5) that ensures the tangency condition 
of the shear stress vector 73 to the boundary of the cross-section. 
See figure 8.18. 

It should also be noted, as made clear in figure 8.18, that for 
any point on the chord | the shear stress vector is always lying on 
the line towards the point O that belongs to the symmetrical axis 
and is determined by the intersection of two tangent lines passing 
through a and b. 

Hence it is very easy to prove that the shear stress component: 
along the neutral axis x is given by 


Tag = — Tgy® (8.310) 


l 
T3x (2) = —T3y tana (8.311) 


I 
di (-2) = Ta tana (8.312) 
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O 


Figure 8.18: 03, distribution for symmetrical cross-sections. 


Finally, for a generic point within a symmetric cross-section we 
know the whole state of shear stress 


TS, 


= 2 8.313 
T3y Jal ( ) 
2t 
T3x = — Tye (8.314) 


where / is a generic chord that splits the section. 


8.6.5 State of strain 


With the same approach followed for the state of stress, the 
strained configuration of a beam under terminal forces, which pro- 
duce shear and bending forces along the whole beam, will be in- 
vestigated separately. Namely, by using the superposition principle, 
the deformation concerning a generic cross-section will be obtained 
as sum of the contribution due to the bending state of strain and 
the contribution due to the shear state of strain. 
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It must be observed that in this section we shall investigate 
only the global strain of the cross-section assumed remaining plane, 
instead of the local strain tensor Ejj. 


Bending strain 


Likewise the pure bending case analyzed in section 8.4.3 for a 
generic cross-section we have 


Mz 
EJ, 


Pr = KZ = z (8.315) 
where Yx is the rotation in the (y,z)-plane at the point of appli- 
cation of the internal couple M, = -ÎI (l — z) that is the bending 
moment produced by the external force T. The key difference with 
respect to the pure bending case is that here the rotation Yx is no 
longer constant, but varies linearly with z. 


2 


i da 
pi Ju: 


Figure 8.19: Bending strain for an infinitesimal beam segment. 


y 


Let us define now the rotation per unit of length as follows 


dgr Ma Îl- 2) 
dz EJ; EJ 


(8.316) 


Shear strain 


The natural consequence of the shear state of stress discussed 
before is the shearing strain that causes a sliding of the cross- 
sections that, initially plane, become warped. See figure 8.20. 

We shall focus our attention only on the sliding of the cross- 
section in order do describe its global deformation. Furthermore, 
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dz 
-T, T, 
A 


| 


Figure 8.20: Shear strain for an infinitesimal beam segment. 


we shall also show that the contribute of the sliding to the deflec- 
tion line defined as v = v(2)|x=y=0 will be sometimes neglected, 
sometimes not, depending on the geometrical features of the cross- 
section. 

With respect to figure 8.20 we can set 


dn = ydz (8.317) 


where, chosen two cross-sections dz apart from each other, dn rep- 
resents the strain due to the shear force. We can easily write the 
shear strain energy ®, as follows 


1 
®, = Fi (2TzyEzy + 2TzxE€zx) dV (8.318) 
y 


so that for a small beam’s portion dz, considering the constitutive 
law, the above energy becomes 


_ dz 


Dre 
Seco 


(72, + 72,) dA (8.319) 
which, recalling equations (8.313) and (8.314), turns into 


_ dz 2 Atan?a s 
d®, = 2G Ta (1 + — 7 ) dA (8.320) 


The shear stress Tz, depends only on y, so that the above integral 


ÉNotice that the subscript s denotes the portion of the energy associated to 
the shear force alone. 
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can be rewritten as 


dz [* 1/2 4tan? a 
aE | Ta ay f (1+ ) de = 
2G Jys i —1/2 1? 


dz M, ( 4tan?a N 
= | ndy|(x+ x = 
2G i y 312 -1/2 
T?dz f" S2 tan? a 
_ _ y x 
>; [ (1+ : ) (8.321) 


Finally, Clapeyron’s theorem allows us to equilibrate the strain 
energy computed by means of the internal stresses with half of the 
work done by the external forces 


1 
d®, = 5T,dn (8.322) 


so that equation (8.322) equals equation (8.321) as follows 


1 T?dz pu 92 tan? 
d®, = —Tydy = yY / x (14 an 2) da 


2 2GJ? ga 3 
Tydz f” S2 tan? a 
= Z (1 .32 
dn i. 7 (1+ 3 ) ay (8.323) 


and finally the relevant result is that the deflection of an infinites- 
imal portion of beam due only to the shear force is given by the 
following expression 


XyTy 
= .324 
dn GA dz (8.324) 


where we have defined the shear factor y1 as follows 


A fa S82 tan? a 
“=F f : (1+ 7 ) ay (8.325) 


Y2 


Equation (8.324) leads in the end to write the sliding angle y as 


IT 


(8.326) 
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8.6.6 Total strain energy 


Now we are able to compute the total strain energy for a beam 
subjected to shear and bending forces. Keeping unaltered the hy- 
potheses of symmetrical cross-section with respect to the y-axis as 
symmetry line and z-axis as neutral axis, we can state that the 
total strain energy is given by two terms 


= ð, +8, (8.327) 


where ®y is the strain energy concerning the bending state of strain 
and ®, the energy related to the shear state of strain. See figures 
8.21(a) and 8.21(b), respectively. 


T; 


| z 
i» 
n, 


(a) Bending state of strain. 


Vy 


(b) Shear state of strain. 


Figure 8.21: Two contributions to the state of strain for a beam 
subjected to terminal forces. 


Clapeyron’s theorem allows us to write easily both the contri- 
butions as 


®, = Tym (8.328) 


®, = 31yMs (8.329) 
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where the first contribution is easily obtained replacing the con- 
stitutive relationship for the normal stress, see equation (8.19) on 
page 147, and making use of equation (8.286) such as 


T2 (l — 2)? 
by = 5 | oe w= Ji =a y dV = 
v 2 Jy 


2 E ER 
1 Ti 2 2 dale 
= sez | aA [0-2 d = TT (8.330) 
__ 


=Jx 

To compute the shear strain energy we recall the expression 
obtained in the preceding section, see equation (8.324), and we use 
Clapeyron’s theorem 


1 Ve bd) 
è, = FT = 5% [dow = 5 no% fac = 3001 (8.331) 


In the end, the strain energy for the linear elastic beam with a 
symmetrical cross-section subjected to forces at its ends is 


273 2 
Tyl 1 XyTy! 
6EJ, 2 GA 


and consequently the total deflection at the point of application of 
the external force, in the direction of the force itself, is 


TAP, DL) 


b= (8.332) 


Tın =T, = 29 = 8.333 
yn y (m + s) 3EJ, GA ( ) 
hence 5 
Tyl XyTyl 
= = 8.334 
7= +s = 351 + GA ( ) 


8.6.7 Rectangular cross-section 


Consider a rectangular cross-section A = wh where w is the 
width and h is the height. The cross-section area is assumed to be 
constant, so that we can easily compute the moment of inertia and 


the static moment with respect to axes x — —y assumed as above 
to be the neutral and flexural axes, respectively. 
wh? 
Jpm 8.335 
p= (8.335) 


Sy == (x —y ’) (8.336) 
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Moreover, the tangential stress 73, = 0, so that in this case the 
stress state is completely defined by 


T,,(l-z 
oie = pea (8.337) 
_ Sg 8 Ty y? 
Met, (1 1° (8.338) 


The shear factor x, can be directly computed by using the ex- 
pression (8.325) which yields yy = £, 
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Appendix A 
Applications of the shell theory 


This appendix contains some applications of the shell theory discussed 
at the ends of the first four chapters. For all cases presented the external 
load ensures a membrane state of stress and consequently some analytical 


closed-form solutions can be reached. 


A.1 Spherical dome 
A.1.1 Geometry 


The spherical dome is a shell modeled on a portion of sphere 
having radius r and aperture 7/2 (hemisphere). Given the geometry, 
the first step is to identify the simplest coordinate system able to 
describe such a geometry. Of course it is a spherical system, see 
section 1.4.3 on page 19. 

Let X be the spherical coordinate system! so that 


X = (0,9,0):E— R? (A.1) 


where E is the affine Euclidean space in which the surface Q is 
embedded. The origin of the system is located at the center of the 
hemisphere. With respect to a Cartesian coordinate system, the 
following transformations hold 


x = psinpsin® (A.2) 
y = psinycosv (A.3) 
Zz = pcos p (A.4) 


The adapted coordinate system X induces the surface coordi- 
nate system Xİ by imposing the constraint p = r. Therefore, the 


Note that this coordinate system has been slightly changed compared with 
that depicted in figure 1.3. 


Claudio Borri, Michele Betti, Enzo Marino (a cura di), Lectures on Solid Mechanics, ISBN 978-88-8453- 
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induced coordinate system is 
Xt = (t, 9) : Q > R? (A.5) 


The covariant and contravariant expressions of the metric tensor 
gi associated with the induced coordinate system are, respectively 


g= r°d? Qd? + r? sin? yd? @ qd? (A.6) 
1 = = 1 = _ 

pr ue A. 

g A ETA (A.7) 


The nonvanishing Christoffel symbols on Q are 


I's = — sin cos p 
T? eS ed COS L 
ee op sing 


The unit normal vector of Q is 
i= Os (A.8) 


The Weingarten tensor and the second fundamental form for Q 
are, respectively 


L= 1 (a @ 9, +4” ® dy) (A.9) 
L= r(d’ ® d? + sin? yd” @ d?) (A.10) 


A.1.2 Displacements and strains 


To compute the in-plane state of stress only the stretching strain 
tensor a is required 


Oy = Vos tru® (A.11) 

agg = va, + sin Y cos y + r sin? pus (A.12) 
1 cos 

agg = = (Vp +09,6) — —V (A.13) 


2 sin Y i 
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A.1.3 Equilibrium and constitutive law 


The equilibrium equations (3.149) to (3.151) for a spherical shell 
assume the following form 


N?? p + cot pN¥? — sin pcos yN? +q? =0 (A.14) 
—N??r — Nr sin? p + g& =0 (A.15) 
N°? p +3cot yN”? + q° =0 (A.16) 


The constitutive equations are 
N? =D 1 E 
DA Yow +ru> ) + 
+D (Fiato + sin y cos pup + rsin? ov) (A.17) 


1 
NÙ =D (vow + sin Y cos pVy + rsin? gwt) + 


rísin y 
— s (vy, +70) A.18 
r4 sin? y ( we ( ) 
lv 1 COS Y 
lap = D| ———- (vy, so) — A.19 
È a? pe (Uy o +U9 g) eu) ( ) 


Load case: self weight 


The dead load due to the self weight provides, of course, a sym- 
metrical action so that the expected solution will not depend on 
ù. 

Suppose the load per unit area is g, uniformly distributed through- 
out the shell. The vector has only the vertical component 


g = -E (A.20) 


whereas, with respect to the basis {d,, 0g, n} the vector load g is 
written follows 


gf? = —qf cos pn + q? sin yd, (A.21) 
By multiplying equation (A.14) by sin y we obtain 


(sin pNP), — sin? y cos yN”” + sin yq? = 0 (A.22) 
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Let us introduce now the physical components of the stress ten- 
sor N, so that 


a N° abg J 
NSO? = agaj = NP el (A.23) 


Hence, equation (A.22) becomes 
(sin pN <P>) p — cos pNP? + rsinyq<?* = 0 (A.24) 


Analogously, by multiplying equation (A.16) by sin? y, consid- 
ering the physical components and noticing that q? = 0, we obtain 


(sin pN<*?*),,, + cos pN <? = 0 (A.25) 
The remaining equilibrium equation becomes 


N<99> N<vo> 
+q =0 (A.26) 


r r 


where, resolving equation (A.26) for N°, equation (A.24) turns 
into 


(sin? pN<?"*),, = (g°&rcosp— q<?"rsiny)sing (A.27) 


which can be integrated as follows 


sin? pN<??? = T r(g° (9) cos è — q? (Q) sin $) sin ddd + K 


Á (A.28) 
Equation (A.28) represents the equilibrium of a spherical cap 
included by latitude p and € [%, 7/2]. In particular the quantity 
2rrK, excepting the sign, equilibrates the resultant acting on the 
cap identified by the aperture ¢. 
Considering now equation (A.21) 


sin? pN<??> = — rg” | — cos ġ j (A.29) 


for the latitude y the whole meridian stress when 6 =0 => K =0 


1S 
z 


N<69> = rad cos £) D rq (A.30) 
sin Y 1+cosy 
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so that equation (A.26) becomes 


+2 
N<98> — ro? sin Y — COS p A31 
rg ( 1 + cosy ) ( ) 


The third equilibrium equation does not depend on the two 
latter results, therefore, since q? = 0, we have 


N> = 0 (A.32) 


Load case: uniform load on the horizontal projection 
of the shell 


This load case keeps unaltered the simplifications regarding the 
symmetry already discussed in the preceding case. Indeed, here too 
we are looking for a solution not depending on %®. 

The load q* is now projected on the horizontal plane 


q = —gq° cos pe, (A.33) 


therefore with respect to the local basis, the physical components 
are 
<> _ Zz 2 Z= 3 5 
q~° =—¢ cos* p+ g°nsin y cos pO, (A.34) 
By means of a procedure similar to that formerly used we obtain 
that equation (A.28) now becomes 


sin? pN<99> = iL CC) cos @ — qg°f?(@) sin d) sin @d@ò + K 


p 
p 
= | —rg’ sin y cos y + K (A.35) 
p 
from which 1 
sin? pNS9?> = — [cos? ls (A.36) 


Next, if = 0 > K = 0, the whole meridian stress is 
1 
NS<P> = a UCA (A.37) 
Finally, from equation (A.26) we obtain 


1 
N< — rid cos 2 (A.38) 
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A.2 Cylindrical shell 


In this example we want to compute the stress state for a cylin- 
drical shell subjected to some of the most typical load conditions, 
e.g. uniform pressure, dead weight, hydrostatic pressure. 


A.2.1 Geometry 


Obviously we choose as an adapted coordinate system a cylin- 
drical one with a little rearrangement compared with the one intro- 
duced in section 1.4.3 on page 18, 


X =(9,2,p): E > IR? (A.39) 


where, as usual, E is the affine Euclidean space in which the cylin- 
drical surface Q is embedded. The relationships between the Carte- 
sian system, with the origin along the axis of the cylinder, and the 
cylindrical coordinates are 


x = psenv (A.40) 
y = pcosì (A.41) 
Bae (A.42) 


The above adapted coordinate system induces the surface sys- 
tem XT due to the constraint p = r, where r is the radius of the 
cylinder. So we have 


Xt=(@,2):Q—- R? (A.43) 


The covariant and contravariant forms of the surface induced 
metric are, respectively 


g=r°d° 9 d° + d° @ d° (A.44) 
1_ > m = 
J = 7200 ® dg + Oz Q Oz (A.45) 


All Christoffel symbols vanish on Q. 
The unit normal vector of Q is 


n=O, (A.46) 
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The Weingarten tensor and the second fundamental form are, 
respectively 


L= 1408 dg (A.47) 
E 
L=rd’ 8d’ (A.48) 


A.2.2 Displacements and strains 


To compute the in-plane state of stress only the stretching strain 
tensor a is required 


Agg = Vy y + rue (A.49) 
1 

Adz = z (Ure +259 ) (A.50) 

Azz = Vzsz (A.51) 


A.2.3 Equilibrium and constitutive law 


For a cylindrical shell subjected to a membrane state of stress 
the equilibrium equations in the scalar form are 


N°? 9 +N™*,, +p" =0 (A.52) 

N”? + N”, +p = 0 (A.53) 

-N®? Log + pî = (A.54) 

N? = N° (A.55) 

The constitutive equations assume the following form 

D(1 

NI? = 3 (= (00. +rve) + vee | (A.56) 

di 

N = D ( oo (wo + 0) (A.57) 
v 

N” =D (30o +rve) + Vase ) (A.58) 


Load case: uniform pressure and self weight 


This load condition is characterized by two load components, 
namely gf and q7. The symmetry around the z-axis permits to 
delate all terms containing the derivatives with respect to ù. 
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The equilibrium equations become accordingly 


3 
nif (A.59) 

r 
N,,=0 (A.60) 
N” „+p =0 (A.61) 


Next, taking into account the boundary conditions (at z = 0) 
related to the particular load condition and using the physical com- 
ponents, we obtain 


N> — qr (A.62) 
N> = 0S NS =0 (A.63) 
VA 
N“ z+ =0 > N° = | -gdt + K > 
0 


N? = N = g (z —h) (A.64) 


Thus, the only nonzero components of 0 are those along € and z 
due to the self load and to the Poisson effect, which are respectively 


EL rg + rvqg? (z — h) 
5 E(2e) 


v = 355 ( yf (5 hz) vr) (A.66) 


Hydrostatic pressure and self weight 


(A.65) 


In this case the load vector g is made up of two components: qË 
and q* and the equilibrium equations are 


3 
wee E (A.67) 

r 
N (A.68) 
N”, ta = (A.69) 


Furthermore, by taking into account the boundary conditions 
(at z = 0) related to the particular load condition and using the 
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physical components, we obtain 


N<%> = Gr (A.70) 
NS (A.71) 


VA 
N” z +0 =0 > N° = f —q dt + K > 
0 
N” = N? =-q(z-h) (A.72) 


The nonzero components of v are 


hi (ŅJLr 
En LEA 
v E° ( De + v7) (A.73) 
z 2 
2_ tvy- d (2 
= A.74 
v Bae (5 hz) (A.74) 


A.3 Hyperboloid of one sheet 


The last example we propose concerns an hyperboloid of one 
sheet, that is the geometry of shell structures usually adopted for 
cooling towers. The structure is supposed to be loaded by the self 
weight lone so that the axial symmetry is preserved. 


A.3.1 Geometry 


The adapted coordinate system is X = (f,V,z). Here, as made 
for the preceding geometries, we will begin computing the metric 
tensors and the Christoffel symbols for such system. Then we will 
consider the surface Q, i.e. the hyperboloid, described by the in- 
duced coordinate system XÎ = (Vt, z') obtained by imposing the 
constraint fig = 0. For this system the metric and the Christoffel 
symbols will also be computed. 

The hyperbolic coordinate system is 


Cai IR? (A.75) 


with the origin in o € E that coincides with the origin of a Cartesian 
coordinate system. See figure A.1. 
We use the coordinate function f to define the surface Q, i.e. 
f=0, that is characterized by the following implicit expression 
ety 2 


po i= (A.76) 
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Figure A.1: Hyperbolic coordinate system. 


The function f is defined as 
f=p- p(z) (A.77) 


where p(z) = 2b? +2 = fa? + y2. 
Whit respect to the Cartesian system the following coordinate 
transformations hold 


z= (f + ZVI +2) cos (A.78) 


y=(f+V0 +2) sind (A.79) 
AS) (A.80) 


The covariant and contravariant expressions of the metric tensor 
are, respectively 


VIE 2) d? @ d'+ 


3 
b 


+ 1) d? 9 d? (A.81) 
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Si azz? = A 
g= (fa SE 2) + i) Of Q Ort 
eee ð; 9 ð; +; 9 ð 
o FT (ð; @ Oz + Oz @ Of) + 
1 
f+ avez) 


The nonzero Christoffel symbols for the adapted coordinate sys- 
tem are 


+ -s Q Dy + 3: @ O: (A.82) 


b 
pres: A.83 
(\ /b2 + 2°)? ( ) 
Py) = (+ VP +2) (A.84) 

1 
eresse = A.85 
fo of ft? ma 2 ( ) 

az 

T% = T9 = (A.86) 


b(f + EV b2 + 22) Vb? + 22 


Consider now the constraint p = ¢Vb? + 22, i.e. f = 0, in such 
a way we pass from the adapted coordinate system X = (f, 0, z) to 
the induced one Xt = (v1, 21)?. 

With respect to the surface coordinate system the expressions 
of the covariant and contravariant metric tensor are, respectively 


t= a E @ d+ a d” ® d° (A.87) 
IT O CCE) 000 

5 b? 1 = = b? (6° + 2?) = _ 
gi = peA 2) Og Q dg + Pee 0,90, (A.88) 


and the nonzero Christoffel symbols Tİ are 


ape = Ls i (6? is 22) (A.89) 
2 202 x 

r= a A. 
ZZ [a2 22 + b2 (b2 + z2)] (b2 + 2) ( 90) 
z a?z(b? + 22 


a?z? + b?(b? + 2?) 


?From now on we will omit the symbol + to denote the entities on Q when 
it is not ambiguous. 
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Next, the unit normal vector is 


ð, (A.92) 


RE az? + b2(62 + 22) = az 
j b? (b? + 2?) : fa2z? + b2(b2 + 22) 


and the Weingarten tensor and the second fundamental form are, 
respectively 


b2 


L= d’ ® ĝy+ 
aya?z2 + b2 (b2 + 22) : 
4 
do -d 98, (A.93) 
(a2z? + b? (b2 + 22))2 
be e ia 
© az +b (b2 + 22) 7 
2 
ay d? ® d° (A.94) 


(b2 +2?) fa2z? + b2 (62 +22). ~ 


The total curvature of the surface, i.e. the Gauss curvature, and 
the mean curvature are, respectively 


p6 
[a2z2 + b2 (b2 + 22)]? 
2021,2 _ 12 A2 ,2 
ya DE b°) + b*(b te) (A.96) 
aļa?z? + b?(b? + 22)]? 


(A.95) 


Moreover, from the Weingarten tensor, the principal curvatures 
are readily obtained 


2(12 2 
L= dt) (A.97) 
a(b? + 22) /a222 + 62 (62 + 22) 
4 
ae a (A.98) 


[a2z2 + b2(b2 + 22)]? 


The surface coordinate system X! = (V, z) is definitely comfort- 
able to describe and identify points on the hyperboloid, however to 
solve the in-plane equilibrium problem for the symmetrical load 
condition it is more convenient to chose instead of the z coordinate, 
the angle y that the segment line along © forms with the verti- 
cal axis z. This new variable is related to the former one by the 
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following relationship 


See a (A.99) 
OOS NPE RE | 


where k is a dimensionless geometric factor 
(A.100) 


A.3.2 Equilibrium 


By making use of the coordinate y, the solution of the equi- 
librium equations allows us to write the expressions of the stress 
tensor as 


N<e> = qav\/k? sin? y — 1 


— A.101 
<99> avk? — 1 NS99> (k? sin? p — 1)? 
N = Pron =| qeosg + 5 
p avk? — 1 
(A.102) 


where q is the dead load per unit of area (assumed to be constant 
along the thickness) while the function ¢ equals 


— cos Y 1 Vk? —1—kcosy 
C= +2 + In 
2(k?sin?y—-1) 4kvk2-1 Vk2—1+kcosg 
(A.103) 
Further details on the derivation of the above expressions are 
available in [17]. 
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